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ABSTRACT 

This  thesis  presents  the  theory  necessary  for  the  concep- 
tual design  of  a  tactical  missile.   The  design  process  begins 
with  the  well  known  linear  aerodynamic  theory  for  initial  sizing 
and  later  includes  nonlinear  effects  to  determine  the  final 
design  of  the  missile.   Where  theory  does  not  apply,  empirical 
methods  are  presented  which  are  known  to  give  accurate  results. 
An  air-to-air  missile  is  designed  for  a  specific  threat  as 
an  example  which  immediately  follows  the  development  of  the 
theory  for  each  section.   Several  small  digital  computer  pro- 
grams are  presented  and  used  for  analysis  of  specific  areas 
of  the  design.   One  large  program  (AEROl)  is  used  for  deter- 
mining the  aerodynamic  coefficients  of  the  final  design. 
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I.   INTRODUCTION 

This  study  was  made  to  present  a  method  for  the  conceptual 
design  of  tactical  missiles.   The  starting  point  for  the  design 
was  a  recently  completed  report  by  General  Dynamics,  Convair 
Division  entitled  Rapid  Approach  to  Missile  Synthesis  (RAMS) .  A 
procedure  was  then  developed  more  akin  to  aircraft  design  and 
which  bears  little  resemblance  to  RAMS  and  which  uses  basic 
equations  to  size  components  rather  than  nomograms  and  table 
look-up.   The  procedure  starts  with  a  threat  description  and 
proceeds  with  the  formulation  of  performance  objectives  and 
a  conceptual  design  of  a  tactical  missile  to  counter  the  threat 
The  design  is  not  unique,  but  as  will  be  shown,  is  a  compromise 
of  parameters  to  give  one   possible  solution  to  the  design  prob- 
lem; therefore,  the  point  design  arrived  at  is  not  necessarily 
the  optimum  design  for  the  presented  threat.   An  attempt  is 
made  to  find  the  optimum  performance  within  specific  areas  of 
the  design  process. 

Throughout  this  study  the  theory  involved  is  explained  and 
specific  examples  are  worked.   A  complete  design  example  is 
worked  out  in  detail.   It  is  an  air-to-air  missile  designed  to 
counter  the  new  Soviet  RAM-K  fighter  aircraft.   This  example 
is  worked  in  each  section  immediately  following  the  develop- 
ment of  the  theory  for  that  specific  area. 


18 


II.   PROBLEM  DEFINITION 

A.   THREAT  ANALYSIS 

1.   Operational  Requirements 

The  design  of  a  new  missile  is  usually  in  response  to 
an  operational  requirement  which  arises  as  the  result  of  one 
of  the  following:   (1)  A  new  technology  provides  the  means  to 
design  a  more  effective  missile  to  meet  a  current  threat.   An 
example  of  this  might  be  an  advance  in  material  science,  which 
allows  higher  inlet  turbine  temperatures  for  a  turbojet  engine; 
therefore,  allowing  higher  missile  flight  speeds.   (2)  Intelli- 
gence indicates  a  new  threat  for  which  existing  missiles  are 
not  effective.   (3)  Operational  reports  indicate  a  current  mis- 
sile is  inadequate  against  a  current  threat. 

Regardless  of  how  the  operational  requirement  is  derived, 
a  statement  of  the  threat  is  required  before  the  design  process 
can  proceed.   Experience  has  proven  that  one  missile  cannot  be 
designed  to  meet  all  types  of  threats  without  seriously  com- 
promising performance  or  effectiveness.   This  can  be  illustrated 
with  the  design  of  the  warhead.   A  contact  fuze,  shaped  charge 
warhead  designed  to  penetrate  and  kill  hard  targets  such  as 
tanks,  would  not  be  effective  against  a  highly  maneuverable 
aircraft  for  which  the  expected  miss  distance  is  several  feet. 
For  this  reason  the  design  of  a  missile  must  start  with  a  de- 
tailed analysis  of  the  threat.   The  more  detailed  this  analysis 
is,  the  more  effective  the  final  design  can  be. 
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2 .  Design  Example  (Operational  Requirement) 

A  design  example  will  be  used  as  a  continuous  thread 
throughout  this  thesis  to  demonstrate  applications  of  the  theory. 
An  air-to-air  missile  will  be  designed  to  counter  the  new  class 
of  Soviet  fighter,  which  is  in  the  advanced  development  stage 
at  the  Ramenskoye  Experimentation  Center.   The  fighter,  as  des- 
cribed in  Aviation  Week  [1, 2]has  been  designated  the  RAM-K.   The 
RAM-K  is  a  twin  engine  fighter  with  variable  geometry  inlets  and 
swing  wings.   The  aircraft  bears  a  resemblance  to  both  the  F-14 
and  F-15.   It  is  expected  to  be  the  recipient  of  a  new  look 
down,  shoot  down  radar  and  the  40  km  range  AA-X-9  missile.   The 
following  unclassified  dimensions  and  performance  data  are  avail- 
able on  the  RAM-K: 

Wing  span  40  ft 

Overall  length  64  ft 

Gross  weight  60,000  lbs 

Maximum  speed  M=2.5 

Service  ceiling  60,000  ft 

Figure  (2-1)  is  a  drawing  of  the  RAM-K. 

3 .  Scenario 

The  scenario  within  which  the  missile  is  expected  to 
operate  should  also  be  described.   If  the  normal  mode  of  operation 
of  the  threat  is  not  known,  an  attempt  should  be  made  to  define 
the  most  demanding  scenario  that  can  be  expected.   For  a  defen- 
sive weapon  the  most  challenging  incoming  threat  will  normally 
be  a  head-on  encounter.   The  threat  profile  may  vary  from  a  high 
level  attack  with  a  terminal  dive  to  a  low  level  attack  with  a 
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Figure  (2-1).  RAM-K  Fighter, 
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terminal  pop-up  maneuver.   In  the  case  of  the  AS-6  (Kingfish) 
two  modes  of  attack  can  be  expected.   In  such  a  case  both  pro- 
files must  be  evaluated  to  determine  the  most  demanding  in 
terms  of  missile  performance  objectives.   For  an  offensive  sys- 
tem, such  as  an  air-to-air  missile  designed  to  intercept  and 
destroy  an  enemy  fighter  before  it  launches  its  weapons,  the 
scenarios  analyzed  should  include  all  possible  encounter 
geometries . 

4 .   Design  Example  (Scenario) 

The  scenario  for  the  above  threat  would  likely  be  an 
intercept  situation  in  defense  of  the  fleet  high  value  unit. 
The  scenario  is  taken  to  be  a  head-on  encounter  with  the  missile 
and  the  target  at  the  same  altitude.   Since  the  combat  specifi- 
cations of  an  aircraft  are  normally  given  at  10,000  ft.,  this  is 
taken  as  the  scenario  altitude. 

B.   HISTORICAL  SURVEY 

Missile  design  is  an  iterative  process,  and  the  first  time 
through  the  design  loop  many  assumptions  have  to  be  made  con- 
cerning component  sizes  and  weights.   One  method  of  approach 
at  the  early  stage  is  to  employ  historical  data  of  existing  mis- 
sile sizes  and  weights;  since  justifications  for  these  parameters 
were  made  duirng  their  design  processes.   An  example  of  the  use 
of  historical  data  in  determining  the  initial  missile  length 
can  be  made  with  the  length  to  diameter  ratio.   The  length  to 
diameter  ratios  of  existing  missiles  of  the  same  class  as  that 
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being  designed  are  collected,  and  an  average  is  computed.   The 
diameter  of  the  design  is  fixed  by  one  of  three  driving  factors 
(propulsion,  warhead,  or  guidance) .   From  the  average  length  to 
diameter  ratio  the  initial  missile  length  is  then  estimated. 
From  this  historical  data,  initial  choices  based  on  the  experi- 
ence of  others  can  be  made  for  many  of  the  missile  parameters. 
These  parameters  define  a  baseline  missile,  which  is  the  initial 
configuration  from  which  design  iterations  and  refinements  can 
be  made. 

Since  missiles  are  designed  for  specific  missions,  speci- 
fic parameters  such  as  length  and  diameter  are  of  little  value 
in  comparing  missiles.   Dimensionless  ratios  such  as  length  to 
diameter,  L/D,  ratio  and  aspect  ratio,  AR,  are  more  meaningful 
when  relating  missiles.   Some  parameters  which  are  useful  in 
defining  the  baseline  missile  are  listed  below: 

L/D  =  Length  to  diameter  ratio 

L  /D  =  Nose  length  to  diameter  ratio 

AR   =  Aspect  ratio  of  the  wing 

AR.  =  Aspect  ratio  of  the  tail 

W/S  =  Weight  to  lifting  surface  area  ratio 

V,   =  S^lVfS   ,-d   ,.)  Tail  volume  coefficient 
t     t  t   ref  ref 

W  /W  .  =  Gross  weight  to  warhead  weight  ratio 
G  wn 

The  tail  volume  coefficient,  V  ,  is  a  dimensionless  para- 
meter used  to  initially  size  the  tail.   For  a  tail  control  mis- 
sile, it  is  a  measure  of  the  relative  control  effectiveness  when 
comparing  missiles.   For  a  wing  control  missile  it  is  a  relative 
measure  of  stability. 
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A  complete  historical  survey  should  not  be  limited  to 
the  parameters  listed  here.   Any  dimensionless  parameter  which 
will  add  information  about  the  proposed  design  should  be  inclu- 
ded for  completeness. 

Table  2-1  is  an  example  of  a  collection  of  such  para- 
meters for  existing  air-to-air  missiles  EL. 31 .  in  this  table 
the  subscript,  c,  is  used  to  indicate  a  canard  control  surface. 

C.   LAUNCH  PLATFORMS  AND  PHYSICAL  CONSTRAINTS 

The  problem  definition  must  also  include  a  description  of 
the  intended  launch  platform  for  the  missile.   The  aircraft  or 
shipboard  launcher  from  which  the  missile  will  be  launched  will 
fix  many  design  features  of  the  missile.   For  instance  the  most 
important  consideration  in  the  problem  definition  phase  is  any 
physical  constraints  imposed  by  the  launcher.   Since  it  is  not 
normally  feasible,  economically  to  design  a  launcher  to  fit  the 
missile,  most  new  missiles  must  fit  existing  launchers.   In  the 
case  of  shipboard  launchers  there  will  be  a  maximum  length  and 
diameter  and  a  maximum  launch  weight  which  can  be  accomodated. 

For  the  case  of  an  air  launched  missile,  there  will  be  a 
maximum  weight,  and  the  dimensions  may  also  be  limited  due  to 
the  performance  requirements  of  the  aircraft. 

For  the  air-to-air  missile  design  example  of  this  study, 

the  launch  platforms  will  be  the  F-16  and  F-18A.   Figure  (2-2) , 

which  is  from  Interavia  [4, 5]shows  the  pylon  weight  limitations 

for  these  aircraft.   From  these  figures  it  can  be  seen  that  the 

maximum  launch  weight  for  this  design  is  limited  to  2500  pounds 

by  the  pylon  limitations  of  the  F-18A. 
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Load  Carrying  Capability  of  F-18A 


Load  Carrying  Capability  of  F-16 


250/9.0 
3500/5.5   250/9.0 
2200/5.5      2500/5.5 


Figure  (2-2) .   Weight  Limitations 
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D.   MISSION  PROFILES  AND  PERFORMANCE  OBJECTIVES 

The  mission  profile  of  a  missile  consists  of  dividing  the 
flight  into  fundamental  segments  which  consist  of  a  single  func- 
tion, such  as  boost  to  cruise  speed  and  altitude,  cruise  to 
the  target,  and  terminal  homing  phases.   The  mission  profile 
will  vary  from  missile  to  missile.   For  a  cruise  missile  it 
may  consist  of  a  series  of  pop-up  maneuvers  and  low  level  cruises, 
For  a  short  range  missile,  the  profile  may  be  entirely  terminal 
homing.   The  Mach  number  and  altitude  are  specified  at  the  begin- 
ning and  end  of  each  mission  segment,  as  well  as  the  range 
covered  by  each  segment.   The  range  covered  by  a  segment  can  be 
considered  in  one  of  two  ways.   If  the  segment  is  short,  such 
as  the  terminal  phase,  the  distance  along  the  intended  flight 
path  is  considered.   For  longer  range  missiles,  the  distance  over 
the  ground  is  of  importance.   The  mission  profile  must  be  defined 
during  the  problem  definition  phase  in  order  to  specify  missile 
performance  objectives. 

The  mission  profile  of  a  missile  normally  consists  of  a 
boost,  cruise  (mid-course)  and  terminal  phase.   The  boost  phase 
accelerates  the  missile  to  its  flight  speed.   This  acceleration 
may  be  large  for  a  surface-to-air  missile  which  must  be  acceler- 
ated from  rest  to  a  high  supersonic  speed  or  it  may  be  small  for 
an  air  launched  missile  which  has  the  speed  advantage  of  the  air- 
craft from  which  it  is  launched. 

The  cruise  segment,  or  mid-course  phase,  primarily  is  used 
to  deliver  the  missile  to  a  point  in  space  where  the  seeker  can 
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acquire  the  target.   The  range  and  speed  of  the  cruise  segment 
is  then  a  function  of  the  required  stand-off  distance  for  the 
target. 

The  terminal  phase  is  somewhat  more  difficult  to  analyze. 
If  the  target  maintains  constant  heading  and  velocity  the  flight 
path  of  the  missile  may  be  modeled  by  circular  arc  segments.   This 
method  will  give  approximate  values  of  range.   If  the  target 
maneuvers,  the  missile  must  follow  and  the  range  and  speed  require- 
ments become  more  complicated. 

In  determining  the  range  and  speed  requirements,  all  expected 
encounter  geometries  should  be  analyzed.   The  most  demanding 
encounter  will  then  fix  the  performance  objectives.   The  most 
demanding  speed  requirement,  in  terms  of  maintaining  a  minimum 
stand-off  distance,  will  normally  be  a  head-on  encounter.   Al- 
though the  required  missile  speed  and  range  are  determined  in 
this  section,  the  missile  velocity  may  be  varied  later  in  the 
design  process  due  to  guidance  considerations. 

1 .   Design  Example  (Mission  Profile) 

From  the  threat  defined  above,  the  ideal  situation  would 
be  to  obtain  a  fire  control  solution  and  launch  such  that  the 
minimum  separation  distance  between  the  launch  aircraft  and  the 
target  is  40  km.   This  can  be  accomplished  in  one  of  two  ways. 
The  launch  aircraft  can  fire  a  semi-active  homing  missile  at 
such  a  range  and  speed  that  intercept  occurs  before  the  minimum 
range  is  reached,  or  an  active  homing  missile  can  be  fired,  and 
once  missile  lock-on  is  achieved,  the  launch  aircraft  can  maneu- 
ver to  maintain  the  minimum  separation  distance. 
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The  active  radar  homing  missile  would  decrease  the 
launch  range  and  the  range  required  of  the  missile,  which  will 
lessen  the  constaints  on  the  launch  aircraft.   The  terminal  por- 
tion of  the  engagement  is  a  function  of  the  guidance  law  and  will 
be  determined  in  Chapter  3. 

Two  cases  are  investigated  to  determine  the  effect  on 
the  range  requirement  of  the  missile  when  a  minimum  separation 
distance  from  the  target  to  the  missile  of  40  km  is  maintained. 
The  first  case  is  a  semi-active  homing  missile,  for  which  the 
launch  aircraft  must  maintain  a  closing  course  until  intercept. 
The  second  case  is  an  active  homing  missile  which  has  a  lock-on 
range  of  10  km.   The  launch  aircraft  may  then  maneuver  to  main- 
tain a  separation  distance. 

a.   Case  1:   Semi-active  homing  missile 


L  M 


^T 


R.  =  Range  at  which  the  missile  is  launched 

M_  =  1.5  =  Launch  Mach  number 

M^  =  2.5  =  Missile  Cruise  Mach  number 

a   =  Speed  of  sound 

V  =  M  a  =  1.5a  =  Launch  speed 

Ll  Li 

V  =  Ma  =  2.5a  =  Target  speed 

Vw  =  M,.a  =  2.5a  =  Missile  speed 
MM 
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The  instantaneous  range  from  the  missile  to  the  target  is 
given  by,  R^ 

*MT  -  R0  =  (VM  +  VT}t  (1) 

The  instantaneous  range  from  launch  aircraft  to  the  target 
is  given  by,  R^ 

*LT  =  R0  "  (VL  +  V*  (2) 

If  the  target  does  not  maneuver,  intercept  will  occur  at 
tf,  when  R^  -  0 

0  =  R0  "  (VM  +  VT}tf 

R        R 
t   =  Q ,  =  -0 

f    (VM+VT)    5a 

If  the  launch  aircraft  is  at  the  minimum  separation  distance, 
R^_  =  40  km,  when  intercept  occurs. 


40  km  =  RQ  -  (VL  +  VT)tf 


substituting  for  V  ,  V   and  t_ 

Li         T  X. 


R 
40  km  =  R0  -  (1.5  +  2.5)a(^Q.) 


5a 


Solving  for  R 
0 


R   =  20  0  km  =  Launch  range 
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The  range  required  of  the  missile  is  then,  R^ 

R^  =  V  t   =  100  km  =  53.96  nmiles 

If  the  missile  speed  is  increased  to  M„„  =  3.0,  then 

M 

RQ  =  148. 15  k  and, 


*M 


=  80.81  km  =  43.61  nmiles. 


b.   Case  2:   Active  Homing  Missile 

The  lock-on  range  is  a  function  of  the  seeker  in  the 
missile  and  will  be  covered  later  in  this  thesis.   If  it  is 
assumed  that  the  launch  aircraft  must  maintain  its  course  until 
lock-on  occurs  at  a  range  of  Rj^r    the  problem  can  still  be  solved, 
The  geometry  is  the  same  as  in  Case  1.   Instead  of  following  a 
constant  course  until  intercept,  the  launch  aircraft  must  now 
only  maintain  a  closing  course  until  R^   =  Ryn.   Then  from 
equation  (1) 

"mT  "  *LO  =  R0  "  (VM  +  V*fl 

Solving  for  tf, 

R0~RLO 


t 


fl    (VM+V 


If  at  the  time  of  target  lock-on,  tf,,  the  target 
and  launch  aircraft  are  at  the  minimum  separation  distance, 
^LT  =  Rmin'  from  equation  (2) , 
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R__    =    R    .       =    Rn   -     (VT     +    Vm)t 
XT  min  0  L  T 


fl 


Inserting    for    tf, , 


Rmin  "   R0  "    <W  (V*LO)/(W 


Solving    for   R    , 


Ro  = 


mm 


1    -      /VV 

\VM+VT, 


,vm+vt;  ^ 


1    -      /VL+VT 

\VVT, 


(3) 


For  the  same  geometry  and  relative  speeds  of  the 


first  case,  with  M„  =  3.0. 

M 


•fl 


5.5a 


A  reasonable  value  of  lock-on  range  is  10  km.   This 
will  be  shown  later  in  the  guidance  section  of  the  study.   The 
time  to  lock-on  then  becomes, 


•fl 


R  -10 
_  _0 

5.5a 


From  equation  (3),  Rn  then  becomes,  R   =  118.52  km 
range  to  lock-on  is  then,  RMw 


The  missile 


R»„i  =  V„t_n  =  59.19  km 
Ml     M  f  1 
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If  the  target  does  not  maneuver  the  time  from  lock- 
on  to  intercept  becomes,  t-2« 


£2  (VM+VT) 


and  the  missile  range  from  lock-on  to  intercept  becomes, 

*M2  =  VMfcf2  =  5-45  km 

The  total  missile  range  is  then  the  sum  of  the  two, 
*M   ~   ^Ml  +  RM2  =  64-64  km  =  34-88  nmiles  . 

As  can  be  seen  from  the  above  analysis,  both  the 
detection  range  of  the  target  and  the  required  missile  range  are 
decreased  significantly  when  an  active  homing  missile  is  used. 
On  the  other  hand,  it  must  also  be  remembered  that  the  complexity 
and  cost  of  the  missile  will  be  increased  as  a  result  of  choosing 
an  active  radar  seeker.   For  the  design  example  in  this  study, 
an  active  radar  seeker  is  chosen;  therefore,  the  maximum  range 
requirement  will  be  35  nmiles  at  a  speed  of  M   =  3.0,  however, 
this  missile  velocity  is  tentative  until  a  guidance  analysis 
is  complete. 

From  the  preceding  analysis  the  mission  profile  is 
determined.   It  must  be  kept  in  mind  that   the  mission  profile 
may  be  changed  during  the  design  process  to  meet  other  design 
objectives.   The  following  profile  assumes  both  the  target  and 
launch  aircraft  at  the  same  altitude. 
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BOOST 


CRUISE 


TERMINAL 
HOMING 


SEGMENT 

begin 

begin 

end 

(1) 

Boost 

1.5 

10,000  ft 

3.0 

(2) 

Cruise 

3.0 

10,000  ft 

3.0 

(3) 

Terminal 

3.0 

10,000  ft 

3.0 

hend  *™GE 
10,000  ft 

10,000  ft  29.6  nmiles 

10,000  ft  5.4  nmiles 
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III.   GUIDANCE  LAW  SELECTION 

Although  the  specifics  of  the  guidance  system  is  beyond  the 
scope  of  conceptual  design,  the  selection  of  a  guidance  law  is 
necessary  for  initial  calculations.   The  warhead  design  depends 
on  the  expected  miss  distance  between  the  missile  and  the  target 
and  the  lifting  surface  area  depends  on  the  maneuvering  require- 
ments of  the  missile.   Both  the  miss  distance  and  maximum  acceler- 
ation required  are  functions  of  the  missile  guidance  law. 

The  guidance  law  for  a  missile  is  the  analytical  formulation 
used  by  the  guidance  system  to  convert  sensed  target  information 
into  missile  steering  commands.   Three  general  guidance  laws  are 
used.   Most  others  can  be  forced  to  fit  into  one  of  these  cate- 
gories.  These  are: 

1)  Pursuit  Guidance 

2)  Line-of-Sight  Guidance 

3)  Proportional  Guidance 

A.   PURSUIT  GUIDANCE 

A  pursuit  guidance  law  is  illustrated  in  Figure  (3-la)  and 
is  one  in  which  the  missile  velocity  vector  is  always  directed 
toward  the  target.   The  target  and  the  missile  velocity  vectors 
must  therefore  be  sensed;  so  this  type  of  guidance  normally 
assumes  an  on-board  tracker.   The  missile  may  have  a  separate 
mid-course  guidance   package  to  increase  range,  but  target  lock- 
on  initiates  the  pursuit  guidance  for  the  terminal  homing  phase. 
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Figure  (3-la).   Pursuit  guidance. 

X 


Figure  (3-lb) .   Beam  rider  guidance 


Figure  (3-lc) .   Proportional  guidance 
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For  this  reason  it  has  the  advantage  of  launch-and-forget  at 
lock-on.   Since  the  signal  processing  is  limited  to  looking  and 
pointing,  the  avionics  are  relatively  simple  and  usually  on- 
board the  missile.   An  option  for  this  type  of  guidance  would 
be  to  include  a  lead  angle  to  accomodate  faster  moving  targets. 

B.  LINE-OF-SIGHT  GUIDANCE 

Line-of-sight  guidance  is  used  in  a  beam  rider  missile.   This 
guidance  scheme  is  illustrated  in  Figure  (3-lb)  and  requires  that 
the  missile  remain  on  a  line  (beam)  joining  the  target  and  a 
control  point.   The  target  tracker  is  located  at  the  control 
point;  therefore,  avoiding  the  necessity  of  an  on-board  tracker. 
Because  of  this,  a  dedicated  fire  control  system  is  needed  from 
launch  to  intercept.   The  range  of  this  type  of  guidance  is 
normally  less  than  with  the  other  types.   A  speed  advantage  is 
required  for  line-of-sight  guidance  since  no  lead  angle  is  incor- 
porated.  The  main  advantage  of  this  type  of  guidance  is  the 
simple  avionics  required  to  maintain  the  missile  in  the  beam. 

C.  PROPORTIONAL  GUIDANCE 

A  proportional  guidance  law  is  one  in  which  the  rate  of  change 
of  the  missile  heading  is  made  proportional  to  the  rate  of  change 
of  the  line-of-sight  between  the  missile  and  the  target.   This 
is  illustrated  in  Figure  (3-lc) .   Since  the  guidance  law  antici- 
pates the  target's  future  position,  it  can  attain  a  higher  degree 
of  responsiveness  than  other  guidance  laws.   In  proportional  gui- 
dance the  rate  of  change  of  the  line-of-sight  must  be  sensed 
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on-board  the  missile.   Because  of  this  requirement,  and  the  need 
to  provide  anticipated  steering  commands,  the  avionics  required 
are  the  most  complex  of  the  three  guidance  systems. 

D.   COMPARISON  OF  GUIDANCE  LAWS 

In  early  design  considerations  two  parameters  of  interest 
are  acceleration  required  of  the  missile  and  the  miss  distance 
attainable.   Of  the  three  guidance  laws  only  the  proportional 
law  can  respond  to  fast  maneuvering  targets.   Since  the  missile 
must  stay  in  the  line-of-sight  for  a  beam  rider  system,  any  tar- 
get maneuver  will  cause  large  excursions  in  the  missile  flight 
path,  resulting  in  large  normal  accelerations.   The  pursuit  gui- 
dance law  causes  similar  large  excursions  near  intercept  due 
to  the  velocity  vector  always  pointing  at  the  target. 

Several  system  parameters  affect  the  miss  distance  attain- 
able with  a  particular  guidance  law.  An  excellent  source  on 
the  effects  of  these  parameters  is  an  article  written  by  Dr. 
Robert  Goodstein [6] .  The  parameters  studied  for  their  effect  on 
miss  distance  were: 

1)  Sensor  Bias  Angle 

2)  Noise 

3)  Target  Heading 

4)  Target  Acceleration 

5)  Target  Speed 

6)  Wind  Gusts 

The  results  have  been  reproduced  and  are  included  in  Figures  (3-2) 
through  (3-7).   Table  3-1  provides  overall  guidance  in  the  selection 
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of  a  guidance  law  and  is  also  reproduced  from  the  above  reference. 
A  first  glance  would  indicate  that  proportional  guidance  is  a 
proper  choice  for  all  cases.   It  must  be  kept  in  mind,  though, 
that  cost  and  simplicity  are  also  driving  factors  in  the  design 
process.   Furthermore,  it  can  be  seen  that,  while  proportional 
guidance  with  a  high  gain  has  good  performance  against  maneuvering 
targets,  any  noise  in  the  system  will  highly  degrade  this  perfor- 
mance.  For  this  reason  another  guidance  law  may  be  desired,  or 
a  compromise  in  the  gain  selection  may  have  to  be  made  in  which 
some  performance  is  given  up  in  order  to  deal  with  a  noisy  system. 
A  reasonable  range  of  proportionality  constants  that  gives  good 
performance  against  both  maneuvering  targets  and  noisy  systems 
is  k  =  2  to  k  =  6 . 

Once  a  guidance  law  is  selected,  a  more  detailed  analysis  has 
to  be  performed  to  determine  if  the  maximum  acceleration  required 
of  that  particular  guidance  system  is  within  the  attainable 
maneuverability  limits  for  the  missile.   The  maximum  acceleration 
required,  in  turn,  determines  the  lifting  surface  area  needed. 
A  good  estimate  of  maximum  acceleration,  which  keeps  the  miss 
distance  less  than  50  feet,  is  to  set  it  equal  to  three  times 
the  target  acceleration  plus  ten. 

a   =  3a^  +  10 
m     t 

Figure  (3-8)  shows  the  miss  distance  sensitivity  to  target 
acceleration. 
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E.   PURSUIT  GUIDANCE  (DETAILED  ANALYSIS) 

As  stated  previously,  a  pursuit  guidance  law  requires  the 
missile  velocity  vector  to  always  point  at  the  target.   For 
this  reason  the  missile  always  ends  up  in  a  tail  chase  situation, 
with  the  maximum  acceleration  occurring  at  the  end  of  the 
encounter.   From  this  description  the  maximum  acceleration  of  the 
missile  can  be  determined. 


M. 


'P   vT 


/ 
/ 

/ 


/ 


V    ' 
M. 


/v 


R. 


Reference 
Direction 


Figure  (3-9).   Pursuit  geometry. 


From  Figure  (3-9)  the  time  rate  of  change  of  the  range,  R,  is 


R  =  VTcos  6  -  VM 


also 


or 


B  =  -VT  sin  3/R 


R  *£  =  VT  sin  6 
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-M^  =  -Is  dt 

sin  6     R 


(1) 


dt   = 


dR 


VT  cos  3  -  VM 


(2) 


Substituting  equation  (2)  into  equation  (1)  gives, 


dB 


V 


sin 


T  dR 

R  (VT  cos  B  -  VM)  '* 


Letting  k  = 


M 


(cos  B  -  k)  An    _   dR 
iTn~B d3 R 


(3) 


Integrating  equation  (3)  yields, 


JinR  =  kAn(tan  =■)  -  Jin  (sin  B)  +  Anc 


AnR  =  An 


c   tarF  8/2 


sin  B 


By  trigonometric  identity, 


(tan  B/2)k  =   (sin  3) 


(1+cos  B) 


Therefore, 


AnR  =  An 


c^    (sin  B) 
(1+cos  3)1 


k-l 


and 


R  = 


c   (sin  B) 
(1+cos  B)} 


k-l 
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From  the  initial  condition  B  =  6   when  R  =  R  , 


R   ( 1+cos B  )k 
o,  =  -^ 0 


(sing  ) 


/l+cosB  \     /    •    o   \ 
R  =  R   f      0  I     /sing  \ 

0  \  1+cosB/     I  sinS  J 


Substituting  the  above  equation  for  R  into  the  equation  for  @ 


6  =  -V   sinB/R 


S  _  Vt  /l+cosB  \k    (sinB)2"k 
~R„  1 1+cos  B"J    ",~;„Q  I  l-k 


D  \      0'     (sing  ) 


The  missile  acceleration  can  be  expressed  as  a  normal  and  a 
tangential  component, 


a  =  VM  B  n  +  VM  t 
M         M 


The  normal  component  is  a  ,  where, 


a   =  V..  B 

m    M 


M  T   / 1+cos  B  \    sin  b  ^7. 

am   ~  P    I  1+cos  B  /    sins 
0    v      0  I  0 

The  terminal  acceleration  for  a  pursuit  guidance  law  will  occur 
at  the  end  of  the  encounter  (B  ■*  0)  .   From  the  above  expression 
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the  terminal  acceleration  can  be  evaluated 


a   =0,     for    1  <  k  <  2 
m 


V  V 

a  =  — §-£■  I  rr  "  a   J   sin^   for   k  =  2 
m     R0   \  1+cosB  /      o 


a   =  °°   for   k  >  2 
m 


Since  pursuit  guidance  always  ends  up  in  a  tail  chase  situation, 
the  missile  will  never  intercept  if  k  <  1  .   Thus,  for  pursuit 
guidance  operating  against  a  non-maneuvering  target,  the  velocity 
ratio  should  be  between  one  and  two.   These  results  indicate 
this  guidance  system  would  not  be  effective  against  air  targets; 
therefore,  results  for  a  maneuvering  target  were  not  pursued. 

F.   LINE-OF-SIGHT  GUIDANCE  (DETAILED  ANALYSIS) 


/ 


/ 


/  Reference 

n  Direction 


Figure  (3-10).   Line-of-sight  geometry. 
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Figure  (3-10)  illustrates  the  geometry  used  to  derive 
the  beam  rider  equations  of  motion.   The  basic  concept  of  beam 
rider  guidance  is  that  the  missile  is  maintained  in  the  line-of- 
sight  of  the  target  and  a  control  point.   This  can  be  expressed 
in  an  equation  as  follows: 


V_  sin  a.     V..  sin  a 
+  =  -2 L_=J5 m  (1) 

rt 


where,  r   =  range  from  point  0  to  the  target 
r  =  range  from  point  0  to  the  missile 

From  equation  (1) 


r.VV,  sin  a   =  r  V_  sin  a.  (2) 

t  M      m    m  T       t 


As  in  the  case  of  pursuit  guidance,  the  missile  and  target  acceler- 
ations can  be  divided  into  normal  and  tangential  components.   If 
the  target  is  limited  to  contant  g  turns,  and  the  normal  component 
of  missile  acceleration  is  of  interest;  then, 

Differentiating  equation  (2)  with  respect  to  time  yields, 


r. V„   sin  a     +   rt   V„   a      cos    a      =   r     V_   sin   a.    +  V_   a.     cos   a. 
t  M  mtMm  mmT  tTt  t 


Solving   for  a      , 
m 


•  «  »  • 

a      =   — =-= [r     Vm   sin   a.    +   r      V_   a.    cos    a,     -    r.    V.,   sin   a    ] 

nt        r.V^cosa       LmT  t  m     T      t  t  t      M  m 

t  M  m 


52 


From  the  original  figure, 


Also 


rt  =  VT   cos   a 


r      =   V.,    cos    a 
m  M  m 


•     «  a     +  <j>  *5     =  at  +   * 


)      ■  a     +   <b 8     «  a     +  $ 

m  m  mm 


The  target  and  missile  accelerations  (normal  components)  become, 


at  "  VT  9t 


a  -  v„  e 

m    Mm 


Collecting  equations; 


't  *  at/Vt 


<J>   =  VT   sin   at/rt 


.  »  i 

»t  -  et  -  ♦ 


rt  =  VT   cos   at 


r      =  V..   cos    a 
m  M  m 

i         r  * 

a™  =   — T7 r      Vm   sin   a.    +   r      V_   a.    cos    a. 

m        r,V„cosa  m     T  t  m      T      t  t 

t  M  m  L 

in    a 

mj 


-    r.    VM    s. 
t      M 


•  .  • 

9      =   a      +    d> 
m  m        r 


a      =  VM    9 

m  Mm 
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The  above  equations  are  the  equations  of  motion  which  describe 
the  target  and  missile  trajectories.   These  equations  cannot  be 
solved  analytically  except  for  highly  specialized  cases.   The  com- 
plete set  of  equations  can  be  solved  using  a  numerical  integra- 
tion technique.   If  Euler's  one  step  method  is  used,  the  algorithm 
is  as  follows; 


>t(i+D  =  et(i 


(i+l)    =      <f>(i 


a,  (i+1)    =   cud 


rt(i+l)    =   rfc(i 


r    (i+l) 

m 


rm(i 
m 


a    (i+l)    =   a    (i 
m  m 


)    (i+l)    =    e(i 
m  m 


+   At    8    (i) 


+    At      (j)(i) 


+   At  a    (i) 


+   At   rt(i) 


+    At   r    (i) 
m 


+    At   a    (i) 
m 


+    At    9    (t) 
m 


With  initial  conditions; 


rt(0)  =  rQ 


r  (0)  =  0 
m 

<j>(0)  =  <J>, 


a  (0) 

m 


=  0 


et(0) 


m 


at(0) 


m. 


=  a 
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The  target  and  missile  positions  can  be  expressed  as  follows; 

x  (i+1)  =  x  (i)  +  At  v  cos  0  (i) 
m         m          Mm 

y  (i+1)  =  y  (i)  +  At  VM  sin  8  (i) 
m         m          M      m 

xt(i+l)  =  x.(i)  +  At  VT  cos  9.(i) 

yt(i+l)  =  yt(i)  +  At  VT  sin  8  (i) 
Where 


x  (0)  =  ym(0J  =  0 

m       m 

Yt(0)    =   rQ    cos    <J>0 
yt(0)    =   rQ    sin    ^ 

The  above  equations  have  been  programmed  on  the  HP  9  8  30  com- 
puter.  Table  3-II  is  a  listing  of  this  program.   The  program 
asks  the  user  for  the  initial  conditions  and  the  target  acceler- 
ation.  It  also  asks  for  the  integration  step  increment,  At. 
It  should  be  kept  in  mind  when  using  the  program  that  the  error 
involved  in  integrating  is  of  order  At.   The  output  is  a  plot  of 
missile  and  target  trajectories  as  well  as  the  missile  maximum 
acceleration  and  time  of  flight.   Three  examples  follow  which 
demonstrate  possible  uses  of  the  program.   (Note:   All  angles 
are  input  in  radians.) 
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TABLE  3-II 


19  PRINT  "THIS  PROGRAM  DETERMINES  THE  RCCELflRRTION  OF" 

28  PRINT  "fl  BEAM  RIDER  MISSILE  AND  PLOTS  THE  TRRJECTOR 

38  PRINT 

48  PRINT  "INPUT  TIME  INCREMENT  FOR  INTEGRATION" 

50  INPUT  ru 

68  HI M  X [ 2 5 9  ] j Y C  2 5 0  3 j U C 2 5 6  ]  j V C 2 5 8  ]  > R C 2 5 8  J ?  R C 2 5 8  ] 

78  PRINT  "INPUT  INITIAL  TARGET  RANGE" 

38  INPUT  Rl 

98  PRINT  "INPUT  TARGET  SPEED" 

188  INPUT  VI 

118  PRINT  "INPUT  MISSILE  SPEED" 

128  INPUT  V2 

138  PRINT  "INPUT  INITIAL  LINE  OF  SIGHT  ANGLE" 

148  INPUT  PI 

158  PRINT  "INPUT  MISSILE  ALPHA" 

168  INPUT  A2 

178  PRINT  "INPUT  TARGET  ALPHA" 

138  INPUT  Hi 

198  PRINT  "INPUT  MISSILE  THETA" 

288  INPUT  T2 

218  PRINT  "INPUT  TARGET  THETA" 

228  INPUT  Tl 

2  3  8  P  R I N  T  "  I N  P  U 1"  T  A  R  0  E  T  A  C  C  E  L  E  R  A  T  1 0  H " 
248  INPUT  Gl 

258  1=1 

268  R2=8 

265  RE  1  ]  =  R1 

278  AC  1 3=0 

238  XC I  3=0 

298  YC I  3=0 

388  UC I 3=R1*CQS<P1) 

318  VC I 3=R1*SIN(P1) 

328  AC  I  3=0 

338  PRINT  "    NH  YM  NT 

3  9  5  P  R I N  T  ' '  D  0  Y  0  U  W  fl  N  T  A  P  R  I N  T  0  F  T  H  E  0  U  T  P  U  T  t    8  =  Y  E  S  s  1  = 

396  INPUT  T9 

397  IF  T9=l  THEN  418 

488    PRINT    NC  I  3?  YC  I  3  s  U[  I  3>VC  I  1 

418  1=1+1 

428  T3=G1/V1 

438  P3=Vl*SIN(fil>/Rl 

448  H3=T3-P3 

458  R3=V1*C0S<A1 ) 

4  g  q  R  4  =  y  2  *  C  0  S  ( H  2  •' 

461  D6=R4*V1*SIN(AD 

462  D7=R2* ' ,'  1  #  A 3 * C 0 S  (  A 1  > 

463  D8=-R3*V2*SIN<A2) 
4  f,  4  i\  q  =  r  i  +  '■/  2  *  C  0  S  (.  fl  2  ) 
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TABLE  3-II  (cont) 


465 
4  6  6 
470 
430 
4  9  0 
495 
496 

500 

510 
520 
530 
540 

600 

6 1 0 
620 
6  3  0 
640 
6  5  0 
6  6  O 
665 
670 
630 
6  9  O 

700 

710 
720 
730 
740 

300 

819    - 

3  2  0 
3  3  O 
340 
350 
3  6  0 
370 
3  8  0 
3  9  0 
9  0  0 
910 
920 
930 
940 
950 
951 


Q 


.1  £. 


fl4=(D6+D7+D3 >^D9 

T4=fl4+P3 

T1=T1+B1*T3 

Pi=pi+ru*p3 

Rl=Rl+m*R3 

R1=R1+D1*R3 

R2=R2+D1*R4 

R2=fl2+Di*fl4 

T2=T2+B1*T4 

RC I  ]=V2*T4 

IF    RB3(R[  I  3XHBS 

(AC  I- 

1  3 ) 

THEN    6 

N  =  RC  I  3 

T2=R2+P1 

RC I  ]  =  R1-R2 

X[  I  3  =  X[ 1-1 3+M*V 

2*  COS 

(T2 

'i 

YC I  ]  =  YC 1-1  3+Bl*V 

2*3  IN 

(T2 

j 

UC I 3=UC 1-1 3+Dl*V 

'l*COS 

(.  T 1 

j 

VC I 3=VC 1-1 3+Dl*V 

1*3  IN 

( T 1 

) 

IF    RC  I  ]>0    THEN    :■ 

t97 

T5=B1*<I-1> 

PRINT    "INPUT    MIH 

II  NUN 

'■:'      '■  i1 

RLUE" 

INPUT    X6 

PRINT    "INPUT    P1H>: 

1INUN 

'•■■'      'i 

RLUE" 

INPUT    X7 

PRINT    "INPUT    MINI  NUN 

1        ',' 

'RLUE" 

INPUT    Y6 

PRINT    "INPUT    Mftt 

IINUN 

1 

'RLUE" 

INPUT    Y7 

S  C  R  L  E    X  6  j  )<  7  i  Y  6  ?  'i 

i" 

PRINT    "HfiS    AXIS 

BEEN 

DRRN 

INPUT    Fl 

IF    F1=0    THEN    360 

XRXI3    0>X7/1Q 

YflXIS    0?  Y7.-10 

PEN 

FOR    N=l    TO    1 

PLOT    XCW3jY[W3 

NEXT    N 

PEN 

FOR    3=1    TO    I 

PLOT    UCS3jVCS3 

NEXT    3 

PEN 

PRINT 

PRINT 

IF    RC  I  3>Fl[  1-1  3    " 

fHEN    S 

'54 

M=RC 1-13 

PRINT    "THE    MAXIMUM    AC 

:CEl 

.ERflTIOh 

PRINT    "THE    TINE 

TO       INTERCEPT 

d    STOP 

9! 

9  5  4  P  R  I N  "  T  H  E  M  fl  X  I M  U  M  fl  C  C  E  L  E  R  fl  T  1 0  N 

960  PRINT  "THE  TINE  TO   INTERCEPT  IS 

1000 
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1.      Example    I    (Non-maneuvering     Crossing  Target) 

^<T 

\ 

\ 

\ 


•VM 


»0 m^ 


Figure  (3-11) .   Non-maneuvering  crossing  geometry 


r  =  359  meters 
0 

<J>   =  109.57°  =  1.9124  rad 
9    =  17.57°  =  .3067  rad 

9    =  109.57°  =  1.9124  rad 
mQ 

a    =  -91.0°  =  01.5882  rad 
fc0 

at  =  0 

v  =37  3  m/sec 
M 

V^  =  221  m/sec 

Table  3-III  is  the  computer  output.   As  indicated  the  missile 
maximum  acceleration  is, 


a 
m 


=  -459.25  m/sec/sec  =   -46.86  g's 
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TABLE  3-III 


INPUT  TIME    INCREMENT    FOR    INTEGRATI 

INPUT  INITIAL    TARGET    RANGE 

INPUT  TARGET    SPEED 

INPUT  MISSILE    SPEED 

INPUT  INITIAL    LINE    OF    SIGHT    ANGLE 

INPUT  MISSILE    ALPHA 

INPUT  TARGET    ALPHA 

INPUT  MISSILE    THETA 

INPUT  TARGET    THETA 


INPUT    TARGET    H 

.JCELERATION 

XM 

YM 

.'■,  1 

T 

DO    YOU    WANT    fl 

PRINT    OF    THE    OU 

TPUTj    0=YESj 1=N0 

0 

8 

-120,,  2644734 

3  3  8 . 

256495 

-5. 154984716 

17.92341813 

-189. 7381283 

341. 

5926479 

-9. 197499767 

36. 1300131 

- 9 9 . 1 9 57681 8 

344. 

5 ':.:  3  8  8  8  7 

-12. 11631381 

54. 55819713 

-88. 66141559 

343. 

2649534 

-13.90  8 8 8 092 

73. 11335036 

-78. 127863 

3  5 1 ,, 

6  8 1 1 6  6  2 

-14.53120543 

91.75172791 

--67  „  5927184 

354. 

937259 

-14. 14756982 

110. 3966859 

-  J  l  .  fcj  ■  J  O  J  -J  t   O  £ 

353. 

2734118 

-12. 62997166 

123.9348330 

-  4  6=52  4  8 3  5  2  3 

361. 

6  8  9  5  6  4  6 

-10. 0574053 

147.456557 

- 3 5 . 9 8 9 65264 

364. 

9457174 

- 6 . 46499644  8 

165.7572973 

-25. 455  3 8885 

3  6  8 . 

2313782 

-1.8  9  3 023052 

133.3332156 

-14. 92894745 

371. 

6  1  3  8  2  3 

3. 614079328 

201.6565344 

- 4 . 3 8 6594  8 6  3 

374. 

9541758 

In. 90880633 

219, 1760013 

6.  147757723 

••-,  -* ^  i", 

2983236 

17.24154285 

2  3  6 .  3  6  6  4  8  6  2 

16.68211032 

381 ! 

6  2  6  4  8 1 4 

25.26153519 

253. 2839385 

27.21646291 

334, 

9  6  2  6  3  4  2 

34. 01774632 

269. 67  0 6  0 1 9 

37.7588155 

298737 

4  3 . 45961 8  9  9 

■"'  O  ^j  ~P  cj  O  C|  ■'"'  "7  "7 

48. 28516389 

39 1! 

6  3  4  9  3  9  8 

53.53767559 

301 . 4464555 

5  8 . 8 1 9 5286  8 

394, 

9  7  16926 

64.2841 194 

316.745131 

69. 35387327 

3  9  8  i 

3872454 

75.41315987 

331.6503313 

-".'  i'j   i'-i  i-i  Q  •-,  •-,  C  O  "7 

i  7  „  O  O  O ,-  C  J  O  i 

481  , 

6  4  3  3  9  8  2 

87. 12134988 

346.  16  7746/' 

90. 42257346 

484, 

r"i  ""'  i'"i  C  £T  -i 

Cj  Q   v  •  j  7  7  7  1  CT  -~i 

360. 3028325 

188. 956931 

488, 

3157833 

111.8741571 

374. 0652959 

111  ,,4912336 

411 

6513566 

124. 3449199 

337. 4661897 

122 „ 8256362 

4  1  4 

9  J 38894 

133. 1671393 

480.5175  8 3 6 

132.5599388 

4  1 8 

3241622 

151.8184937 

413.232983 

143.8943414 

421 

668315 

INPUT    MINIMUM 

X    VALUE 

INPUT    MAXIMUM 

X    VALUE 

INPUT    MINIMUM 

Y    VALUE 

INPUT    MAXIMUM 

Y    VALUE 

HAS    AXIS    BEEN    DRAWN,    0=YESj     1=N0 


THE    MA X I M U M    A C C E L E R A T 1 0 N    I  3         - 4 5 9  .  2 5 2 6872  M E T E R S  ■ " 3. E C ••■ ' 3 E i 

THE    TIME    TO       INTERCEPT    IS  1.3  SEC 
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The  large  acceleration  is  typical  of  line-of-sight  due  to  the 
missile  requirement  to   stay  in  the  beam.   The  trajectories  are 
plotted  in  Figure  (3-12) . 

2.   Example  II  (Effect  of  VM) 

This  example  is  presented  to  study  the  effect  of  missile 
velocity  on  a  crossing  target  (non-maneuvering) . 

a, 


/ 


VM    / 


/ 


M 


*o  =  V 


Figure  (3-13).   Crossing  target  geometry 


The  initial  conditions  are  as  follows; 


rn  =  4000  m 


cf)   =  45°  =  .7854  rad 


9    =  45   -  .7854  rad 


m, 


a, 


=  135   =  2.3562  rad 


=  180   =  3.1416  rad 


0 


Vm   =  200  m/sec 
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TABLE  3- IV 


THIS  PROGRAM 

A  BEAM  RIDER 


DETERMINES  THE  AC 

MISSILE  HMD  PLOTS 


^ELA RATION  OF 
THE  TRAJECTORY 


JT 

JT 

JT 

IT 

IT 

IT 


INPUT 
INPUT 
INP 
INPI. 
INP  I 
INP  I 
INPI. 
INPI. 
INPUT 
INPUT 
XM 
DO  VOL 
INPUT 
INPUT 
INPUT 
INPUT 
HAS  A',; 


TIME  INCREMENT  FOR  INTEGRATION 

INITIAL  TARGET  RANGE 

TARGET  SPEED 

MISSILE  SPEED  =  Yoc   sn/s&c 

INITIAL  LINE  OF  SIGHT  ANGLE 

MISSILE  ALPHA 

TARGET  ALPHA 

MISSILE  THETA 

TARGET  THETA 

TARGET  ACCELERATION 

YM 
I  WANT  A  PRINT  OF  THE  OUTPUT.  8 
MINIMUM  X  VALUE 
MAXIMUM  X  VALUE 
MINIMUM  Y  VALUE 
MAXIMUM  Y  VALUE 
as  BEEN  DRAWN)  0=YESj  l=NO 


'ESs  l=NO 


THE 
THE 


MAXIMUM 
TIME  TO 


ACCELERATI 
INTERCEPT 


IS 


03129548 

SEC 


METER 


3Ei 


THIS  PRi 
A  BEAM  I 


GRAM 
IDER 


DETERMII 
MISSILE 


IES-THE  AC 
AND  PLOTS 


JELARATION  OF 
THE  TRAJECTORY 


INPI 

INPI 

INPI 

INPI 

INPI 

INP 

INPI 

INPI 

INPI 

INPI 

DO 

INP 

INP 

INP 
INP 
HAS 


JT 
JT 

JT 
JT 
JT 
JT 
LIT 
UT 
UT 
UT 
XM 
YOI 
UT 
UT 
UT 
UT 
A 


TIME  INCREMENT  FOR 
INITIAL  TARGET  RAM 
TARGET  SPEED 
MISSILE  SPEED  = 
INITIAL  LINE  OF 
MISSILE  ALPHA 


INTEi 

IE 


IRATIC 


C  OO     rr\  /  S£~C 

SIGHT  ANGLE 


TAR 
M I S 
TAR 
TAR 


GET  ALPHA 
SILE  THETA 
GET  THETA 
GET  ACCELERATION 
YM 

PRINT  OF  THE 

X  VALUE 

X  VALUE 

Y  VALUE 

Y  VALUE 
DRAHN?  U 


J  WANT  A 
MINIMUM 
MAXIMUM 
MINIMUM 
MAXIMUM 

aS  BEEN 


TPUTj  0: 


:t 

'ES 


r'T 


l=Nl 


=YESj  l=NO 


THE 
THE 


MAXIMUM 
TIME  TO 


ACCELERATIC 

INTERCEPT 


IS 


774168 

SEC 


METERS/SEC/SEC 
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TABLE  3-IV  (cont) 


THIS  PROGRAM 
R  BERN  RIDER 


D  E  T  E  R  M I H  E  S    T  H  E    fl  C  C  E  L  R  R  R  T  1 0  N 
MISSILE    RND    PLOTS    THE    TRflJEi 


:TOR' 


INPUT 

TIME  INC 

JREMENT  FOR  INTEl 

IRRTIC 

INPUT 

INI 

TIRL 

TARGET  RANGE 

INPUT 

TRR 

GET  '■ 

SPEED 

INPUT 

M I  £ 

SILE 

SPEED  =   800 

no 

/sec 

INPUT 

INITIAL 

LINE  OF  SIGHT 

Rl 

•IGLE 

INPUT 

M I  £ 

SILE 

ALPHA 

INPUT 

TfiF 

GET  ALPHA 

INPUT 

M I S 

SILE 

THETA 

INPUT 

TRR 

GET  " 

rHETfl 

INPUT 

TRF 

GET  1 

ACCELERATION 

XM 

VM 

DO  YOl 

J  WF 

NT  R 

PRINT  OF  THE 

DUTPUTj 

INPUT 

M I  h 

I  MUM 

X  VALUE 

INPUT 

MA> 

:imum 

X  VALUE 

INPUT 

M I  h 

II  MUM 

Y  VALUE 

INPUT 

MR> 

IIMUM 

Y  VALUE 

HRS  r; 

<  I S 

BEEN 

DRAWN)  0-YE8? 

1 

=NO 

THE  MAX  It 

1UM  R 

:CELERRTIUN  IS 

O  (      u 

THE  TIME 

TO 

INTERCEPT  IS 

4.4 

)N 


9= YES* i=N0 


22950S71 
SEC 


METE! 


oci 
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The  program  was  run  three  times  for  a  missile  velocities  of 
400,  600,  and  800  m/sec.   Table  3-IV  contains  the  program  out- 
puts and  Figure  (3-14)  is  the  plot  of  the  trajectories.   As 
can  be  seen  from  the  output,  the  missile  maximum  acceleration 
increases  with  increasing  missile  speed.   For  example, 

V,./Vm  =  2  gives  the  smallest  acceleration,  although  the  maximum 
M  T 

acceleration  for  V  /V_  =  4  is  not  exceedingly  large  for  this 
scenario. 

3 .   Example  III  (Maneuvering  target) 

In  this  example  the  effect  of  a  target  maneuver  is  inves- 
tigated.  If  at  the  time  of  launch  the  target  initiates  a  7  g 
(68.6  m/sec/sec)  turn,  the  following  encounter  would  result: 


V 


M  ^ 


\< 


*o  - 


m. 


Figure  (3-15).   Maneuvering  target 
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TABLE  3-V 


THIS  PROGRAM 
fl  BERM  RIDER 


DETERMINES 
MISSILE  fih 


THE  RCCELRRflTION  OF 
D  PLOTS  THE  TRRJECTOR' 


I  HP  I 
INPI 
I  HP l 
I  HP! 
INPI 
INPI 
I  HP 
I  HP 
I  HP 
I  HP 

DO 

I  HP 
I  HP 
I  HP 
I  HP 
HAS 


I  j 

JT 

JT 

JT 

JT 

JT 

LIT 

LIT 

LIT 

LIT 

XM 

VOL 

UT 

LIT 

UT 

UT 


TIME  INCREME 
INITIRL  TARGET 
TARGET  SPEED 
MISSILE  SPEED 
INITIRL  LI  HE  OF 
MISSILE  ALPHA 
TARGET  ALPHA 
MISSILE  THETA 
TARGET  THETA 
TARGET  RCCELERR 
YM 


T  FOR  INTEGRRTIO 
T  RANGE 


IGHT  ANGLE 


I  OH 


YT 


1  WANT 

M I H I M 
MAXIM 
M I H I M 
MAXIM 
as  BEc 


R 


THE 


PRINT  OF 
X  VALUE 
X  VALUE 

Y  VALUE 

Y  VALUE 

DRAWN?  0=YESh  l=NO 


ITPUTj  0= 


I  CO 


1=NC 


THE  MAXIMUM 
THE  TIME  TO 


ACCELERATION 
INTERCEPT  I: 


321? 
SEI 


METER 3 /SEC:/ SEC: 
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The  initial  conditions  are; 


r   =  10000  m 


VT  =  821.436  m/sec 


*   =  90   =  1.5708  rad 


v.,  =  985.723  m/sec 
M 


"b 


=  90   =  1.5708  rad 


c^   =  -180 


=  -90 


"0 


a   =  68.60  ra/sec/sec 

From  the  output  (Table  3-V) ,  notice  the  large  missile  acceler- 
ation required  (211.55  m/sec/sec)  to  intercept  a  maneuvering 
target.   The  trajectories  are  plotted  in  Figure  (3-16). 

G.   PROPORTIONAL  GUIDANCE  (DETAILED  ANALYSIS) 


V. 


Reference 
Direction 


Figure  (3-17).   Proportional  guidance  geometry 
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Proportional  guidance  automatically  establishes  a  lead  angle 
and  reacts  to  a  changing  line-of-sight.   The  basic  guidance  law 
used  equates  the  rate  of  change  of  the  missile  heading  to  a 
constant  times  the  rate  of  change  of  the  line-of-sight.   From 
the  above  figure  this  law  can  be  expressed  as , 


)   =  k  a 
m 


From  the  figure,  the  rate  of  change  of  the  line-of-sight,  a,  is 

given  by, 

Vm  sin  B   -  v„  sin  8 
•  _  _T t    M _m 

r 

As  in  the  case  of  pursuit  and  line-of-sight  guidance,  the  para- 
meter of  interest  here  is  the  normal  acceleration;  therefore,  the 
missile  and  target  tangential  accelerations  are  assumed  to  be 
zero.   In  this  case, 


m    Mm 


at  =  VT  9t 


From  the  guidance  law, 


and 


a  =  VM  k  a 

m    M 


a  r  =  V,,  k(Vm  sin  B.  -  VM  sin  B  ) 
m     M    T       t    M      m 
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Since  V  =  V.  =  0 ,  the  time  derivative  of  this  equation  is. 
m     t  -a.  i 


r   a     +  r  am  =    k  Vm     [(VS.    cos    B.    -   V   8      cos    &)  3  (1) 

m  m  mtt  tmm  m 


From  the  original  figure, 


Also, 


e     = 

m 

6 
m 

+   o 

• 
6m  = 

• 

0 
m 

• 

-   o  = 

» 

e     - 

m 

0 
m 

k 

sm  = 

• 

6 
m 

"4> 

h~ 

* 

9t 

» 
-   a 

K- 

K 

k 

Making  these  substitutions,  equation  (1)  becomes, 


m  ■  1 

r   a      =   -a      r   +  k-  VM  [  V_  O.-t—)  cos    B.    -   VAyl0    (l-f)cos    8m 
m  m  -MTtk  t  Mmk  m 


Since, 


a   =  VK„  0 
m    Mm 


»t  -  VT  9t 


r  a   =  -a   r  +  k  VM  a.  cos  6.  -  k  V..  a_  cos  6^ 
m     m      -Mt      t      Mm      m 


-  a   (V_  cos  B.  -  VM  cos  6  ) 
m   T       t    M      m 
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From  Figure  (3-17) , 


r  =  V_  cos  B   -  v   cos  & 
T      t    M      m 


r  a   =  -2  am  r  +  k  VM  a.  cos  6,  -  k  VM  a   cos  6 
m       m        Mt       t      Mm      m 


k  Va. cosB.    am 

a   = —  [  2r  +  k  VMcos6 

m       r  r  w    m 


Collecting  equations; 


r  =  VT  cos  8t  -  VM  cos  6m 


Bt  =  6t  -  a 


m    m 


a  = 


VTsinSt  -  VMsinSm 


at 
9^  =  ^ 


9   =  k  a 
m 


The  above  equations  are  the  equations  of  motion  for  a  missile 
using  proportional  navigation  assuming  constant  missile  and 
target  speeds.   As  with  the  line-of-sight  equations,  the  motion 
is  quite  complex.   The  equations  cannot  be  solved  analytically 
except  for  special  cases.   One  such  case  will  be  investigated 
here.   That  is  for  a  non-maneuvering  crossing  target. 
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1.   Example  IV  (Non-Maneuvering  Crossing  Target) 


V, 


T 


V. 


M 


(at  =  0) 


m 


Figure  (3-18).   Non-maneuvering  target. 


For  this  case,   r  =  -V„  cos  3     from  equation  (1) 

M  m 

S      =   -!S      [2r   -    kr] 

m  r 


3-1    (k-2) 
m 


&na      =    (k-2) £nr  +  flic, 
m  1 


If 


a      =   a»      at      r  =   r_ 
m  0  0 


k-2 


a      =    aA(~) 
0   rQ 


From  this  equation, 
if  k  >  2 


m 


0   as   r  -*  0 
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if  k  =  2    a   =  a_  =  constant 
m    0 


if   k  <  2    a  — *  °°   as   r  -*  0 

m   ' 


From  the  above  example  it  can  be  realized  that  the  pro- 
portionality constant,  k,  must  be  greater  than  two.   A  more 
general  analysis  of  the  equations  of  motion  can  be  obtained  by 
solving  the  equations  numerically.   The  same  Euler's  one  step 
method  is  used  here  with  initial  conditions,  at  t  =  0, 


a  (0)  =  a 

m        in 
0 


r(0)  =  rQ 


3.  (0)  =  3. 
0 

m       mQ 


o(0)    =      aQ 


6.(0)  =  9 


6  (0)  =  9 
m       nu 


The  algorithm  used  is  as  follows: 


r(i+l)  =  r(i)  +  At  r(i) 


a  (i+1)  =  a  (i)  +  At  a  (i) 

mm  m 


a(i+l)  =  a(i)  +  At  a(i) 

)  (i+1)  =  9  (i)  +  At  9  (i) 
mm  m 
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it(i+l)  =  6t(i)  +  At  9t(i) 


6(i+l)  =  8  (i)  +  At  8  (i) 
mm  m 


8t(i+D  =  Bt(i)  +  At  8t(i) 


The  trajectory  of  the  missile  and  target  can  be  deter- 
mined assuming  the  missile  is  at  the  origin  at  t  =  0.   From  the 
original  figure  the  missile  and  target  positions  are  given  by, 

x  (i+1)  =  x  (i)  +  At  V„  cos  9  (i) 

m         m  Mm 

y  (i+1)  =  y  (i)  +  At  VM  sin  9  (i) 

m        u  m  M      m 

xt(i+l)  =  Yt(i)  +  At  VT  cos  0  (i) 
yt(i+l)  =  yt(i)  +  At  VT  sin  0t(i) 

With  initial  conditions, 

x  (0)  =  ym(0)  =  0 

m      ■*  m 

xt(0)  =  rQ  cos  aQ 
yt(0)  =  rQ  sin  aQ 

Evaluating  the  missile  initial  acceleration  can  be  more  complicated 

One  procedure  which  is  both  realistic  and  of  interest  is  to  have 

the  missile  and  target  on  a  constant  bearing  -  decreasing  range 

course  (a   =  0)  when  the  target  initiates  a  constant  g  turn  at 

t  =  0.   In  this  case  a  (0)  =  0  and  the  subsequent  motion  can  be 

m 

found. 
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TABLE  3-VI 


Pf 


18 
29 
36 
35 

36 

40 
58 
68 
78 
88 
99 

1  8  0 

118 

128 

130 

140 

150 

168 

178 

188 

288 

218 

228 

238 

248 

258 

268 

278 

288 

298 

295 

296 

297 

388 

318 

328 

338 

348 

345 

358 

360 

365 

370 

3  8  0 

3  9  0 

408 

4 1 8 

428 


RINT  "THIS 
RINT  "MISS 
RINT  "PROP I 
RINT 

IM  XL  £50]  j  YC250  3 
IM  flC250]iT[258] 


PROGRAM  FINDS  TH 
LE  ACCELERATION 
RTIONAL  HAVIGATI 


NUN 


ON 


I  Cl 


T 


RINT  "INPUT 

NPUT  Dl 

RINT  "INPUT 

NPUT  K 

RINT  "INPUT 

NPUT  VI 

PRINT  "INP 

INPUT  V2 

PRINT  "INP 

INPUT 

PRINT 

INPUT 

PRINT 

INPUT 

PRINT 

INPUT 

PRINT 

INPUT 

PRINT 

INPUT 

PRINT 

INPUT 

PRINT 

INPUT 

1  =  1 

PRINT 

PRINT 

PRINT  " 

XC I  3=0 

VC I  3=0 

UC I  ]  =  RC 1 3*COS 

VC I 3=RE 1  3*31  N 

PRINT  "DO  YOU 

INPUT  Ql 

IF  Ql=l  THEN 


3?VC 


:nen 


TINE  INi 

NAVIGATION  CONSTANT 

MISSILE  VELOCITY" 


"I 

fll 

"INPi 

AC  1  ] 

"INPI 

RC  1  ] 

"INPUT 

61 

"INPUT 

B2 

"INPUT 

Tl 

"INPUT 

T2 

"INP 

SI 


;et 


[HI 


T 


TAR 
TAR 

I N I  T 

I N I  T 

BETH 

BETH 

THETA 

T  H  E  T  A 

SIGMA 


HL 


'Al 


ELOCITY" 

CCELERATION" 

MISSILE  At CELERP1 ION" 

RANGE" 


GET" 
■■1  ISSUE" 
TARGET" 

MISSILE 


MISSILE 


IT 


"ARGET    POS 


T  C 


S 1  > 
S 1  ) 
WANT 


A    P  R I  N  T    0  F    T  H  E    0  U  T  P  U  T  »    8  =  Y  E  S  j  1  =  N  0 ' ' 


XC  I  3jYC  I  3»UC  I 


PRINT 

1  =  1  +  1 

XC I ]=XC 1-1 3+ 

YC I 3=YC 1-1 3+ 

UC I 3=UC 1-1 3+ 

VC I 3=VC 1-1 3+D1 

R2=V2*C0S(B1 >-Vl* 

A8=K*V1*A1*C0S(B1 


'C  I  ] 


1*C0S<T2)*D 

1*D1*SIN<T2 

COS<Tl 

S I  !-•!  (  T 1 
0SCB2) 
/RC  I--1 


l 
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TABLE    3-VI     (cont) 


43@  fi  9 = R  [  I  - 1  ]  *  '•  2  *  F:  2  +  K  *  V 1  * 

2  0  S  (  B  2 

,  •■  '  • 

440  fi:j  =  Ho-fiy 

445  TC I 3=<I-1)*D1 

450  S2=<V2*SIN<B1>-V1*SIN 

-■  B  2  '  '  ■• 

-R[ 

460  T8=fll/V2 

470  T9=K*S2 

430  B3=T3-S2 

490  B4=T9-S2 

500  R[ I 3=RE 1-1  ]+Dl*R2 

510  AC  I  ]=flC I-i  3+Bl*A3 

520  S1=S1+B1*S2 

530  T1=T1+IU*T8 

540  T2=T2+D1*T9 

550  B1=B1+B1*B3 

560  B2=B2+B1*B4 

570  IF  fit  I 3<AC 1-1  ]  THEN  6 

y  u 

530  Z1=AC  I  ] 

600  IF  PC  I  ]>0  THEN  350 

601  IF  ABSCAC I 3><flBS<flC I - 

1  3)  Tl 

HEN 

602  Z1=RC  1-1  ] 

603  Z2=TC I  ] 

613  PRINT 

619  PRINT 

620  PRINT  "MAX I MUM  MISSILE  flCC 

IS 

640  PRINT  "MISSILE  TIME  C 

F  FLI 

:ht 

700  PRINT  "INPUT  MINIMUM 

VALUE 

of 

710  INPUT  X6 

720  PRINT  "INPUT  MAXIMUM 

VALUE 

OF 

730  INPUT  X7 

740  PRINT  "INPUT  MINIMUM 

VALUE 

OF 

750  INPUT  Y6 

760  PRINT  "INPUT  MAXIMUM 

VALUE 

OF 

770  INPUT  Y7 

7 3 0  S C  A LE  X 6  j  X  7  ?  Y  6  >  Y  7 

790  PRINT  "HAVE  THE  AXIS 

BEEN 

DRA 

791  INPUT  Fl 

792  IF  F1=0  THEN  330 

800  X fl X IS  0 ?  X 7 ■■•'  1 0 

810  VAX  IS  0jY7.-"10 

330  PEN 

340  FOR  J=l  TO  I 

350  PLOT  XCJ3»YCJ3 

860  NEXT  J 

870  PEN 

830  FOR  W=l  TO  I 

390  PLOT  UCNIj  VC  M  ] 

900  NEXT  W 

910  PEN 

1O00  STOP 

RC  I- 

1  ] 


1  ] 


METERS--  SEC-" SEC 


IWNj    0=YES»1=NC 
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The  above  algorithm  was  programmed  on  the  HP  9  8  30  computer 
The  listing  is  included  in  Table  3-VI.   The  inputs  and  output 
of  the  program  are  the  same  as  the  line-of-sight  guidance  pro- 
gram. Several  examples  follow  which  demonstrates  the  use  of  the 
program. 

2.   Example  V  (Crossing  Maneuvering  Target) 

afc  =  16 .0  g's 


a 


0 


Figure  (3-19).   Crossing  maneuvering  target. 


V%.  =  373  m/sec 
M 

VT  =  221  m/sec 


a  =  156.8  m/sec/sec 


rQ  =  359  m 


m, 


=  11.51      =  1.2491  rad 


9    =  11 .51      =  0.3067  rad 

$    =  -91.0°  =  -1.5882  rad 
t0 


8    =  -38.0   =  -0.6632  rad 

mo 

a   =  109.57°  =  1.9124  rad 
o 


The  output  is  listed  in  Table  3-VII.   An  important  aspect  of  this 
problem  is  the  maximum  acceleration  required  (202.92  m/sec/sec). 
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TABLE  3-VII 


THIS  PROGRAM  FINDS  THE  MAX  I  NUN 

MISSILE  ACCELERATION  FOR  fl 

P  R  0  P  0  R  T 1 0  H  A  L  N  A  V I G  A  T 1 0  N  S  Y  S  T  E  M 


INPU 
IHPU 
INPU 

INPU 
INPU 
INPU 
INPU 
INPU 
INPU 
INPU 
INPU 
INPU 


T  TINE  INCREMENT 

T  NAVIGATION  CONST ANT 

T  MISSILE  VELOCITY 

T  TARGET  VELOCITY 

T  TARGET  ACCELERATION 

T  INITIAL  MISSILE  ACCEI 

T  INITIAL  RANGE 

T  BETA  TARGET 


T  BETA  MISSILE 
T  THETA  TARGET 
T  THETA  MISSILE 
T  SIGMA 
MISSILE  POSIT 


IATIC 


TARGET  POSIT 


DO  YOU  WANT  A 

0 

■j  •  y  9 1'  t'  y  y  o  y  2 
11.72925942 
17. 50036490 
23.21064114 
28.35360473 
34.44056325 
39. 95057932 
45.33038997 
50. 71944940 
55. 95494106 
61 . 07134968 
66. 053071 1 
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a  O  r 


9563' 


75.53054295 

79. 93371892 

34,21557618 

38. 20169238 

91 . 91709127 

9  5 . 3  3  6622  3 3 

93. 43536132 

101 . 1890036 

103. 5742899 

185. 5693244 

107. 1539546 

103. 310O167 

109. 0215346 

109. 2743269 

109.0535142 

108.3634271 

107. 1324200 

105.51O1054 

103.3425413 

100. 6769169 


Y  1 

PRINT  OF  THE  0 

0 

17.  69  326338 
35.48776 6  2  6 
53. 14233894 

7  0  .  3  :  •'  6  C  9  1  1  0 

88. 67891421 
106. 4659  r  49 
124.2834453 
142. 1255214 
159. 9  9  4  9  6  3 
177, 3950225 
195.32934 
213.8018172 
231.8164631 
iT^u  y  r  r'  y  i  y  t 
267.9877605 
236. 1512922 
384.378331 
322. 6464  9  8 9 
340. 9303234 
359. 3710956 
377.3166893 
396.3135251 
414.3565113 
4  3  3 .  4  3  9  8  6  8  9 
452. 0  5  3  2  0  3 9 
470. 63  9  6  2  6  3 

4  i  9 .  3  3  7  i  O  6  2 
507. 9866517 
526. 623 6  9 4  3 
545, 2362632 
563.811135 

5  ::  c  „   .:.  4  F  4  6  1 

6  0  m  .  7  9  3  2  6  5  7 


I  c 


jTPIJTj  9= YES;  1=N0 

-120. 2644734 
-109, 73O1208 
-99,32872172 
39. 84937433 
-  ?  8  »  9  2  9 k 1 8  5  2  4 

•  6  8 .  9  ?  2  3  4  7  7  6 
•59.  19193143 

•  4  9 . 6  0  O  0634  5 
•40. 20831377 

•  3  1  »  0  2  i  ,~l  9  9  "■_■  S 
-22.07517195 
-13.3555933 
-4.382251242 

3. 334206  794 
11.2834366 
18. 95543527 
26. 34854372 
33. 42943337 
48. 21332036 

4  6 , 6  8 3 52174 

5  2 . 3  3 14  46  2 
58.650  8 5  6  8  9 
64. 13293155 
69.27127132 
74.05941184 
73. 49132646 
32.56143373 
86.26462715 
89. 59623167 
92'  ,.  55206 
95.  12339263 
97. 32198776 

9  9  .  1 3  8  O  8  4  9  4 

1 0  0 . 5  5O409 


3  3  8 , 

256495 

341, 

5926479 

345, 

3  0  0  3  3  U  6 

349, 

S  r'  4'c'  r'  f  f 

3  5  3 , 

31 1 1620 

6  0  3 601 

3  6  3 . 

7461642 

3  6  9 . 

2323885 

y  7  5  „ 

0553462 

381. 

207734O 

y  y  r' . 

6  8  1  3  O  2 1  j 

394. 

4  6  9  4  O  3  6 

481  . 

561997  6 

408, 

9  5  O  6  5  9  O 

416. 

6  2  6  0  9  0  3 

424,, 

5  7  3  6  3  3 

432, 

'?  Q  O  ~.<  C>  d  \.< 

1   ."  '._'  C.  '_'  U  u_ 

441,, 

2746385 

449, 

9971918 

458. 

9548139 

468. 

136233 

477. 

5380453 

437. 

1 2  4  2  8 1 5 

496. 

9  O  6917 

586, 

8  6  5  6  426 

5 1 6 , 

9  8 7926  € 

c  -i  -:' 
•J  i  f  ■ 

2610316 

tr  -i  "7 
■J  O  i  « 

6720  3  0  4 

548. 

2073223 

558, 

3551496 

569 

60O6142 

530 

4  3 06946 

591 

3317626 

682 

2  9 01O10 
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TABLE  3-VII  (cont) 


MAXIMUM  MIS 
MISSILE  T I  r- 
INF' LIT  MINir 
INPUT  MAX  It' 
INPUT  MI  Nit' 
INPUT  MAX If 


SILE 

E  OF 

UM 

UN 

UM 

UM 


HCC  IS 
FLIGHT 


'FlLUE 
'flLUE 
'RLUE 

'flLUE 


OF 
OF 
OF 
OF 


2  0  2 

is 


9212: 
1.7 


METERS-- SEC-  S 


z<  c 


SEC 


HAVE    THE    AXIS    BEEN    DRAWN)    0: 


la  NO 
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w 
u 

Q) 
■U 


■P 

a 
a) 
o 

a> 
-p 

c 

•H 

0^ 
C 
•H 
M 
<D 
> 

c 
id 

6 


c 

•H 
CO 
03 
O 
U 

o 
o 

CO 
CD 

•H 
M 
0 

+J 

o 

<D 
•n 
<8 

u 

Eh 


o 

CN 
CO 


•H 
fa 
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The  same  problem  was  run  for  line-of-sight  guidance  with  no 
target  acceleration  in  Example  I.   The  maximum  acceleration 
was  459.25  m/sec/sec.   This  points  out  the  advantage  of  pro- 
portional navigation  over  line-of-sight  guidance.   Figure  (3-20) 
is  a  plot  of  the  trajectories. 
3.   Example  VI  (Effect  of  k) 

This  example  demonstrates  the  effect  of  varying  the  pro- 
portionality constant,  k.   The  scenario  is  as  follows. 


7  g  turn 


VL  = 


V   = 
m 


Figure  (3-21).   Initial  geometry. 

20  8m/sec 
413  m/sec 


a.  =  68.60  m/sec/sec 


=  -106   =  -1.85  rad 


m, 


=  -29   =  -.5061  rad 
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TABLE  3-VIII 


INPUT  TIME  INCREMENT 

I N P U T  N fl V IGflTI 0 N  C 0 N S T fl H T  =  3 

INPUT  MISSILE  VELOCITY 

INPUT  TARGET  VELOCITY 

I N P U T  TARGET  R C  C  E  L E  R fl T 1 0  H 

INPUT  INITIAL  MISSILE  ACCELERATION 

INPUT  INITIAL  RANGE 

INPUT  BETA  TARGET 

INPUT  BETA  MISSILE 

INPUT  THETA  TARGET 

INPUT  THETA  MISSILE 

INPUT  SIGMA 

MISSILE    POSIT  TARGET    POSIT 

,■•',  1  i  1  ■•■•,  d  i  c 

D  0    Y  0  U    W ANT    A    P R I N  T    0  F    T  H  E    0  U T  P U T  j     @  =  Y E  S  > 1  =  N 0 

MAXIMUM  MISSILE  ACC  IS  S3,,  70521471     METERS/ SEC-SEC 

MISSILE  TIME  OF  FLIGHT  IS    3.2     SEC 

INPUT  MINIMUM  VALUE  OF  X 

INPUT  MAXIMUM  VALUE  OF  X 

INPUT  MINIMUM  VALUE  OF  Y 

INPUT  MAXIMUM  VALUE  OF  Y 

H  A  V  E  T  H  E  A  X  I S  BEE  N  D R A  W  N j  O  = Y  E  S  ,  1 = N 0 

THIS  PROGRAM  FINDS  THE  MAXIMUM 

MISSILE  ACCELERATION  FOR  A 

P  R  0  P  0  R  T  1 0  N  A  L  H  A  V I G  A  T  1 0  N  S  Y  S  T  E  M 

INPUT  TIME  INCREMENT 

INPUT  NAVIGATION  CONSTANT  ■  V 

INPUT  MISSILE  VELOCITY 

INPUT  TARGET  VELOCITY 

INPUT  TARGET  ACCELERATION 

I N  P U T  INITIAL  M I S  S I L  E  A  C C E  L  E  R A  T I  0  N 

INPUT  INITIAL  RANGE 

INPUT  BETA  TARGET 

INPUT  BETA  MISSILE 

INPUT  THETA  TARGET 

INPUT  THETA  MISSILE 

INPUT  SIGMA 

MISSILE  POSIT  TARGET  POSIT 


T  U 


XI  Yl  X2 

W  ANT  A  P  R  I N  T  0  F  1"  H  E  0  U  T  P  U  T  >  Q  =  Y  E  S  ?  1  =  N  0 


M A X I M U M  M I S S I L E  A C C  IS  71.11 3 3 3 1 9 7     METE P S / S E C • ' S E C 
MISSILE  TIME  OF  FLIGHT  IS    3.2     SEC 
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TABLE  3-VIII  (cont) 


THIS  PROGRAM  FINDS  THE  MAX  I  MUM 

MISSILE  ACCELERATION  FOR  R 

P R 0 P 0 R T 1 0 N A L  HA V I G A T 1 0 N  S Y S T E M 

INPUT  TIME  INCREMENT 

I N  P  U  T  N  A  V I G  A  T 1 0  N  C  0  H  S  T  A  N  T  -  5" 

INPUT  MISSILE  VELOCITY 

INPUT  TARGET  VELOCITY 

INPUT  TARGET  ACCELERATION 

INPUT  INITIAL  MISSILE  ACCELERATION 

INPUT  INITIAL  RANGE 

INPUT  BETA  TARGET 

INPUT  BETA  MISSILE 

INPUT  THETA  TARGET 

INPUT  THETA  MISSILE 

INPUT  SIGMA 

MISSILE  POSIT  TARGET  POSIT 

XI  Yl  X2  Y2 

D  0  Y  0  U  N  A  N  T  A  P  R  I N  T  0  F  T  H  E  0  U  T  P  U  T  t    O  =  Y  E  S »  1  =  N  0 

M  A  X I M  U  M  M I S  S I L  E  A  C  C  I S  6 5 „ 3 3 5772  4  8     METE  R  S / S  E  C / S E  C 
MISSILE  TIME  OF  FLIGHT  IS    3.2     SEC 
•  INPUT  MINIMUM  VALUE  OF  X 
INPUT  MAXIMUM  VALUE  OF  X 
INPUT  MINIMUM  VALUE  OF  Y 
INPUT  MAXIMUM  VALUE  OF  Y 


H  A  '•/  E  T  H  E  A  X I S  B  E  E  N  D  R  A  N  N :-  O  =  Y  E  S  t  1  =  N  0 
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9    =  48.6°  =  .8483  4ad 
m0 

9    =  -28.4°  =  -.4956  rad 
t0 

aQ      =    77.6°  =  1.3544  rad 

The  problem  was  run  for  k  =  3,4,  and  5. 

Table  3-VIII  is  the  output.   It  can  be  seen  that  the 
effect  of  increasing  the  proportionality  constant  is  to  decrease 
the  maximum  acceleration  required.   Figure  (3-22)  is  a  plot  of 
the  trajectories. 

H.   DESIGN  EXAMPLE  (GUIDANCE  LAW  SELECTION) 

From  the  examples  given  in  this  chapter,  it  can  be  seen  that 
a  missile  designed  to  encounter  a  highly  maneuverable  target, 
such  as  a  fighter,  requires  a  proportional  guidance  law  to  limit 
the  maximum  acceleration  required  of  the  missile.   To  select 
the  proportionality  constant,  it  was  assumed  that  the  threat 
could  maintain  a  constant  7g  turn  at  M,  =  1.5  and  an  altitude 
of  10,000  feet. 

Three  cases  were  investigated,  (1)  A  head-on  encounter  with 
the  target  initiating  a  turn  at  10,000  meters  range,   (2)  A 
crossing  encounter  in  which  the  target  turns  into  the  missile  at 
10,000  meters  range,  and   (3)  An  oblique,  closing  encounter  in 
which  the  target  turns  into  the  missile.   The  scenario  and 
computer  outputs  are  shown  in  Figures  (3-23),  (3-24),  and  (3-25). 
From  this  analysis  the  crossing  encounter  requires  the  largest 
acceleration  (263.2  m/sec/sec) . 

For  the  crossing  scenario  then,  the  missile  speed  was  varied 
and  the  results  indicated  that  as  the  speed  increased  the  maximum 
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V 
m 


V. 


V   = 
m 

V.  = 


821  m/sec  (M  =  2.5) 
492  m/sec  (M  =  1.5) 


rQ  =  10,000  m 


m. 


=  0 


=  180 


7  g  turn 


m, 


=  0 


=  180 


ao  = 

k  = 


THIS  PRO GRAM  FINDS  THE  MAXIMUM 

MISSILE  ACCELERATION  FOR  fl 

P  R  0  P  0  R  T 1 0  N  fl  L    N  fl  V  I G  fl  T 1 0  N    S  Y  S  T  E  M 


INPL 
IN  PI 
IN  PI. 

I  MP  I 
INPl 
INPI. 
I  HP  I 
INPI. 
INPl 
INPI. 
INPI 
INPI 


TI 
Nfl 
MI 
Tfl 
Tfl 
IN 
IN 
BE 
BE 
TH 
TH 
SI 
ISS 


ME 
V I 

RG 
RG 

IT 
IT 
Tfl 
Tfl 
ET 
ET 
GM 
IL 


INC 
GAT  I 

ILE 

ET  V 

ET  fl 

IAL 

IAL 
TAR 
M I S 

fl  Tfl 

fl  MI 

fl 

E  PO 


EMENT 

H  CONSTANT 

ELOCITY 

LOCITY 

CELERATION 

ISS ILE  ACCEL 


:rati 


ET 

ILE 

GET 

SILE 


SI 


XI  Yl 

DO  YOU  WANT  fl  PRINT  OF 


HE  OUTPU" 


TARGET  POSIT 


0=YESj I=NC 


M  fl  X  I M  U  M  M I S  SILE  fl  C  C  I S  126.1911  9  9  4     METE  F:  S  •  S  E  C  ■•  S  E  i 
MISSILE  TIME  OF  FLIGHT  IS    9       SEC 
INPUT  MINIMUM  VALUE  OF  X 

Figure  (3-23).   Head-on  scenario. 
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r«  = 


V.  = 


V   = 


m 


10,000  m 

492  m/sec 

821  m/sec 
o 


m, 


=  -120 


=  -46.2 


fc0   = 

0 

= 

73.8 

mo 

0 

120° 

k     = 

4 

™t-:'T^UI:?RHM  FINDS  THE  MAX  I  MUM 
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Figure  (3-24)  .   Crossing  scenario. 
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Figure  (3-25).   Oblique  scenario. 
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TABLE  3-IX 
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acceleration  decreased.   The  results  for  M  =  2.0  and  M^  =  3.0 
are  shown  in  Table  3-IX.   This  result  indicates  the  desirablity 
of  retaining  the  missile  Mach  number  originally  selected. 

For  the  crossing  case  and  a  missile  speed  of  M^  =  3.0,  the 
proportionality  constant  was  then  varied  from  k  =  2  to  k  =  6 . 
The  results  are  plotted  in  Figure  (3-26).   If  the  maximum 
acceleration  is  limited  to  31  g's  (3a   +  10),  the  required  pro- 
portionality constant  is  k  =  3.75.   This  is  well  within  the 
desirable  range  of  2  -  6  indicated  earlier. 

From  this  analysis  the  required  performance  objectives  are: 

M      =    3.0 
m 

k  =  3.75 


(a)     =  31  g's 
m  max 
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IV.   SIZING  THE  DIAMETER 

The  missile  diameter  is  determined  by  one  of  three  driving 
factors.   For  relatively  short  range  missiles  the  diameter  will 
be  fixed  by  either  the  warhead  or  the  seeker  requirements.   As 
might  be  expected,  for  longer  range  missiles  the  diameter  will 
more  likely  be  fixed  by  either  the  warhead  or  the  seeker  require- 
ments.  As  might  be  expected,  for  longer  range  missiles  the 
diameter  will  more  likely  be  fixed  by  the  propulsion  require- 
ments in  order  to  prevent  excessive  propulsion  system  lengths. 
An  initial  estimate  of  the  missile  diameter  must  be  made  at 
this  point  in  order  to  proceed  with  the  design.   The  initial 
seeker  requirement  can  be  determined  from  a  knowledge  of  the 
lock-on  range  requirement  found  in  Chapter  2.   The  warhead 
necessary  to  inflict  a  "kill"  can  also  be  estimated  from  infor- 
mation about  the  target  and  characteristic  explosives.   The 
propulsion  requirement  cannot  be  determined  because  of  the  lack 
of  any  aerodynamic  drag  or  weight  information  at  this  point.   For 
this  reason  the  missile  diameter  will  be  now  sized  for  seeker 
or  warhead  requirements.   The  missile  propulsion  requirements 
will  be  determined  later  in  the  design  process,  and  it  may  be 
necessary  at  that  time  to  resize  the  missile  to  meet  these  pro- 
pulsion requirements. 

Selection  of  the  type  of  seeker  depends  upon  the  operational 
arena    of  the  missile.   The  seeker  of  a  shoulder  fired,  battle- 
field missile  would  not  be  the  optimum  seeker  of  a  shipboard 
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missile  where  antenna  and  component  sizes  are  not  limiting 
factors.   All  types  of  guidance  use  some  portion  of  the  electro- 
magnetic spectrum.   The  three  primary  areas  of  use  are  the 
electro-optical,  infrared,  and  radio  frequencies.   The  milli- 
meter wave  section  of  the  spectrum  is  also  of  current  interest 
in  the  design  of  missiles  due  to  small  component  size  and  will 
also  be  discussed.   The  following  table  lists  some  of  the  major 
advantages  and  disadvantages  of  the  three. 


Optical 


Infrared 


RF 


Advantages 

Target  resolution  (de- 
tail) 

Real  time  information 
Three  dimension  effect 

Improved  resolution 
over  RF 

Longest  range 

Least  absorption  and 

attenuation 


Disadvantages 

Bad  weather  degrades 
Night  use  degrades 


Attenuation  due  to 
aerosols  and  atmosphere 

Larger  components 


A.   THE  RADAR  RANGE  EQUATION 

An  omnidirectional  antenna  is  one  that  radiates  power  in  all 
directions  equally.   If  the  power  radiated  by  an  antenna  is  P  , 
the  power  density  at  a  distance  R.  from  the  source  is  given  by, 


Power  density  =  P  /  ( 4tt  R   ) 


4tt  R   =  area  of  a  sphere  of  radius  R. 


Since  antennas  are  normally  directive  instead  of  omnidirectional, 
most  of  the  power  is  radiated  in  a  particular  direction.   The 
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gain,  G  ,  is  a  measure  of  the  increased  power  from  a  directive 
antenna  as  compared  to  an  omnidirectional  antenna.   Therefore  the 
power  density  from  a  directive  antenna  can  be  expressed  as, 

P   G 
t   t 
Power  density  =   *■ 

4TT  Rt 

This  is  the  power  density  which  arrives  at  the  target.   The  tar- 
get intercepts  a  portion  of  this  energy  and  reradiates  it  in  the 
opposite  direction.   The  radar  cross  section,  a,  is  a  measure  of 
the  effective  area  of  the  target.   The  power  radiated  by  the 
target  is  P   .  , 

P  G  a 

echo     .   _  2 
t 

This  energy  propagates  as  if  it  were  radiated  by  an  omni- 
directional antenna.   Therefore  if  the  receiving  antenna  is  a 
distance,  R  ,  away  the  power  density  at  the  receiver  is 

P   G   a 
t   t 
(Power  density)   =   ^ 5— 

r     (4tt  RtZ)  (4tt  RrZ) 

If  the  energy  is  intercepted  by  the  receiving  antenna,  which  has 
an  effective  area  as  seen  by  the  returned  energy  of  A  ;  then  the 
power  received  by  the  radar  ,  P  ,  is 
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P.  G.  a  A 

pr=  — H — *-*-  «« 

(4tt  r  z)  (4tt  R  *) 

I—  X. 


This  is  the  simplest  form  of  the  radar  equation  and  can  be  used 
to  determine  the  size  of  antenna  required. 

B.   ACTIVE  RADAR  HOMING 

Active  homing  is  the  method  of  missile  guidance  in  which  the 
radar  transmitter  and  receiver  are  located  on-board  the  missile. 
In  this  case  the  same  antenna  is  used  for  both  transmitting  and 
receiving.   The  radar  equation  then  becomes, 


Pfc  Gt  g  Afc 
r      (4tt  R^) 


where,  R.  =  R  =  R  and  A  =  A.  . 
t    r  r    t 

The  minimum  power  for  which  the  target  can  be  detected,  P    ,  is 
r  3  nan 

a  function  of  many  variables.   A  full  development  of  this  term 
can  be  found  in  reference  (8) . 

P  .   =  k  Tn  B   F   (~) 

mm      0   n   n  Nn'  . 

0  min 

Boltzmans  constant  is  k  =  1.38  x  10   "  joule/  k.   The  value  of 

-21 

kT   at  room  temperature  is  4  x  10   '  watt/cps  of  bandwidth.   The 

bandwidth,  B  ,  noise  figure,  F  ,  and  minimum  signal  to  noise 

ratio,  (S_/N  )     are  all  functions  of  the  receiver.   Typical 
min 
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values  are  listed  below. 


B   =  1  MHz 

n 

F   =  7.5  db  =  2.37  (for  crystal  mixer) 

(Sq/Nq)     =  14.7  db  =  5.43  (for  probability  of  detection, 
min 

PD  =  0.9  and  probability  of  false  alarm, 

P,.   =  1/15  minutes) 

The  above  values  give  P  .   =  5.15  x  10~   watts.   This  is  the 

mm 

value  which  will  be  used  throughout  this  section. 

From  antenna  theory  the  gain  is  related  to  the  effective 
antenna  area  by, 


n     -        4tT  & 


The  maximum  radar  range  can  then  be  shown  to  be 


R 
max 


Pt  a  At 

47T  X2  P   . 

mm 


2  n   V4 


(1) 


Equation  (1)  for  radar  range  does  not  include  any  system  losses. 
It  also  does  not  include  the  statistical  nature  of  several  of 
the  parameters.   Because  of  these  assumptions  the  actual  range 
of  a  radar  may  be  as  small  as  one-half  of  what  the  radar  range 
equation  predicts  for  laboratory  conditions.   For  this  reason, 
twice  the  required  range  should  be  used  when  using  the  above 
equation. 
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1.   Example 

The  AN/APQ-15  3  is  the  airborne  attack  radar  system  used 
on  the  F-5E  aircraft.   The  following  parameters  apply  to  this 
radar. 


fQ  =  8-10  GHz 

A 


/•f   =  3  x  10   m/sec 

c/  0  9 

u    9  x  10  /sec 


=  .0333  m 


c  =  speed  of  electromagnetic  propagation 

o 

=  3  x  10   m/sec 

P   =  80  kW 

Antenna  =  Parabolic  dish  30.5  x  40.6  cm 

Afc  =  0.12383  m2 

for  a  target  of  1  square  meter  of  radar  cross-section, 


R    =  36 . 16  km  =  19 . 5  nmiles 
max 

For  the  case  of  an  active  homing  radar,  the  size  of  the 

transmitting  and  receiving  antenna  is  the  parameter  of  interest 

The  antenna  size  may  very  well  drive  the  design  diameter  of  the 

missile.   The  antenna  diameter  can  be  expressed  as  follows. 

1/2 


Let 


then 


At " 


C  = 


4tt  A2  P  .   R 

mm   max 


9      o 


4tt 

X2 

P      . 

min 

pt 

a 

1/2 


(2) 


ndt  : 

t      4       max 


4C  R 


max 


therefore.   d   =  2  R 

'    t      max 
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Equation  (2)  has  been  plotted  in  Figure  (4-1)  for  various 
values  of  the  transmitter  frequency.   From  this  plot  and  a  know- 
ledge of  the  maximum  lock-on  range  required,  the  antenna  size 
can  be  determined. 

From  equation  (2)  and  Figure  (4-1)  there  are  two  obvious 
ways  to  decrease  the  antenna  size  required.   (1)  Increasing  fre- 
quency is  the  best  way  to  reduce  antenna  and  electronic  component 
sizes.   The  current  trend  is  toward  higher  frequencies.   (Milli- 
meter waves.)   One  problem  is  that  the  equations  developed  in 
this  section  do  not  include  atmospheric  attenuation.   For  fre- 
quencies above  about  30  GHz  the  absoprtion  due  to  atmospheric 
gases  increases.   This  is  shown  on  Figure  (4-2).   As  indicated 
on  this  figure  there  are  "windows"  where  the  attenuation  is  less. 
These  "windows"  occur  at  frequencies  of  34  GHz,  94  GHz,  140  GHz, 
and  220  GHz.   These  are  the  frequencies  where  most  of  the  cur- 
rent research  and  development  is  going  on.   As  the  frequency  in- 
creases, the  wavelength  approaches  the  size  of  rain  droplets. 
For  this  reason,  radar  performance  is  greatly  reduced  in  inclement 
weather.   (2)  Increasing  transmitter  power  will  also  decrease 
the  size  of  antenna  necessary.   The  limiting  factor  in  this  area 
is  the  lack  of  high  power  sources.   In  the  millimeter  range  the 
available  power  from  current  traveling  wave  tubes  is  50-100 
watts.   Increasing  power  is  obviously  confined  to  size  and  weight 
limitations  of  missile  components. 
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Figure  (4-2).   Atmospheric  Absorption  [8] 
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C.   SEMI-ACTIVE  HOMING 

The  advantage  of  semi-active  homing  is  obvious  when  the  radar 
range  equation  is  investigated.   From  equation  (1) , 

9   0    P   G.  a   A 

R  2P  2  -     t       t  r 

Rt  Rr   "  77T2Z 

(4tt)  P  . 
mm 

In  the  above  equation  the  missile  range  from  the  target  is  R  . 
The  transmitting  and  receiving  antennas  are  at  different  ranges 
and  have  different  characteristics  in  this  case.   As  before, 

4tt  a. 


Gt  =  ~ 


therefore, 


0   ~   P.  A,  A  a 

R  2r  2  =    t  t  r 

t      r      4 it  X2  P  . 

min 


The  main  advantage  is  in  the  transmitter  characteristics.   Since 
the  transmitter  is  not  located  in  the  missile,  it  is  not  normally 
limited  in  size  and  weight  requirements.   In  the  above  equation 
if  a  transmitter  power  and  standoff  range,  R  ,  is  chosen  the 
receiving  antenna  can  be  sized  for  a  maximum  homing  range,  R    , 


of  the  missile. 
1.   Example 


Pt  =  100  km 

R  =  100  nmiles 

f   =  10  GHz 

At  =  4  m 

P     =  5.15  x  10"14  W 
min 
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R     (nmiles)  d   (in) 

max  r 


10  5.8 

20  11.6 

50  29.1 

It  can  be  seen  from  comparing  these  numbers  to  those  of  Figure 
(4-1)  that  the  required  antenna  size  is  less  than  one  half  of 
that  required  for  an  active  homing  radar  of  the  same  frequency. 

D.  DESIGN  EXAMPLE  (ANTENNA  SIZING) 

An  active  radar  was  assumed  in  Chapter  2  to  decrease  the 
required  missile  range.   The  lock-on  range  was  10,000  m  or 
5.4  nmiles.   As  stated  earlier,  twice  this  number  should  be  used 
for  determining  antenna  size.   From  Figure  (1)  for  a  range  of 
10.8  nmiles,  and  a  transmitter  power  of  10  kW  at  f   =  20  GHz, 
the  required  antenna  size  is  d  =  10  inches. 

E.  INFRARED  SEEKERS 

In  the  design  of  missile  seekers  two  parameters  of  primary 
importance  are  range  and  size.   The  idealized  range  for  an  infra- 
red tracker  relates  these  two  factors.   The  idealized  range  is 
the  range  at  which  the  signal-to-noise  ratio  is  unity  and  is 
given  by, 

1/2 


R0    = 


D*T      TT„         D    2    J 
a      IR  a 


4     y/bTH 


(1) 
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The  derivation  of  equation  (1)  and  its  use  are  the  subjects  of 
this  section.  Some  references  (9)  may  give  the  above  equation 
in  terms  of  the  Noise  Equivalent  Intensity,  NEI. 


where 


Ro  ■ 


NEI    = 


1/2 


NEI 


/Af~T 


D*    A      T      TIR 
a      a 


1.      Planck's    Law 

The  radiant  emittance  of  a  body  is  a  measure  of  the  radiant 
power  per  unit  area  emitted  from  the  surface. 


P      ,   2 

W  =  —  watt/cm 


(2) 


The  spectral  radiant  emittance  is  the  radiant  emittance  per  unit 
wavelength  interval, 


103 


\      =  TT  watt/cm  I-1 

—  6 
y   =  micron  =  10    meters. 

Planck's  law  gives  the  blackbody  spectral  radiant  emittance  as 
a  function  of  wavelength  and  temperature, 

(W  )    =  2_JL_h_cJl  1 () 

K    X;BB      .5     exp(hc/AkT)-l  lJ; 

A 

-34 
h  =  Planck's  constant  =  6.6238  x  10     Joule-sec 

g 

c  =  speed  of  light  =  3  x  10   m/sec 

X  =  wavelength 

k  =  Boltzmann's  constant  =  1.38  x  10~    Joule/  K 

T  =  Absolute  temperature,   K 

Figure  (4-3)  shows  equation  (3)  for  various  absolute  tempera- 
tures.  As  can  be  seen  the  wavelength  at  which  maximum  radiant 
emittance  occurs  varies  with  temperature.   This  maximum  occurs 

at  a  wavelength  given  by  Wien's  displacement  law,  X 

r    3       J  r  max 

X  T   -  2897. 8y  °K 

max 

2.   Emissivity 

Actual  bodies  do  not  emit  radiation  according  to  Planck's 
law.   A  more  typical  plot  of  radiant  emittance  is  shown  in 
Figure  (4-4).   Spectral  emissivity,  E,,  is  defined  as  the  ratio 
of  the  actual  spectral  radiant  emittance  to  the  blackbody  spectral 
radiant  emittance, 
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Figure  (4-3) .   Spectral  radiant  emittance. 
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Figure  (4-4).   Actual  spectral  radiant  emittance. 


As  shown  in  the  above  figure,  £,,  may  vary  with  wavelength.   A 
grey  body  is  defined  as  one  which  has  constant  spectral  emissivity, 

e  =  e,  =  constant 

IR  systems  normally  use  filters  to  limit  the  accepted 
radiation  to  a  specific  wavelength  band.   The  radiant  emittance 
of  a  body  between  wavelengths  X,  and  X^   becomes, 


A, 


W 


■   / 


e  (W, )    dX 


X 


BB 


In  Figure  (4-5)  the  surface  at  the  orgin  emits  a  total 
energy  WA  into  a  hemisphere  normal  to  A. 

The  radiance  is  defined  as  the  radiant  power  per  unit 
solid  angle  per  unit  projected  area, 

Radiance  =  N  = 


d    P 


cos98A3ft    cos6  9ft 


dA 
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Solid  angle  ft. 


-^•Surface  area 
A,  e,    T 


Figure  (4-5).   Radiant  emittance  of  a  body. 


From  equation  (2) 


N  = 


cos9  9ft 


watts/steradian  cm 


The  radiant  intensity  is  defined  as  the  radiant  power  per  unit 
solid  angle  from  a  point  source. 

Radiant  intensity  =  J  =  -^r  watt/steradian 

The  following  is  a  summary  of  the  definitions  of  radiant  energy 
quantities, 

,   2 

W  =  radiant  emittance,  watts/cm 

J  =  radiant  intensity,  watts/steradian 

2 
N  =  radiance,  watts/steradian-cm 
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From  the  definitions  the  relationship  between  W,  J  and  N  are, 

W  =  ttN  =  tt  J /A 

From  the  above  definitions  and  Figure  (4-5)  the  energy  into  solid 
angle,  ft,  is, 

W 

NA  ft,  cos  9  =  —ft,  A  cos  0 
1         tt   1 

The  energy  into  solid  angle  ft_  is  given  by, 

W 

NAft.  =  -  ft.  A 
2     TT   2 

3 .   Energy  into  a  Hemisphere 

From  the  definition  of  radiance  and  the  radiant  emittance, 
the  radiance  in  terms  of  the  radiant  emittance  can  be  found, 


N  ■ 


W  = 


N  = 


32P 

cos98A8ft 

8P 

dA 

1 

9W 

cos9 

8ft 

rsin9d<j> 


Figure  (4-6).   Energy  emitted  into  a  hemisphere. 
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From  Figure  (4-6)  the  incremental  solid  angle  is 


dfi 


r   sin  6djr   d8_ 


=   sin  9  d  8  d  <J> 


If  the  above  surface,  A  ,  is  considered  a  Lambertian  surface, 
the  radiance,  N,  is  independent  of  the  direction  of  radiation 

dW  =  cos  6  N  dft 


dW  =    cos  0  N    sin  8d8d|i 

The  total  radiant  emittance  into  a  hemisphere  above  the  surface 
is  then, 

277       -77/2 

W  =       J         f       N   cos  9  sin  9  d  9  d  <£ 
o        o 


W  =    2tt   N 


1  .     2 

2  Sin 


77/2 


W    =     77    N 

4.       Targets 

Infrared  targets  include  a  wide  variety  of  radiation 
sources.   The  radiance  of  most  bodies  can  be  divided  into  that 
due  to  self-emission  and  that  due  to  reflection  of  incident 
radiation, 


N  =  N   +  N 
e    r 


The  relative  magnitude  of  these  contributions  depends  on  a  number 

of  factors  and  varies  from  target  to  target  and  operating 

environment. 
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a.  Self-emission. 
Self-emission,  also  referred  to  as  thermal  emission, 

depends  primarily  on  the  temperature  of  the  body  and  the  emis- 
sivity.   The  most  often  used  appraoch  is  to  consider  the  body  as 
a  grey  body  which  emits  radiation  according  to  the  Stefan- 
Bo  It zmann  law. 

4 

N   =  watts/steradian-cm 

e     tt 

where         a  =  5.67  x  10~12  watts/cm2 (°K) 4 

This  term  is  the  total  radiance  (over  all  wavelengths)  and  is 
not  the  same  as  used  previously. 

b.  Reflection. 
Radiance  due  to  reflectance  depends  on  the  illumi- 
nating source.   This  source  may  be  the  sun,  active,  or  semi- 
active  sources.   It  is  obvious  that  at  night  for. a  passive 
infrared  system,  the  radiance  due  to  reflection  is  not  a  con- 
tributing factor.   For  this  reason  only  the  radiance  due  to  self- 
emission  is  considered  in  this  section. 

5.   Target  Temperature 

Since  the  self-emittance  of  a  target  depends  on  the 
temperature  of  the  target,  a  method  for  determining  this  tem- 
perature is  needed.   The  temperature  of  an  aerial  target  varies 
depending  on  the  aspect  of  the  target.   The  propulsion  system 
has  hot  surfaces  such  as  the  nozzle  and  exhaust  plumes.   There 
may  also  be  hot  surfaces  due  to  aerodynamic  heating  and/or 
solar  radiation. 
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Air  breathing  engines  normally  have  exhaust  plumes 
ranging  from  600  to  1000  K.   Rockets  typically  have  much  hotter 
plumes.   The  flame  temperatures  for  liquid  propellants  range 
from  2500  to  7500  K.   Solid  propellants  flame  temperatures  range 
from  1700  to  3500  K.   The  plume  temperature  can  be  estimated 
from  the  relation. 


Tflame   =  JO  =  i  +  Izl  Me2 
T  T  2 

plume     e 


where    T   =  stagnation  temperature 
T   =  static  temperature 


M  =  Mach  number  at  the  nozzle  exit 
e 


a.   Example  I 

For  a  flame  temperature  of  2700  K  and  an  exit  Mach 
number  of  3.0  the  plume  temperature  can  be  found, 


Tflame   =  1270.6°K 


plume    1+I-1M  2 
2   e 


For  many  missile  encounters  the  exhaust  plume  may 
be  shielded  from  the  infrared  sensor.   For  a  head-on  encounter 
the  temperature  of  interest  is  the  skin  temperature  of  the  target 
This  temperature  is  due  to  aerodynamic  heating  and  is  a  function 
of  the  target  speed  and  the  target  material.   One  approach  to 
finding  this  temperature  is  through  the  use  of  the  recovery 
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factor,  which  requires  some  knowledge  of  the  material  of  the 
target.  The  recovery  factor,  r,  of  a  material  is  defined  as 
follows : 


T         -  T 
surface     ambient 


T  -  T 

stagnation    ambient 


The  skin  temperature  of  the  target  then  becomes, 


T         =T         +  r  (T      -T) 
surface    ambient       stag    amb 


The  stagnation  temperature  is  found  from  the  relation, 


T  .    -  T   .   (1  +  ^L  M2) 
stag     amb        2 


The  Mach  number,  M,  is  that  of  the  target  and  the  specific  heat 
ratio,  y,  is  for  air. 
b.   Example  II 

A  target  flying  at  M  =  2.5  where  the  ambient  temper- 
ature is  300  K  has  a  recovery  factor  of  0.75. 

T  .    =  575°K 
stag 

T    .       581°K 
surface 

This  is  the  temperature  used,  along  with  the  emissivity  of  the 
target,  to  find  the  radiant  emittance  of  the  target  from 
equation  (4) . 
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6 .   Simple  IR  System 


Target 
A  /  e,  T 


Detector 


filter 


I 


chopper 


preamp 


Figure  (4-7) .   Simple  IR  system. 

Figure  (4-7),  above,  shows  a  simple  IR  system  and  a 
target  at  a  range  R.   If  the  system  is  sensitive  to  radiation 
in  the  3  to  5  micron  region  the  radiant  emittance  becomes, 


X2=5 


W  = 


X  =3 


f  e  (WJ    dX 

\  A  BB 


W-  c  / 

3 


27rhc' 


X 


5    exp(hc/ 


1 1 

c/  kT)-l 


dX 


The  above  integral  is  best  evaluated  on  the  computer.   If  the 
ambient  temperature  and  target  speed  is  known,  the  skin  temper- 
ature of  the  target  can  be  determined. 
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The  radiance  from  the  target  then  becomes, 


W 
N  =   - 

7T 


If  Figure  (4-7) ,  the  solid  angle  of  the  aperture  as  seen  from 

the  target  is,    0   _a 

R 

The  power  seen  at  the  detector  surface  is  then, 
IR  Power  =  N  Q   A 

The  above  formula  assumes  no  attenuation  by  the  atmosphere  or 
the  IR  system  optics.   This  attenuation  is  significant  in  actual 
IR  systems.   These  factors  are  normally  accounted  for  through 
the  use  of  atmospheric  and  IR  optics  transmission  coefficients. 

T   =  transmission  of  the  atmosphere 
a  r 

TTT5   =   transmission   of   the   IR  optics 

IK 

The  total  power  at  the  detector  then  becomes, 

T   TV   A^  A   W 
P  =  g   I*  t  a (5) 

ttR 

7 .   Detectors 

Detectors  are  devices  which  are  radiation  transducers. 
It's  purpose  is  to  change  the  incoming  radiant  power  to  an 
electrical  signal,  which  can  then  be  amplified.   Detectors  can 
be  divided  into  two  main  categories.   (1)  Thermal  detectors  - 
The  responsive  element  of  a  thermal  detector  is  sensitive  to 
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temperature  changes  brought  about  by  the  incident  radiation. 
(2)  Photodectectors  -  Responsive  elements  of  photodetctors  are 
sensitive  to  the  number  of  incident  photons. 

Detectors  also  are  made  up  of  windows,  apertures  and 
Dewar  flasks.   The  window  restricts  the  bandwidth  to  which  the 
detector  is  sensitive.   The  aperture  may  limit  the  field  of  vie- 
in  order  to  limit  photon  noise.   The  Dewar  flask  cools  the 
detector  which  improves  the  detectivity. 

Detectivity  of  a  detector  is  defined  as, 

_    signal/noise    S/N  tc. 

D    =  -. - —  — — —  \b) 

input  power      P 

The  specific  detectivity  is 

D*  =  D  [  Af  Ad]  1/2  (7) 

Af  =  Bandwidth 


A,  =  Sensitive  area  of  detector 
d 


For  a  tracking  system  the  bandwidth  is  that  of  the  preamplifier 
in  Figure  (4-7) .   The  input  to  the  preamplifier  is  proportional 
to  the  incoming  IR  energy,  which  has  been  modulated  to  give  tar- 
get resolution  from  the  background  and  provide  line-of-sight 
information. 

A  simple  chopper  is  shown  in  Figure  (4-8)  .  It  consists 
of  an  opaque  material  which  has  a  wedge  cut-out  of  angle  a.  The 
rotation  causes  the  input  from  a  point  source  to  be  modulated, 
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Figure  (4-8) .   Simple  Chopper. 


while  that  of  the  background  is  not.   The  input  signal  to  the 

preamplifier  would  look  like  Figure  (4-9) .   The  frequency  content 

of  the  signal  in  Figure  (4-9)  can  be  found  from  a  Fourier  Analysis 

If  the  pulses  are  assumed  to  be  sinusoidal  of  period  T,  the 

optimum  bandwidth  is,        -> 

Af    -I 


T    = 


W 


TTCtW 


.Target^ 


T     = 


2jr 

w 


Figure    (4-9) .      Preamp    input 
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The  specific  detectivity  is  characteristic  of  the 
detector  used.   Reference  [9]  is  an  excellent  source  of  infor- 
mation on  operational  detectors. 
8 .   Idealized  Range 
Equation  (5)  is, 


P  = 


Ta  TIR  At  Aa  W 
ttR 


From  equations  (6) and  (7)  this  becomes, 


S/N  tt  R' 


[Af  Ad]1/2   Ta  TIR  At  Aa  W 


The  idealized  range,  where  the  signal  to  noise  ratio  is  unity,  is 

„  2    D*  Ta  TIR  At  Aa  W 
0   =    [Af  Ad]  1/2, 


To  simplify  this  equation  the  radiant  intensity  is  given  by, 


J  = 


At  W 

TT 


Since  the  parameter  of  interest  in  missile  design  is  the  aperture 

ttD 


diameters  A  ,  it  is  replaced  with 
a'         v 


so    that, 


Ro  ■ 


D      = 
a 


2       1  1/2 

D*    T      T__     tt    D         J 
a     IR  a 


4  /AfA 


RQ       4/AF-A7 

D*    T       TTTD     TT    J 
a     1K 


-I  1/2 


117 


a.   Example  (Idealized  Range) 

From  example  II  a  target  flying  at  M  =  2.5  has  a  skin 

temperature  of  581  K.   If  this  target  has  a  presented  area  of 

2 

1  m  ,  the  detector  size  needed  to  detect  the  target  at  a  range 

of  10,000  m  can  be  determined.   From  Wien's  displacement  law  the 
maximum  radiation  occurs  at, 


X  =  4.99  y 

max 


If  the  system  is  designed  to  accept  radiation  from  3  to  5  microns, 
and  the  emissivity  of  the  target  is  0.7. 

5        2 

2-rrhc         1  ,, 

,5   exp (hc/XkT) -1 


If  this  equation  is  integrated  on  the  computer,  the  radiant 
emittance  becomes, 

W  =  1410  watts/m 

The  radiant  intensity  becomes, 

A   W 
J  =  — - —  =448.82  watts 

Typical  values  of  the  parameters  in  the  idealized  range  equation 
are  (9) 


T   =  0.75 
a 

TIR  -  0.95 

Af  =  1000  cps 

A,  =  1  cm 
d 

D*  =  1  x  1010  (cps)1/'  cm/watt 


Substituting  these  values  into  the  idealized  range  equation  gives, 


D   =  0.0355  m  =  1.4  in. 
a 


As  can  be  seen  the  size  of  the  IR  seeker  is  relatively  small 
compared  to  other  seekers.   The  above  analysis  is  for  the 
"idealized"  range.   Attenuation  of  IR  radiation  can  be  quite 
high  thereby  increasing  the  seeker  size  required. 

F.   WARHEAD  SIZING 

The  conditional  kill  probability  of  a  missile  is  the  prob- 
ability that  the  target  is  destroyed  given  that  the  warhead  is 
delivered  to  a  point  in  space  and  the  fuze  detonates  the  war- 
head at  a  miss  distance  r. 


Figure  (4-10)  .   Encounter  geometry 
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The  fragment  distribution,  D(cf>),  is  the  number  of  fragments  per 
unit  solid  angle,  Q.   The  total  number  of  fragments  within  the 
cone  is, 

N  =  D(4>)ft 

The  fragment  density,  p,  is  the  number  of  fragments  per  unit  of 
area  normal  to  the  path. 

N    D(4>)fl 
P    A      A 


however, 


and 


p  =  Dji!  =  D(cj>)sin2(f)  (8) 

r        r 
m 


1.   Target  Vulnerability 

The  vulnerability  of  aircraft  or  missile  components  is 
normally  determined  experimentally.   Fragments  of  a  specified 
size  are  fired  at  the  component,  and  the  damage  is  assessed  to 
determine  if  the  fragment  would  cause  a  kill.   If  a  large  num- 
ber of  fragments  are  fired,  the  ratio  of  killing  fragments  to 

hits  can  be  determined.   This  ratio  is  defined  as  the  probability 

N 
that  given  a  hit  a  kill  will  result,   ^       K    ,     .    . . 

3  '   PT,  /TT  =  S5—  •   Assuming  the 

K/H    NH 

distribution  of  hits  is  uniform  over  the  target, 

A 

P     =   — 
K/H  "A  (9) 

P 

120 


A   =  presented  area  of  target 
A  =  vulnerable  area  of  target 

As  would  be  expected,  PK/H#  depends  greatly  upon  the 
encounter  geometry.   It  will  be  dependent  upon  the  aspect  of 
the  aircraft,  and  also  depends  upon  the  type  of  kill  specified. 
If  the  target  is  assumed  to  be  a  spherical  target,  the  prob- 
ability of  kill  given  a  hit  can  be  assumed  constant. 
2.   Conditional  Kill  Probability 

The  fragment  density  is  given  by  equation  (8).   From 
this  density  expression,  the  average  number  of  hits,  a,  on  the 
vulnerable  area  of  the  target  is, 

2 

_  D  (<j>)  sin  (ft 

2  v 

rm 

It  is  customary  to  assume  that  the  distribution  of  hits  on  the 
presented  area  follows  a  Poisson  distribution.   The  conditional 
kill  probability  then  becomes, 

PD  =  1  -  e"a  (10) 


PD  =  1  -  exp 


/  x  .  2 

-D  (<ft)  sin  (ft  Av 

2 

r 

m 


The  above  expression  depends  upon  the  fragment  distribution, 
D((ft).   If  the  warhead  casing  is  scored  such  that  it  produces  N 
fragments  of  uniform  size  and  mass,  m,  the  problem  is  simplified 
by  formulating  an  alternate  expression  for  a. 
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£w  +  2rm  tan  6 


w 


Figure  (4-11) .   Static  encounter 


The  warhead-target  encounter  geometry  is  shown  in 
Figure  (4-11) .   As  shown  the  area  of  the  fragment  ring  depends 
on  the  miss  distance,  r  .   The  area  of  the  ring  is  given  by,  A 


m 


A£  =  2tt   r  (£.   +  2r   tan  6) 
f        m  w     m 


If  the  N  fragments  are  distributed  evenly  in  Af/  the  number  of 


fragments  per  unit  area  is 


P  = 


N  =  N 

A£  *  2  it  r  U   +  2r   tan  6) 
f       m  w     m 


The  average  number  of  hits  on  a  vulnerable  component  is  then, 


a  =  p  A 


v 
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From  equation  (9),  A   =  P  .      A 

V      K/  rl   p 


Therefore,    a  =  p  p__  .„   A 

K/H   p 


a  = 


N  PK/H  % 


2tt  r  U   +  2r   tan  3) 
m  w     m 


If  in  Figure  (4-11) ,  the  average  target  width  into  the  paper  is 
given  by  W,  the  presented  area  is, 


A^  =  w(£   +  2r   tan  3 ) 
p      w     m 


The  average  number  of  hits  is  then, 


N  p     W 

ar  =  ->*'  (ID 

I     2tt  r 

m 


The  distance  at  which  the  target  just  fills  the  fragment  ring 
is  the  critical  miss  distance,  and  the  maximum  distance  for 
which  a   applies.   This  is  found  by  setting 

L  =  £   +  2r   tan  3 
w     c 

From  which 

L  -  % 

w 
r   = 


c    2  tan  3 


If  the  miss  distance   in  Figure  (4-11)  is  such  that  the  entire 
target  is  always  presented  to  the  fragment  ring;  i.e.,  r   >  r  , 
the  presented  area  becomes , 


A   =  W  L 
P 
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The  average  number  of  hits  then  becomes,  a 


N  PK/H  W  L 

(12) 


II    2TrU   +  2r   tan  6) 
m  w     m 


3.   Sizing  the  Warhead  Radius 

The  parameter  of  interest  in  this  chapter  is  the  diameter 
of  the  warhead  required  to  achieve  a  specified  kill  probability. 
The  development  thus  far  is  the  conditional  probability  of  kill, 
P  .   It  has  been  assumed  that  the  guidance  system  delivers  the 
warhead  to  the  point  of  interest,  and  the  fuze  detonates  the 
warhead  at  this  point.   Since  only  the  conditional  probability 
of  kill  is  determined,  the  purpose  of  this  section  will  be  to 
maximize  P   or  to  find  the  warhead  diameter  which  sets  P   =  1. 

From  the  threat  to  be  encountered,  the  target  presented 
area,  A  ,  and  the  vulnerability,  P  ,„,  can  be  determined.   Also 
from  an  alaysis  of  the  threat,  the  siz  e  and  impact  velocity  of 
the  fragments  necessary  to  kill  the  target  can  be  determined. 
The  initial  velocity  required  to  obtain  the  impact  velocity  is 
a  function  of  the  explosive  used  and  the  charge  to  mass  ratio, 
C/M. 

C  _  Mass  of  explosive/unit  length 

M   Mass  of  warhead  casing/unit  length 

The  initial  velocity  is  given  by  Gurney ' s  equation, 

1    1/2 


V.  =  J2E 
1    v 


C/M 


1+C/2M 


(13) 
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Gurney ' s  constant,  v2E,  for  various  explosives  are  given  below  [11 


Explosive 

Density,  kg/m 

TNT 

1590 

RDX 

1650 

HMX 

1840 

PETN 

1730 

Tetryl 

1620 

Composition  B 

1680 

Octol 

1800 

n/2E, m/sec 
2316.5 
2834.6 
3118.1 
2834.6 
2500.0 
2682.2 
2895.6 


From  equation  (13) ,  the  charge  to  mass  ratio  necessary  to  attain 
the  specified  initial  velocity  can  be  determined. 


c   vi  /2E 


M        2 

1-V.  /2(2E) 

l 


Figure  (4-12) .   Warhead 
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From  Figure  (4-12) ,  C/M  can  be  expressed  as, 


C     ^e2  pe 


where 


M    ,  _.  2       2. 

^rw   ~  ^e  )pc 


p   =  explosive  density 
p   =  casing  density 


2 

2       re  ^e 

M  "  (r  +t)2  -  r  2    pc 
e        e 


2 

C  _  pe 

M    2r  t  +  t2   pc 
e 


If  t  <  <  r 

e 


2 

r    p     p   r 
e    e    e   e 


M    2r  t  p     p   2t 
e    c    c 


C   pe   (r„  "  tJ 


MP       2t 
c 


(14) 


From  equation  (14)  the  casing  thickness  in  terms  of  the  war- 
head radius  can  be  determined. 


r 
t  = 


2  S  £n  +  i 


M  Pe 
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With  the  casing  thickness  fixed  as  a  function  of  warhead  radius 
the  number  of  fragments,  N,  can  be  determined.   In  order  to 
achieve  the  desired  velocities  of  fragments,  it  is  essential 
to  have  enough  length  for  the  diameter.   An  acceptable  length 
to  diameter  ratio  for  a  cylindrical  warhead  is  from  two  to  three 
Most  air-to-air  missiles  have  a  length  to  diameter  ratio  of  2.5. 
For  this  analysis  a  value  of  2.5  is  used. 
The  warhead  casing  volume  is  then  given  by,  V  , 

V  =    27rr    t£ 
c  WW 

2 

V  =    lOirtr 
c  w 


The  mass  of  the  casing  is,  m  , 

3    '   c 


2 

m  =  lOiTtr   P 
c        w   c 


The  total  number  of  fragments,  N,  is  obtained  by  dividing  the 
case  mass  by  the  individual  fragment  mass. 

m    lOTTtr   p 

N  =  _£  =  JtL-S. 

m       m 

From  the  equation  (11)  or  (12)  the  average  number  of  hits  can 

be  determined  as  a  function  of  warhead  radius 

5tr  2P  PWUW 
_     w   c  K/H 

ai  ~      mr 
m 

5tr  2p  PV/UWL 
w   c  K/H 


11   mr  (i   +2r  tanB) 

ni   w    m 
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The  conditional  probability  of  kill  can  be  determined  now  for  a 
given  warhead  radius  using  equation  (10) . 

a.   Design  Example  (Effect  of  radius  on  P  ) 

If  the  above  equations  are  programmed  for  the  con- 
ditions listed  below,  a  plot  of  the  conditional  probability  of 
kill  versus  warhead  radius  can  be  obtained. 
Threat:   RAM-K 

L  =  19.51  m 
w  =  6.10  m 

PK/H  =  -10 

Fragments : 

m  =  105  grains  =  0.0068  kg 

V.  =  2133.6  m/sec 
l 

p   =  7000  kg/m 
S  =  20  degrees 

Explosive:   Composition  B 
p   =16  80  kg/m 

Miss  distance:   r   =  50  ft  =  15.24  m 

m 

Figure  (4-13)  is  a  plot  of  the  output.   From  this 

figure  it  can  be  seen  that  a  warhead  radius  of  r   =  .06  m  is 

w 

required  to  achieve  a  conditional  kill  probability,  P   =  1.0. 
Therefore  the  missile  diameter  required  for  warhead  considera- 
tions is,  d  =  4.72  inches. 

The  warhead  radius  also  varies  with  the  required 
initial  velocity.   From  the  equations  for  charge  to  mass  ratio 
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and  casing  thickness  the  radius  required  to  achieve  a  specified 

kill  probability  can  be  determined. 

Since, 


t  = 


w 


C  Pc 

M  p 
Ke 


and 


C 

M 


V.  /2E 

l  ' 


1  -  V.  /2(2E) 
1  ' 


the  average  number  of  hits  becomes, 


ai = 


_    3 

5p   p  .   W  r 
Kc   K/H     w 


V.  /2E 
l 


,1  -  V.  /2(2E)/ 


-£+  1 

pe 


mr 


m 


The  above  equation  assumes  r   <  r 
^  m    c 


Letting 


b  = 


5  c  PK/H  W 

• 

2 

/V.2/2E        \ 
\1  -  Vi2/2(2E)/ 

pc 
pe 

+  1 

mr 
m 

the  probability  of  kill  becomes, 

2 


-br 


PD  =  1-e 


w 


(15) 


If  the  conditional  probability  of  kill  is  selected  as  P   =  0.999, 
equation  (15)  can  be  solved  for  the  warhead  radius. 
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'w=   I^M1/3 


The  initial  velocity,  V. ,  required  to  achieve  a  target  kill  is 
a  function  of  the  miss  distance,  fragment  size,  expected  en- 
counter altitude  and  target  characteristics.   For  this  simple 
analysis  the  effect  of  initial  velocity  on  warhead  radius  will 
be  studied.   The  effect  of  varying  the  initial  velocity  in 
equation  (16)  is  plotted  in  Figure  (4-14) .   The  results  were 
determined  for  various  explosives  to  show  their  effect  on  the 
warhead  radius.   If  an  initial  velocity  of  7000  ft/sec  is  chosen 
with  Composition  B  as  the  explosive  the  required  warhead  radius 
is  2.6  inches. 

From  this  analysis  the  warhead  weight  can  be  found.   The 
casing  thickness  is  given  by, 


t  =  .115  r    (from  required  C/M) 
w 

i      =  5  r 
w      w 


Therefore,  the  explosive  weight  is 


W  =  ir(r   -  t)2Z      p 
e      w       we 


W  =  5.9530  kg 
e 


The  casing  weight  becomes, 


2  2 

W      =    (irr        -    nr      )  I      p 
c  w  ewe 


W      =    6.8651   kg 
c  r 
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The  total  warhead  weight  becomes, 


WTT„  =  12.8181  kg  =  28.26  lb 
WH  3 


The  required  diameter  for  the  radar  antenna  was  10  inches 
and  exceeds  that  required  for  the  warhead.   If  the  warhead  is 
made  hollow  and  kept  at  the  same  weight,  equation  (14)  can  be 
reformulated  to  give 


C 

M 


.  K-^-i2] 


2r  t  -  t' 
w 


The  hollow  portion  of  the  warhead  is  then  found  from 


We     /   2 

—  =  it   r   -  r. 
p       \  w      1 


fcw 


Figure  (4-15)  is  the  resulting  warhead 


r .  =  2 . 76  in 

l 


t  =  .12  in 


Figure  (4-15) .   Hollow  warhead, 
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V.   BASELINE  DEFINITION 

A.   CONTROL  CONCEPTS 

The  lifting  and  control  surfaces  of  a  missile  may  be  of 
monowing,  triwing  or  cruciform  configuration.   Figure  (5-1)  shows 
these  three  arrangements. 


Monowing 


Triwing 


Cruciform 


Figure  (5-1).   Control  configurations. 

The  monowing  arrangement  is  typical  of  most  cruise  missiles, 
which  require  long  range  and  low  drag.   For  this  type  of  arrange- 
ment the  missile  must  bank  to  orient  the  lift  vector  for  a  maneuver 
Because  of  this,  the  monowing  missile  is  not  as  rapid  in  maneu- 
vering as  the  cruciform  configuration,  which  can  produce  lift 
in  any  direction  instantaneously.   The  cruciform   control  also 
has  identical  pitch  and  yaw  characteristics  which  results  in 
a  simpler  control  system.   The  triwing  configuration  is  used 
very  seldom  for  conventional  missiles.   It  can  be  shown  that 
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the  triwing  missile  requires  larger  wing  size;  therefore,  there 
is  little  drag  savings  even  though  there  is  one  less  wing 
than  with  the  cruciform  missile  [10].   From  this  discussion  it 
can  be  seen  why  the  cruciform  configuration  is  the  most  com- 
monly used  for  tactical  missiles  of  short  or  medium  range. 

The  positions  of  the  lifting  surfaces  on  the  missile  body 
depends  on  the  method  of  control  used  for  the  missile.   There 
are  three  conventional  methods  of  control  for  tactical  missiles. 
These  are;   1)  Wing  control;   2)  Tail  control;   3)  Canard 
control. 

1.   Wing  Control 

Wing  control  missiles  normally  have  large  movable  wings 
located  slightly  behind  the  missile  center  of  gravity.   Because 
of  the  small  moment  arm  the  wing  surface  must  be  relatively 
large  to  provide  control  effectiveness.   As  would  be  expected, 
larger  activators  are  required  for  moving  the  wings.   A  positive 
deflection  of  the  wings  causes  a  positive  normal  force;  there- 
fore, the  missile  reacts  almost  instantaneously;  thus  making 
wing  control  the  fastest  reacting  method  of  control.   Because  of 
the  smaller  moment  arm,  the  resulting  smaller  pitching  moment 
and  instantaneous  lift  result  in  smaller  angles  of  attack.   This 
feature  makes  wing  control  missiles  attractive  for  applications 
where  theincidence  angle  must  be  kept  small.   Air  breathing 
applications  often  use  wing  control  because  of  inlet  perfor- 
mance degradation  at  higher  angles  of  attack.   This  type  of 
control  is  also  good  for  fixed  seekers. 
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2 .  Canard  Control 

Canard  control  is  normally  the  method  of  control  where 
the  movable  surfaces  are  placed  as  far  forward  as  possible. 
Because  of  the  resulting  larger  moment  arm,  smaller  surfaces  are 
required  to  provide  control  effectiveness.   Lift  on  the  missile 
is  still  developed  primarily  by  the  aft  (wing)  lifting  surface. 
Response  is  slower  than  wing  control  because  of  the  need  to 
pitch  the  missile  to  an  angle  of  attack  before  lift  is  developed. 
Higher  angles  of  attack  are  needed  to  generate  the  required  lift. 
One  advantage  of  this  type  of  control  is  convenience  of  packaging. 
Since  the  controls  and  avionics  are  forward  of  the  propulsion 
system,  the  need  for  connectors  is  elminated.   For  stability 
reasons  the  wing  of  a  canard  control  missile  must  be  located 
farther  aft  than  a  conventional  wing-tail  configuration.   The 
zero  lift  drag  is  normally  lower  than  wing  control  missiles 
due  to  smaller  surfaces. 

3.  Tail  Control 

Like  the  case  of  canard  control,  a  tail  control  missile 
has  the  movable  surface  as  far  from  the  center  of  gravity  as 
possible.   This  also  results  in  a  larger  moment  arm  and  there- 
fore smaller  surface  required.   Control  deflection  is  the  oppo- 
site of  that  for  canard  or  wing  control  since  a  negative  control 
deflection  results  in  a  pitching  moment  that  pitches  the  nose 
up  and  therefore  a  positive  lift  on  the  main  lifting  surface  (wing) 
Tail  control  is  normally  the  slowest  method  of  control.   One 
advantage  is  that  the  flight  controls  are  at  the  end  of  the 


136 


missile  requiring  the  propulsion  system  to  be  located  forward 
of  other  types  of  control.   This  results  in  less  center  of  gravity 
movement  as  the  propellant  grain  burns.   This  becomes  a  definite 
advantage  for  longer  range  missiles. 

B.   GROSS  WEIGHT  AND  CENTER  OF  GRAVITY 

Since  the  component  weights  and  their  precise  locations  can- 
not be  determined  at  this  point  in  the  design  process,  some 
method  of  estimating  the  total  weight  and  center  of  gravity  is 
needed.   There  are  two  approaches  commonly  used  to  find  the  gross 
weight  of  missiles  in  the  conceptual  design  phase.   One  is  through 
the  use  of  the  historical  data  discussed  in  Chapter  2.   Since 
the  warhead  weight  is  known,  an  estimate  of  the  gross  weight 
is  now, 


W 

:[/      WH 

7AVG 


Another  method,  which  is  used  extensively  in  the  design  of 
aircraft,  is  the  use  of  regression  formulas  to  find  gross  weight 
or  component  weights  in  terms  of  parameters  that  are  known  early 
in  the  design.   Reference  (7)  has  derived  such  a  formula  for  the 
gross  weight  of  a  missile.   It  is, 

WQ  =  KG(L)2*13  (D)1'14 

K   =  Constant  to  be  determined  ,,, 

L  =  Total  missile  length  (inches) 
D  =  Missile  diameter  (inches) 
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The  constant  K   in  equation  (1)  is  derived  from  a  baseline 
G 

(generic)  missile.   In  this  case 

(W  J 

_  baseline 

K  _  — 


^        2  13  1  14 

(L)^*1J       (D)1*14 

baseline    baseline 

The  accuracy  of  a  regression  formula  such  as  equation  (1)  depends 

upon  how  close  the  synthesized  missile  is  to  the  baseline  for 

which  K  was  determined.   If  the  parameters  used  vary  more  than 
G 

20  percent,  the  accuracy  of  the  regression  equation  decreases 
rapidly;  therefore,  if  the  parameters  length  and  diameter  vary 
significantly  from  those  of  the  baseline,  equation  (1)  may  give 
an  erroneous  estimate  of  the  gross  weight.   As  an  aid  in  deter- 
mining the  gross  weight,  the  following  values  of  K  were  derived 

G 

from  data  given  in  reference  [7]. 

Missile        Figure        K 

G 


SRAAM  (5-2)  .00128 

MRAAM  (5-3)  .00118 

LRAAM  (5-4)  .00093 

SAM  (5-5)  .00108 

Figures  (5-2) ,  (5-3) ,  (5-4) ,  and  (5-5)  show  the  generic  missiles 
from  which  these  values  were  derived. 

The  center  of  gravity  of  the  baseline  missile  for  this  section 
is  taken  to  be  60  percent  of  the  total  length.   At  this  point 
sufficient  information  has  been  developed  to  define  the  baseline 
missile  from  which  design  iterations  will  be  made.   The  lifting 
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surface  planform  for  the  baseline  is  taken  as  a  delta  planform. 
As  will  be  shown  later,  this  wing  planform  will  be  very  close  to 
an  optimum  wing. 

C.   DESIGN  EXAMPLE  (BASELINE  DEFINITION) 

The  threat  for  which  this  missile  is  designed  is  a  highly 
maneuverable  fighter.   For  this  threat  a  canard  control,  cruci- 
form configuration  is  chosen.   The  rationale  is  that  this  configur- 
ation will  provide  the  fast  response  necessary  at  minimum  drag. 
The  diameter  was  fixed  at  10  inches  due  to  antenna  considerations 
in  Chapter  3.   From  historical  data  in  Chapter  2, 

L  m    (k)    D .  -    (15.89)     10    =    158.9    inches 
L 
LN    =     (    ET     °    =    22'3    inches 


The  MRAAM  of  Figure  (5-3)  is  used  as  a  generic  missile  for 
selecting  K_, .   Inserting  length  and  diameter  into  equation  (1), 
the  first  estimate  of  gross  weight  becomes, 

W„  =  .00118  (158. 9)2'13  (10)1*14  =  794.83  lb 


The  total  lifting  surface  required  is  then,  from  historical 
average, 

s  "  (I'AVG  (V  "  'sFTog''794-33' 

S  =  9.02  ft2 
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Since  the  canard  to  wing  area  ratio  is  known  from  historical 
data, 


S  /ST7  =0.20 
c/  W 


ST7  =  7.52  ft' 
W 


S   =  S  -  Sr  =  1.50  ft' 
C         W 


Lifting  surfaces  (delta  planform) .   From  historical  data, 


AR.,  =  1.61 


lRw 


AR   =3.74 
c 


b/2 


AR  =  — 

s 


Figure  (5-6) .   Lifting  surface 


Wi 


ng:   (b)W  "  ^ARW  SW   =  3'48  ft 


(C  )   =  2  S  /b  =  4.32  ft 
W 


Canard: 


(b)   =VAR   S    =  2.37  ft 
c      c  c 


(C  )   =   1.27  ft. 
r 
c 
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The  baseline  missile  is  now  defined.   The  canards  are  placed 
as  far  forward  as  possible.   The  wings  are  placed  as  far  aft 
as  possible  to  ensure  the  center  of  pressure  is  behind  the 
center  of  gravity.   The  exact  location  of  the  wings  will  be 
modified  in  the  next  chapter.   Figure  (5-7)  is  a  drawing  of 
the  baseline  missile. 
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VI.   LINEAR  AERODYNAMICS 

The  total  drag  on  the  missile  is  used  to  size  the  propul- 
sion needed  for  the  cruise  segment  of  flight.   The  wing  size 
depends  upon  the  maximum  lift  required  by  the  missile.   The 
wing  and  tail  are  normally  placed  to  provide  minimum  drag  during 
cruise  or  a  certain  stability  margin  at  launch.   To  make  any 
of  the  above  calculations,  values  of  the  aerodynamic  coefficients 
are  needed.   This  chapter  presents  the  background  aerodynamic 
theory  necessary  for  these  initial  calculations.   The  theory 
used  is  linear  aerodynamic  theory  and  slender  body  theory,  from 
which  simple  calculations  can  be  made.   Where  linear  theory  did 
not  apply,  an  attempt  was  made  to  find  existing  empirical  expre- 
sions,  which  yield  results  accurate  enough  for  initial  calcula- 
tions.  The  full  nonlinear  theory  will  be  presented  in  Chapter  9. 
The  reference  area  for  all  coefficients  in  this  report  is  the 
missile  maximum  cross  sectional  area.   The  reference  length  is 
the  maximum  missile  diameter. 

A.   MISSILE  DRAG 

The  total  drag  of  a  missile  consists  of  zero  lift  drag,  C   , 

D0 
and  induced  drag,  C 

i 

C^  =  C^   +  C^  (1) 

D     Dn  D. 

0      l 

The  zero  lift  drag  can  be  found  from  a  component  build  up  method 
in  which  the  contributions  due  to  the  nose,  body  and  lifting 
surfaces  are  added  together  to  obtain  the  total  zero  lift  drag. 
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Care  must  be  taken  to  reference  the  appropriate  areas  when 
using  this  method.   The  total  is  then  multiplied  by  1.25  to 
account  for  interference  effects  and  variations  in  skin  roughness 
[12]. 


C„   =1.25 
D0 


(V  +  (v  +  (v  +  (cd  > 

0  N       °  B       °  W       °  T 


(2) 


The  method  used  to  find  the  components  drag  depends  upon  the 

speed  regime  in  which  the  missile  is  operating.   Since  this 

report  is  concerned  primarily  with  supersonic  tactical  missiles, 

supersonic  zero  lift  drag  will  be  discussed  here.   The  component 

supersonic  zero  lift  drag  can  be  divided  into  skin  friction, 

C   ,  and  wave  drag,  C   . 
f  W 


(C   )   =  (C   ,   +  (C   ) 

0  w     f  w     w  w 


1.   Supersonic  Skin  Friction 

The  flow  over  a  body  traveling  at  supersonic  speeds  is 
likely  to  be  turbulent;  so  the  incompressible  skin  friction  co- 
efficient is  given  by, 


C         =      lill (3 

f-  m  d    >2.58  ^ 

i  dog10    Re) 


The  Reynolds  number  in  equation  (3)  is  based  upon  the  cruise 
altitude  and  speed  and  upon  the  characteristic  length  for  the 
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component  being  determined.   The  Reynolds  number  is  given  by, 

P  VM  * 
Re  ~  M 

where    x  =  L.,,  the  length  of  the  nose 

L   ,  the  length  of  the  body  without  the  nose 

(c) ,  the  mean  aerodynamic  chord 

The  compressibility  correction  to  equation  (3)  is, 

-  -0.58 
Cf  -  C-   (1  +  0.15  MMZ)  (4) 

i 

From  equation  (4)  the  skin  friction  drag  coefficient  for  each 
component  can  be  found  when  referenced  to  the  appropriate  area. 

(c_  i    =  (c.)    iSwfted)N 

Df  N     f  H    Sref 

(s      ) 
1  wetted' 

(CD  »   =  (Cf>   S " 

Df  B       B    Sref 

The  lifting  surface  skin  friction  drag  is  determined  in  a  similar 
manner.   Care  must  be  taken  to  include  all  surfaces  in  the 
wetted  area  calculation  for  the  lifting  surfaces. 
2.   Supersonic  Wave  Drag 

The  supersonic  wave  drag  consists  of  components  contri- 
buted by  the  nose  and  lifting  surface.  Nose  wave  drag  depends 
on  the  shape  of  the  nose,  and  the  most  common  nose  shapes  are 
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conical,  ogival  and  hemispherical.   Reference  [10]  lists 
empirical  formulas  for  finding  the  form (wave)  drag  of  various 
nose  shapes  at  zero  angle  of  attack. 

(1)  Conical;   c   =  (0.083  +  0 .  097/M  2)  (a/10)  la  69 

W 

a  =  tan  -^= —  =  Nose  semi-vertex  angle  (Degs) 
2LN 

The  center  of  pressure  for  a  conical  nose  is  at  the  centroid  of 
the  nose  planform  or  two  thirds  the  length  of  the  nose. 

(      2[196(L  /D)2  -  16] 

(2)  Ogival;   C    =  P<  1 = j 

UW     (       28 (M  +  18) (LN/D)Z 


P  =  (C   )  for  conical  nose 
W 


The  center  of  pressure  for  an  ogive  noise  is, 


N 


50  (M  +  18)  +  7M  P  (5M  -  18) 
40  (M  +  18)  +  7M2P  (4M  -  3) 


The  semi-vertix  angle  for  an  ogive  is  twice  the  equivalent  cone 
angle. 

(3)  Hemispherical;  The  drag  on  a  hemispherical  nose  is  extremely 
high  compared  to  other  nose  shapes,  and  is  difficult  to 
estimate.   An  initial  estimate  of  the  wave  drag  can  be 
found  from  Figure  (6-1) . 
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The  wave  drag  due  to  the  lifting  surfaces  can  be  found 
using  the  methods  of  reference  [12] .   For  a  double  wedge  air- 
foil with  sharp  leading  edges  as  shown  in  Figure  (6-2) ,  the 
wave  drag  is  given  by  the  following  formulas: 

(1)   Supersonic  leading  edge, 


C    =i(t2Jw 

Dw     *    c      s 


ref 


(2)   Subsonic  leading  edge, 


2  S 

cn   =  B  cot  A    (i)  — 5L- 

DTT  LE   c  S   - 

W  ref 


-*-*t 


Figure  (6-2) .   Double  wedge  wing  [12] . 


where 


B  = 


c/xt 
1  -  xt/c 


A    =  leading  edge  sweep 


8  " 


rM   -  1 


STT  =  planform  area 
W 
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3.   Base  Pressure  Drag 

The  drag  contribution  of  the  blunt  base  for  non-boattailed 
bodies  is  given  in  reference  [4]  as, 


c '   =  -C 
UAB     ^P 


B 


2 


2  \1.14/  x  \  2.8 


Y+l 


*M 


2yM  -  (  -1) 

M 

Y+l 


-1 


This  term  assumes  no  jet  thrust  from  the  base  of  the  missile, 
or  that  the  missile  is  operating  in  the  power  off  condition. 
This  term  is  not  included  in  equation  (2) ,  which  is  not  a  bad 
assumption  for  powered  flight  where  the  nozzle  exit  area  is 
approximately  equal  to  the  base  area  of  the  missile.   For  the 
case  where  the  nozzle  exit  area  is  much  less  than  the  base 
area  as  in  Figure  (3),  the  base  pressure  contribution  should  be 
included. 


Figure  (6-3).   Base  pressure  areas  [10]. 
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In  this  case  the  base  drag  is, 


S, 

c_  =  c   '  - 


AB     AB    S   , 
ref 


S,  =  shaded  area  of  Figure  (6-3) 


4 .   Induced  Drag 

The  induced  drag  on  a  missile  is  the  drag  due  to  lift. 
This  drag  is  caused  by  the  component  of  the  lift  vector  in  the 
drag  direction.   For  supersonic  flow  the  induced  drag  is  given 

by, 


c   =  -i-  c  2 

l     Not 


where; 


c.  =     w 


L    q  S   - 
^   ref 


The  lift  curve  slope,  C   ,  will  be  developed  later  in  this  chapter 

B.   DESIGN  EXAMPLE  (ZERO-LIFT  DRAG  CALCULATION) 

The  thickness  to  chord  ratio  of  the  wing  and  tail  have  not 
yet  been  determined;  however,  it  is  desirable  to  construct  the 
lifting  surface  as  thin  as  structurally  possible  to  minimize  the 
wave  drag.   Since  structures  have  not  been  covered  the  minimum 
thickness  to  chord  ratio  is  estimated  at  3  percent.   The  flight 
and geometric  conditions  determined  thus  far  are; 
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h   =    10,000    ft 

*M    =    3'° 

S       _   =    0.5454 

ref 

(t/c^    (t/c)c-0.03 


L._   =    1.8  58  3 

N 

LCB  =  11.3833 
(c)w  =  2.8816 
(c)c   =    0.8467 


ST7   =    7.52    ft 
W 

Sc   =    1.50    ft 


From  these  conditions  the  following  table  can  be  constructed. 


Nose 

Afterbody 

Wings 

Canards 

X 

1.8583 

11.3833 

2.8816 

0.8467 

wetted 

3.3394 

29.8014 

15.0400 

3.0000 

R 
e 

2.984 

18.282 

4.6280 

1.360 

cfi 

0.0025 

0.0020 

0.0024 

0.0029 

cf 

0..0015 

0.0012 

0.0015 

0.0017 

\ 

0.0095 

0.0652 

0.0827 

0.0187 

\ 

0.1545 

____ 

0.0400 

0.0050 

c^ 

0.1640 

0.0652 

0.1227 

0.0237 

'D 


The  zero  lift  drag  for  the  wings  and  canards  in  the  above  table 
take  into  account  that  there  are  two  sets  of  wings  (4  panels) . 
The  total  missile  zero  lift  drag  then  becomes, 


(c  )  +  (c  )   +  (c  )  +  (c  ) 

0  N        °  AB        °  W        °  C 


C^  =  0.4695 
D0 
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C.   MISSILE  LIFT  CURVE  SLOPE 


Figure  (6-4).   Wing-body-tail  lift  [13]. 

For  the  purposes  of  this  section  lift  on  a  wing-body-tail 
combination  can  be  taken  as  the  sum  of  the  components  in  Figure 
(6-4).   These  consist  of, 


w 


B 


B 


L„  =  lift  on  the  nose 

=  lift  on  the  wing  in  the  presence  of  the  body 

=  additional  lift  on  the  body  due  to  the  presence 
of  the  wing 

=  lift  on  the  tail  in  the  presence  of  the  body 

=  additional  lift  on  the  body  due  to  the  presence 
(T)    of  the  tail 


N 


(B) 


(W) 


(B) 


The  lift  of  only  the  wing-body  combination  can  be  defined  as,  L  , 


where, 


Lc  =  KcLw 


(5) 
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The  lift  of  the  wing  alone,  L  ,  is  that  obtained  from  thin 
airfoil  theory  or  experiment  considering  the  exposed  area  only 
The  constant  K   in  equation  (5)  is  defined  as 

K   =  K     +  K 
C     B(W)     W(B) 


From  this  equation,  K     is  the  ratio  of  the  additional 

B(W) 
body  lift  in  the  presence  of  the  wing  to  the  wing  alone  lift  for 

zero  control  deflection,  5=0. 

t  (ct  > 

Lr  La  B 

K      =    (W)  m  °(W)        5  m    Q 

B(W)     LW         ^Lo1,, 

w 

K      -     (B)  -         (B)         fisB  0 
W(B>      Lw        (CLa'w 


The  interference  factors,  K     and  Y-  have  been  deter- 

B(W)       ^(B) 
mined  from  slender  body  theory  for  wing-body  combinations. 

Figure  (6-5)  is  a  plot  of  these  values.   The  wing  alone  lift 

curve  slope,  (C_  )   is  determined  from  thin  airfoil  theory  or 

La  W 
experiment.   The  lift  curve  slope  of  delta  wings  with  super- 
sonic leading  edges  is  given  by, 


CLa  =  J       where    8  -aK2-1     k  K   £ 


For  subsonic  leading  edges  this  becomes. 


c   m    2TTtane/tany  k  >  £ 

La       E6 
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Figure  (6-5).   Interference  factors  [12] 
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Where  E  is  the  elliptic  integral  of  the  second  kind  for 


v  1  -  (tane/tany)  , 


Figure  (6)  shows  the  Mach  angle  y  and  sweep  angle,  z 


Supersonic  LE 


Subsonic  LE 


Figure  (6-6).   Wing  leading  edge. 


The  lift  of  the  nose  is  that  obtained  from  slender  body 
theory.   For  small  angles  of  attack 

«W    =  2/rad- 

1   N 


The  tail-body  combination  lift  is  determined  in  the  same 
manner  as  the  wing.   The  above  equations  are  defined  for  the  case 
where  the  incidence  angle  is  zero  and  the  angle  of  attack,  (a), 
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is  varied.   Analogous  terms  can  be  defined  for  the  case  where 

the  angle  of  attack  is  zero  and  the  incidence  angle  (6)  is  varied, 

LR  L<5     R 

v  (T)       _  B(T) 

k  =   -j =   — 7^ < ;  a   =    0 

B(T)  LT  {CLa'T 

L  (C.J 

k  =  (B)     =  (B)        .  „    =    „ 

T(B,  LT  <CLa>T  ' 


The  interference  factors  k     and  k     were  also  found 

B(T)       i(B) 
from  slender  body  theory  and  are  plotted  in  Figure  (6.5) . 

From  the  above  defnitions  the  total  missile  lift  curve  slope 

can  be  found.   If  the  analysis  is  for  small  angles  of  attack 

so  that  L~N. 


(N)    =  (N)„  +  (N)__     +  (N)n     +  (N)_       ,  . 

N       W(B)       B(W)       T(B)  +  (N)B(T) 

in  coefficient  form, 


CN  «  Sreg  "  (CN>   «  Sref  +  (CN>     <*  Sw  +  (cn'     q  SW 

N  W(B)  B(W) 


+  (C„)      qm  Sm  +  (C„)      qm  Sm 
N  _     ^T   T      N '        T   T 

X(B)  B(T) 

If  the  downwash  is  neglected  so  that  q   =  q,  and  the  above 
equation  is  differentiated  with  respect  to  angle  of  attack, 
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s  s 

(C„    )  =    (C„    )       +     (C„    )  ^L.      +     (c       ) 


NacM  Na    N  Na  ^  Sref  Na   B^     Sref 


S  S 

+    (CNa>  S^7+    (CN«L  — 


T(B)       ref  B(T)       ref 


From  the  definition  of  the  interference  factors. 


(C„J    =  (C„  J   +  (C_,   (1L,    +  K     ) 


Nff        Na'       Na)     w  B ,._/   S   , 

CM        N        W     (B)      (W)     ref 


+  (CNa}   (KT     +  KB    }  S"2-  (6) 

Na  T   T(B)     B(T)   Sref 


If  the  tail  is  the  control  surface,  a  similar  development  for 

the  control  effectiveness.  C„ r  yields ~ 

N  o  J 


(cM.)   =  (cXTJm  (km    +  k„   ) 


N6  CM      Na  T    T(B)     B(T)   Sref 


D.   MISSILE  PITCHING  MOMENT 

From  Figure  (6-7)  the  moment  about  the  center  of  gravity 
can  be  found.  The  centroid  of  the  wing  is  the  location  of  the 
center  of  pressure  of  a  wing  alone  in  supersonic  flow.   The 
effect  of  the  wing-body  combination  is  to  move  the  center  of 
pressure  aft  and  as  can  be  seen  in  Figure  (6-7)  the  center  of 
pressure  for  the  additional  lift  on  the  body  due  to  the  presence 
of  the  wing  is  aft  of  this  location.   Reference  [13]  has  an 
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N 


N 
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-W  *W(B) 


N 
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x 


T(B) 


'B(T) 


Figure  (6-7).   Forces  Acting  on  the  Missile. 

excellent  discussion  on  how  to  find  these  centers  of  pressure. 

For  the  purpose  of  this  chapter,  which  is  initial  sizing  and 

placement  of  the  lifting  surfaces,  the  center  of  pressure  for 

both  of  these  forces,  N     and  N     is  taken  as  the  centroid 

W(B)       B(W) 
of  the  wing  planform.   With  this  assumption,  the  moment  about 

the  center  of  gravity  is, 


(M)„M  =  N„  SC_.  +  (N      +  N 

CM     N   N      W ,  _  s 


B(W)'  ** 


)  x,,  +  (N 


+   N        )   Xm 

r(B)      B(T)     T 


Following  the  same  development  as  for  C 

Ma  Na  ..      D  Na   TT        W,_.»  B ,._.  D      S       - 

N  W  (B)  (W)  ref 


+    «L,    )        (Km  +    K  )  (-i)  (tt-^-) 

Na '  T.„x  B/mN  D      S       - 

T  (B)  (T)  ref 


(7) 
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Also  if  the  tail  is  the  control  surface. 


Xm     S 

M6      Na  T    T(B)     B(T)    D   Sref 


Ome  must  be  careful  in  defining  the  moment  arms  in  the 
above  equations.   If  a  nose  up  pitching  moment  is  developed  in 
Figure  (6-7),  the  moment  arm  is  positive.   Conversely,  a  nega- 
tive moment  arm  means  a  nose  down  pitching  moment  is  developed. 

E.   WEIGHT  AND  CENTER  OF  GRAVITY  VARIATIONS 

From  the  preceding  sections  of  this  chapter,  it  can  be  seen 
that  the  analysis  is  normally  for  one  point  in  the  flight  pro- 
file.  Since  the  launch  condition  has  been  defined,  the  tail 
sizing  can  be  accomplished  for  this  initial  condition.   The 
variation  of  missile  weight  and  center  of  gravity  is  due  to 
propellant  burning.   For  a  solid  propellant  this  variation  can 
be  quite  large.   Since  no  information  is  available  on  the  pro- 
pellant at  this  point  in  the  design,  the  following  guidelines 
will  serve  for  initial  calculations.   For  air-to-air  missiles 
the  propellant  is  approximately  35  percent  of  the  launch  weight. 
The  center  of  gravity  travel  is  5  percent  of  the  length.   For 
surface-to-air  missiles  which  must  be  boosted  to  flight  speed, 
the  propellant  weight  can  be  taken  to  be  48  percent  of  the 
launch  weight  with  a  corresponding  center  of  gravity  travel  of 
approximately  8  percent  of  the  length.   These  are  initial  approxi- 
mations taken  from  historical  data  and  can  be  refined  later  in 
the  design  process. 
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VII.   LIFTING  SURFACE  DESIGN 

From  the  analysis  of  wing  lift  and  drag  in  Chapter  6,  it  can 
be  seen  that  the  performance  of  a  wing  will  vary  with  planform. 
Up  to  this  point  the  lifting  surfaces  have  been  considered  delta 
planforms  with  zero  taper  ratio.   This  is  not  necessarily  the 
optimum  planform  since  it  was  taken  from  a  historical  average 
and  does  not  apply  to  a  specific  missile.  This  chapter  deals  with 
the  sizing,  placement  and  planform  definition  of  the  missile. 

A.   WING  PLACEMENT 

The  wing  placement  on  the  missile  depends  upon  the  type  of 
control  used.   For  canard  control  the  wing  (aft  surface)  is  nor- 
mally fixed  as  far  aft  as  possible  for  stability  purposes.   For 
a  wing  control  or  tail  control  missile,  where  the  wing  is  near 
the  center  of  gravity,  the  wing  placement  becomes  more  critical 
and  depends  upon  the  stability  margin  required.   Since  the  drag 
during  the  cruise  segment  includes  the  drag  due  to  lift,  one 
method  to  minimize  drag  would  be  to  place  the  wing  such  that 
zero  lift  is  produced  on  the  tail  during  cruise.   The  moment  about 
the  center  of  gravity  is  zero  for  trimmed  flight,  thus  equation 

(7)  yields,  for  (C  )   =  0. 

T 

„    =     [ftfaja  *n  Sref 

W    "        (CM    )        ST7     (K  +    K  ) 

Na    W      W         W(B)  B(W) 
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B.   MANEUVER  LOAD  FACTOR 

Regardless  of  the  type  of  missile  being  designed,  it  will  be 
required  to  maneuver  in  order  to  intercept  its  intended  target. 
For  an  air-to-air  or  surface-to-air  missile  the  level  of  this 
maneuver  may  be  quite  large.   The  maneuver  load  factor  required 
was  found  in  Chapter  3.   The  maximum  maneuver  the  missile  can 
sustain  depends  upon  the  maximum  trimmed  normal  force  the  missile 
can  develop. 


Figure  (7-1) .   Sustained  maneuver  of  a  missile. 

From  Figure  (7-1)  the  force  developed  by  a  missile  in  a  con- 
stant acceleration  turn  is 


L  =  nW 

n  =  maneuver  load  factor 


(1 
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Small  angles  of  attack  are  assumed  for  which  the  lift  and 
normal  force  are  approximately  equal.   This  approximation  is 
good  for  angles  of  attack  of  up  to  10  degrees,  and  above  this 
value  of  maximum  trimmed  angle  of  attack  the  linear  theory 
becomes  inaccurate.   With  this  assumption  for  now,  equation  (1) 
becomes , 

L  =  N  =  nW 

(CNa)0  q  Sref  "  nW  <2> 

(C„J  nW 


Na         ,    (a)     q  S   ,. 
required    *  'max  *   ref 


Equation  (2)  gives  the  lift-curve  slope  required  to  develop 

the  required  normal  force  at  a  trimmed  angle  of  attack,  (a) 

max 

The  lift  curve  slope  developed  is  given  by  equation  (6)  in  Chapter 
6. 

(CNa>    =  ((W   +  «W  ,(KW    +  KB    >  S  ' 


CM        N        W   "(B)     (W)    ref 

(3) 

S 

+  (CW  (KT    +  KB    }  S"^ 
Na  T   T(B)     B(T)   Sref 


As  discussed  earlier  the  values  of  the  interference  factors 

depend  upon  the  wing  planform.   Typical  values  of  (IC.    +  K     ) 

W(B)     B(W) 

and  (K     +  K     )  are  from  1 . 5  to  2 . 0 .   As  a  conservative 
X(B)     B(T) 

estimate  a  value  of  1.5  is  assumed.   Equation  (3)  then  becomes, 

S  S 

(C„  )    =  (C„  )   +  1.5  (C„  )  ^—   +  1.5  (C„  )   =-i-     (4) 
Na  _,.,     Na  XT         Na  rT  S   -         Na  m  S   , 
CM        N  W   ref  T   ref 


166 


Selecting  supersonic  leading  edges, 


«WW=  <CNa»T  =  4^ 


Therefore 


(CNa>    "  (CNa>   +  TT—    <SW  +  V 

CM        N      ref 


If  the  complete  missile  lift-curve  slope  is  set  equal  to  the 
required  lift-curve  slope,  the  lifting  surface  area  required  to 
maintain  the  maneuver  load  factor  can  be  found,  using  equation  (4) 

<sw+  V      -  [(cNa>      "  <cnc>  ]  B-^irl 

req         req        N 


Let 


(SW  +  V 
K  =  ^ 


<sw  +  V 

baseline 


If  the  same  ratio  of  wing  to  control  surface  area  is  kept  to 
minimize  stability  perturbations,  the  new  wing  and  tail  area  are, 


SW=K(SW>     .. 

baseline 

ST  =  K(S  ) 

baseline 

The  above  analysis  assumes  a  linear  variation  of  C„  which  is 
good  for  small  angles  of  attack.   A  more  precise  analysis  will 
be  performed  in  Chapter  9.   The  above  analysis  should  be  performed 
at  the  expected  encounter  conditions. 
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1.   Design  Example  (Maneuver  Load  Factor) 

Since  a  canard  control  was  chosen  for  the  design  example 
the  canard  and  wing  position  are  fixed  as  fqr  forward  and  aft 
as  possible.   The  analysis  is  done  for  a  conservative  Mach  number 
of  2.5  after  the  missile  has  slowed  from  its  cruise  velocity  due 
to  the  maneuver.   From  previous  results, 

n  =  31  gs 

(S   +  S  )  =9.02  ft2 

baseline 

W_  =  794.83  lbs 

If  the  missile  is  required  to  maneuver  at  one-half  its  powered 
range,  and  the  propellant  weight  is  35  percent  of  the  launch 
weight,  the  maneuver  weight  becomes 

W  =  655.73  lbs 


The  dynamic  pressure  for  M„  =  2.5  and  at  10,000  feet  altitude  is 

M 


1       T7   2 

q  =  2  °  VM 

q  =  6369.84  lbs  ft2 

The  required  lift-curve  slope  becomes  for  (a    )  =  10  degrees 
=  .1745  rad. 


(C   )     =  33.5311  per  rad. 
req 

The  required  lifting  surface  area  is; 

6  S  , 

(ST7  +  S  )     =  [33.5311-2]  ~  =  6.57  ft 

req 
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Therefore    K  =  0.7284  . 

The  required  wing  and  canard  sizes  to  achieve  the  maneuver  are; 
ST7  =  (7.52)  (.7284)  =  5.48  ft2 

w 

Sc  =  (1.50)  (.7284)  =  1.09  ft2 

C.   TAIL  SIZING 

For  a  wing-tail  combination,  the  primary  concern  in  the 
sizing  of  the  tail  is  the  static  stability  of  the  missile.  The 
missile  becomes  more  stable  as  the  mission  proceeds,  since  the 
center  of  gravity  moves  forward  as  the  propellant  burns.   As 
the  missile  becomes  more  stable,  control  of  the  missile  becomes 
sensitive.   If  the  tail  is  sized  for  zero  static  stability  at 
launch,  the  missile  control  will  remain  more  effective  during 
flight.   This  is  the  best  condition  possible  without  the  use  of 
some  form  of  stability  augmentation  system  at  launch.   Therefore, 

at  launch  C.T  =  0 . 

Na 


0  "  (CNa>   XN  +  (CNa'   *W  S-^7  !KW,B,  +  KB,„  1 
N  W      ref     (B)      (W) 


(5) 


ST 

+  (CNa}   XT  S^^  (KT     +  KB    ) 
Na  T   T   ref    T(B)     B(T) 


With  the  lifting  surface  area  fixed  due  to  the  maneuver  load 
factor,  the  tail  can  be  sized  to  satisfy  equation  (5) .   As  can 
be  seen  the  position  depends  highly  on  the  moment  arms  that  the 
lift  forces  act  through.   Since  the  missile  length  has  not  been 
fixed  at  this  point  in  the  design  and  may  vary  due  to  propulsion 
requirements,  this  analysis  will  be  completed  later. 


D.   WING  PLANFORM 

The  wing  planform  is  specified  by  the  leading  edge  sweep, 
A   ,  taper  ratio,  A,  aspect  ratio,  AR  and  planform  area  S  .   The 
planform  area  was  fixed  due  to  maneuvering  requirements  in  a 
previous  section.   This  section  is  concerned  with  defining  the 
remaining  planform  parameters.   Figure  (7-2)  is  the  wing  plan- 
form  and  the  equations  used  to  define  these  parameters. 


i 


AR  = 

X  = 

c   = 

r 

c  = 
c  = 


b2/s 
ct/cr 

2  s 

b(l+X) 


1  cr  (  I+T 

2  tan  ALE 


i+x+y 


Figure  (7-2) .   Planform  geometry. 

The  value  of  the  lift-curve  slope  used  previously  was  derived 
from  linear  theory  and  is  applicable  only  to  simple  planforms. 
For  the  purpose  of  this  section,  which  will  include  low  aspect 
ratio  wings  with  supersonic  and  subsonic  leading  edges,  Figures 
(7-3)  and  (7-4)  are  used.   These  figures  have  been  corrected  for 
3-D  effects. 

1.   Effect  of  Taper  Ratio  and  Leading  Edge  Sweep 

The  lift  and  drag  characteristics  of  the  wing  are  the 
primary  parameters  of  interest.   The  objective  in  wing  design  is 
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Figure  (7-3).   Theoretical  wing  lift  curve  slope  [12]. 
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Figure  (7-4) .   Lift  curve  slope  continued  [12] 
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is  to  obtain  maximum  lift  with  minimum  drag.   It  will  be  shown 
that  these  objectives  are  usually  conflicting;  therefore,  some 
compromise,  or  optimum,  planform  must  be  found.   From  Chapter  6 
the  equations  for  the  drag  of  a  wing  in  supersonic  flow  are  found 
The  skin  friction  drag  depends  upon  the  mean  geometric  chord. 
Equation  (3)  in  Chapter  6  indicates  a  larger  mean  geometric 
chord  would  result  in  less  skin  friction  drag  for  a  fixed  area. 
The  mean  geometric  chord  can  be  expressed  in  terms  of  taper 
ratio  as  follows, 


3  b(l+A)  V   1+A 


5  =  I  \/ar 


1+X+X2 
(1+A)2 


This  equation  leads  to  a  zero  taper  ratio  to  maximize   c  and 
reduce  skin  friction  if  the  surface  area  and  aspect  ratio  were 
fixed.   The  wave  drag  is  constant  for  supersonic  leading  edges, 
and  decreases  when  the  leading  edge  goes  subsonic  or  the  leading 
edge  is  behind  the  mach  line.   The  lift  capabilities  also  de- 
crease as  the  leading  edge  goes  subsonic. 

As  stated  earlier  this  leads  to  conflicting  performance  since 
the  objective  is  to  minimize  drag  while  maximizing  lift.   At 
this  point  an  example  best  illustrates  the  results  of  varying 
taper  ratio  and  leading  edge  sweep.   The  lift  -  curve  slope  is 
derived  using  the  methods  of  reference  (1).   Figures  (7-3)  and 
(7-4)  are  from  reference  [12]  and  are  used  to  find  the  supersonic 
linear  lift  corrected  for  3-D  flow  effects.   The  drag  methods 
of  Chapter  5  are  used  to  determine  the  drag  characteristics. 


2.  Example 

2 
For  a  fixed  wing  area  of  4  ft   and  an  aspect  ratio  of 

2  the  lift  curve  slope  and  drag  were  determined  for  a  flight 
Mach  number  of  M=  2.   This  Mach  number  corresponds  to  a  sonic 
leading  edge  sweep  of  60  .   The  wing  leading  edge  sweep  was 
varied  from  zero  to  75   for  taper  ratios  of  0,  1/2  and  1.   The 
results  are  plotted  in  Figures  (7-5)  and  (7-6) . 

From  this  example  certain  generalizations  can  be  made. 
From  Figure  (7-5)  it  can  be  seen  that  the  drag  is  relatively 
insensitive  to  taper  ratio.   There  is  a  reduction  in  drag  for 
increased  leading  edge  sweep.   For  this  example  there  is  approxi- 
mately a  6  percent  drag  reduction  for  every  five  degree  increase 
beyond  the  sonic  value.   Figure  (7-6)  indicates  that  the  effect 
of  decreasing  the  taper  ratio  is  to  delay  the  drop  in  the  lift 
curve  slope  of  the  wing.   From  this  example  a  general  guideline 
would  be  to  fix  the  leading  edge  sweep  5  degrees  beyond  the  sonic 
value  and  the  taper  ratio  at  zero.   This  would  provide  a  6  per- 
cent reduction  of  wave  drag  while  maintaining  the  maximum  lifting 
capabilities  of  the  wing. 

3 .  Effect  of  Varying  Aspect  Ratio 

The  result  of  increasing  aspect  ratio  is  an  increase  in 
the  lift-curve  slope  of  the  wing  [12]  .   The  aspect  ratio  is 
given  by, 

AR  =  b2/s 

For  a  zero  taper  ratio  wing  this  becomes; 


AR  =  2b/cr  (6) 
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From  equation  (6)  it  can  be  seen  that  increasing  aspect  ratio 
results  in  a  corresponding  increase  in  wing  span  for  a  given 
root  chord.   Missiles  are  normally  span  limited,  due  to  launcher 
constraints;  therefore,  there  is  a  maximum  span  which  can  be 
accomodated.   The  drag  of  the  wing  again  is  in  conflict  with 
the  lift  since  increasing  aspect  ratio  decreases  the  mean  geo- 
metric chord  for  a  given  span  and  therefore  increases  skin  fric- 
tion drag.   A  compromise  wing  AR  must  be  found  which  considers 
both  the  lift  and  the  drag  characteristics.   For  a  discussion 
of  optimum  aspect  ratio,  the  following  functions  are  defined 
for  convenience: 

F  =  F   +  F2  s Lift-drag  function  (7) 

where  C 

— ^Normalized  drag  function 


1 "  (V 

0  max 

1/CNa 
F?  =  TTTr T  =  Normalized  lift  function 

A         X/  ltNaJ  . 
mm 


From  equation  (7)  if  F  is  plotted  over  the  allowable 

range  of  aspect  ratio,  a  minimum  value  of  the  lift-drag  function 

fixes  the  desired  aspect  ratio.   This  aspect  ratio  corresponds 

to  minimizing  the  drag  function,  F, ,  while  maximizing  the  lift 

function,  F_.   For  convenience  the  abscissa  is  plotted  as  c/c  ,    » 
2.  \  max ) 

as  shown  in  Figure  (7-7) . 

From  Figure  (7-7)  the  optimum  mean  geometric  chord  cor- 
responds to  point  A. 
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(max) 


Figure  (7-7) .   Lift-drag  function 


5=    s- 


max 


'max/  A 


The  aspect  ratio  is  then  given  by 


c 

r 

= 

3 
2 

c 

b 

= 

2s/c 

AR 

= 

b: 

) 

7s 

As  mentioned  earlier  a  missile  is  normally  span  limited.   The 
plot  of  the  lift-drag  function  F  is  normally  fairly  flat  on  the 
left,  or  for  c/c .         .  approaching  zero.   For  this  reason  point 
B  may  be  chosen  as  the  optimum  since  F  varies  very  little  up 
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to  this  region.   Point  B  corresponds  to  increased  chord  and 
decreased  span. 

4 .   Design  Example  (Wing  Planform) 

From  the  previous  analysis  the  following  parameters  were 
defined. 

M„=  3.0 

Sw  =  5.48  ft2 

If  the  leading  edge  sweep  is  fixed  5  degrees  behind  the  mach 
line. 

y  =  sin"1  ^  -  19.5° 

ALE  =  95  -  y  =  75.5° 

\   =  0 

The  planform  table  becomes 


AR 

b 
3.31 

c 

r 

c 
2.21 

1.17 

KB(W) 

2.0 

3.31 

.26 

1.5 

2.87 

3.82 

2.55 

1.19 

.30 

1.0 

2.34 

4.68 

3.12 

1.22 

.35 

.75 

2.03 

5.40 

3.60 

1.24 

.40 

.50 

1.66 

6.62 

4.41 

1.29 

.48 

.25 

1.17 

9.36 

6.24 

1.36 

.62 

The  wave  drag  is  constant  and  is  given  by, 
C..   =  4  cot  AT„  (^) 


DT7  LE   c    S  c 

W  ref 


C    =  .0094 
W 
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AR 

% 

\ 

AR  tan  A 

6 

tan  ALE 

2.0 

.0303 

.0397 

7.73 

.73 

1.5 

.0296 

.0390 

5.80 

.73 

1.0 

.0288 

.0382 

3.87 

.53 

.75 

.0282 

.0376 

2.90 

.73 

.50 

.0274 

.0368 

1.93 

.73 

.25 

.0260 

.0354 

.97 

.73 

(C._  )  (KT1         +K_    ) 
W    (B)     (W) 


AR 

c/c 

'  max 

Fl 

F2 

F 

2.0 

.35 

1.0000 

.3755 

1.3755 

1.5 

.41 

.9824 

.3999 

1.3823 

1.0 

.50 

.9622 

.4344 

1.3966 

.75 

.58 

.9471 

.4456 

1.3927 

.50 

.70 

.9270 

.5360 

1.4630 

.25 

1.00 

.8917 

1.0000 

1.8917 

2.1137 
1.9845 
1.8271 
1.7813 
1.4809 
.7937 


From  Figure  (7-8)  it  can  be  seen  that  the  lift-drag 

function  remains  relatively  constant  up  to  c/c    =  .58.   Since 

2  ^  '    max 

the  object  is  to  make  the  span  as  small  as  possible  this  point 
is  taken.   The  wing  planform  becomes/ 
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c  =  .58  c      =  3.62  ft  =  43.44  in 
max 

c   =  5.43  ft  =  65.16  in 
r 

b  =  2.02  ft  =  24.24  in 
c  =  3.90  ft  =  46.80  in 
AR  =  0.74 


5.   Design  Example  (For  Canard  Planform) 

From  previous  analysis  and  development,  the  following 
parameters  have  been  defined: 

M  =  3.0 

S   =  1.09  ft2 
c 

(t/c)   =  0.03 
c 

If  the  leading  edge  sweep  is  set  5  degrees  behind  the  sonic 
value  and  the  taper  ratio  is  set  to  zero,  the  following  planform 
table  results. 


*LE  =  ?5-5 


X   =  0 


AR 

b 

c 

r 

c 

*W(B) 

KB 

(W) 

2.5 

1.63 

1.34 

.89 

1.29 

.52 

2.0 

1.48 

1.47 

.98 

1.31 

.55 

1.5 

1.28 

1.70 

1.14 

1.35 

.59 

1.0 

1.04 

2.10 

1.40 

1.47 

.76 

.5 

.74 

2.94 

1.96 

1.51 

.90 

.25 

.52 

4.19 

2.79 

1.59 

1 

.02 
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AR 

% 

\ 

AR  tan 

*a 

8 

CNa(k) 

tan  ALE 
.  73 

2.5 

.0069 

.0087 

9.67 

2.0 

.0067 

.0085 

7.73 

.73 

2.7900 

1.5 

.0067 

.0085 

5.80 

.73 

3.6339 

1.0 

.0064 

.0082 

3.87 

.73 

2.5952 

.5 

.0061 

.0079 

1.93 

.73 

2.0880 

.25 

.0058 

.0076 

.97 

.73 

1.0124 

AR 

c/c 

max 

Fl 

F2 

F 

2.0 

.35 

1.0000 

.3628 

1.3628 

1.5 

.41 

1.000 

.3844 

1.3844 

1.0 

.50 

.9647 

.3901 

1.3557 

.5 

.70 

.9294 

.4849 

1.4143 

.25 

1.00 

.8941 

1.0000 

1.8941 

From  Figure  (7-9)  it  is  seen  that  there  is  a  minimum  of 

the  lift-drag  function  at  c/c    =  0.5.   The  optimum  canard  is 
r  max  r 

then, 

C  =  0.5  c     =  1.40  ft  =  16.80  in 
max 

c   =  2.09  ft  =  25.08  in 

r 

b  =  1.04  ft  =  12.48  in 
AR  =  0.99 
c  =  2.01  ft  =  24.12  in 
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E.   DESIGN  EXAMPLE  (REVISION  OF  ZERO  LIFT-DRAG  AND  LIFT-CURVE 

SLOPE) 

Wing:         c  =  3.62  ft 

c   =  .0014 

(C   )   =  .0750 
u0  W 

Canard:      c  =  1.4  0  ft 

C-  =  .0016 

(C   )   =  .0166 
UQ    c 

Body:    (C   )   =  .2292 
D0  B 

The  complete  missile  C_   including  the  interference  factor 

0 

of  1.25,  is  now  0.4010.   From  the  previous  calculation  of  the 
baseline  zero  lift-drag  coefficient,  the  drag  has  been  reduced 
by  14.6  percent. 

From  equation  (6)  in  Chapter  6  the  lift-curve  slope  is  now 


C..   =  25.17/rad 

Na 


Figure  (7-10)  is  the  missile  design  to  this  point 
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VIII.   PROPULSION  REQUIREMENTS 

The  following  discussion  presents  a  method  for  preliminary 
sizing  of  a  solid  rocket  motor  for  a  boost-sustain  trajectory 
of  an  air  launched  or  a  surface  launched  tactical  missile.   The 
analysis  consists  of  sizing  the  booster  from  incremental  velo- 
city considerations  and  sizing,  the  sustainer  for  the  maximum 
range  required  at  the  operational  altitude.   The  method  assumes 
a  constant  acceleration  boost  and  a  constant  altitude  cruise. 

A.   BOOSTER  INITIAL  SIZING 

Since  the  control  system  cannot  respond  properly  while  in 
the  boost  phase  of  flight,  it  becomes  important  to  make  the  boost 
time  as  short  as  possible.   The  limiting  factor  for  boost  time 
is  the  maximum  axial  acceleration  the  airframe  and  components 
can  withstand.   This  maximum  acceleration  is  normally  around 
30  gs.   The  boost  time  then  becomes, 

%  -   )       * 

vi 

For  constant  acceleration, 

V2"V1    AVB 
b     a      a 


The  incremental  velocity  during  boost,  AV_. ,  can  be  derived 

o 

from  the  equations  of  motion. 
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/ 


Figure  (8-1).   Forces  acting  on  the  missile. 

If  during  the  boost  phase  the  thrust,  drag  and  launch  angle  are 
considered  constant,  the  axial  acceleration  of  the  missile  is 
constant  and  may  be  written  as : 


dV 

m  jt=  T  -  D  -  W  sin  yT 
dt  'L 


T-D 
dV  =  g(— — )dt  -  g  sin  y      dt 
W  L 


If  the  launch  velocity  is  V,  and  the  velocity  at  the  end  of 


boost  is  V_ , 


V. 

J' 


dV  = 


g(-^-)dt  -  /    g  sin  yL 


(2 


The  vehicle  weight  in  equation  (2)  is  a  function  of  time.  If  the 
propellant  weight  is  given  by  W  ,  and  the  propellant  grain  burn 

IT 

is  linear  with  time,  the  missile  weight  can  be  expressed  as, 
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W  =  W_  -  W  (-^J 

L     p  tb 


where  W  =  Launch  weight. 
Equation  (2)  then  becomes, 


AVB  =  V2  -  V1  =  g(T  -  D) 


t  dt 
(WTt  -W  t)  -3  sin  YL  tb 


Lb  p 


Integrating  the  first  term  yields, 


AVB  = 


g(T  -  D)t} 
w 


£n 


(1  -  V  fcb 


-  g  sin  yL  tb 


Since  the  empty  weight  is  given  by 


W_  =  Wx  -  W 
E     L    p 


The  incremental  velocity  can  be  written  as 

g(T  -  D)t     W 

AVB  =  W ^  WT  "  g  Sln  YL  % 

p  E 


(3) 


A  first  estimate  of  propellant  weight  can  be  made  if  the  drag  is 
assumed  zero  in  equation  (3) ,  and  I    =  Tt. 


sp 


"b/w 


w. 


AV_  =  g  I    in   77-  -  g  sin  yT  t*. 
B    3   sp     W^  'Lb 

E 


(4) 


In  the  above  expression  the  specific  impulse,  I   ,  is  charac- 
teristic of  the  propellant  chosen  and  can  be  found  from  histori- 
cal data,  or  a  specific  propellant  value  may  be  used. 
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WL 
Equation  (4)  can  then  be  solved  for  ^—  and  the  first  esti- 

WE 

mate  of  propellant  weight,  W  ,  can  be  found  as  follows, 

WE 
w   =  W  (1  -  — ) 

p     Ll     W_; 

Li 


With  the  propellant  weight  and  specific  impulse  known,  the 
total  impulse  and  average  thrust  required  for  boost  can  be  found. 


I_  ■  I    W 

T     sp   p 

T  =  r- 


With  the  thrust  known  and  an  average  value  of  drag  assumed, 
equation  (3)  can  be  iterated  for  an  improved  estimate  of  the  pro- 
pellant weight  needed. 

The  final  result  is  used  to  calculate  the  booster  combustor 
volume  and  length. 


vb  "  wp/ppnp 

where  n  =  volumetric  packing  factor 

B      B/trD 
'      c 

The  combustor  diameter  is  limited  by  the  missile  maximum 
diameter  and  is  a  design  choice.   The  propellant  density,  p  , 
and  volumetric  packing  factor,  n  ,  are  determined  from  histori- 
cal data  or  given  for  a  specific  propellant. 
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B.   SUSTAINER  INITIAL  SIZING 

The  sustainer  thrust  required  to  maintain  cruise  is  the 

driving  factor  in  the  sustainer  motor  sizing.   For  initial  sizing 

purposes,  a  level,  constant  velocity  cruise  segment  is  assumed. 

In  this  analysis  the  thrust  required  is  equal  to  the  aerodynamic 

drag  developed  by  the  missile.   The  performance  requirements 

determine  the  maximum  range,  operating  altitude,  and  velocity  of 

the  missile.   From  these  requirements  the  sustainer  burn  time, 

t  ,  can  be  determined. 
s 

.   _   max B 

s  ~       VM 
M 

where  S   =  distance  traveled  during  boost. 

B 

Since  sustainer  thrust  is  equal  to  drag,  T   =  D,  the  total 
impulse  required  becomes, 


The  sustainer  nozzle  is  sized  for  the  operational  altitude 

of  the  missile.   The  thrust  coefficient,  C  , ,  can  be  expressed  as; 

F  d 


2 


Y+l  Y-l 


CFd   cdA  V  7-T  (Y+1} 


po 

pc 


The  constant,  Cd,  is  the  nozzle  efficiency  and  has  been 
determined  from  historical  data  to  be  0.96.   The  nozzle  half 
angle  correction,  A,  is  given  by; 
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A  =  (1  +  cos  a)/2 


where  a  is  the  nozzle  half  angle.   A  nominal  value  used  in  many 
designs  is  15  .   Larger  values  of  a  give  significant  non-axial 
flow  components  in  the  nozzle.   Values  less  than  15   give  large 
nozzle  lengths  and  therefore  contribute  to  excessive  missile 
weights.   The  combustion  chamber  pressure  is  also  a  design  choice 
A  lower  usable  value  of  P   is  200  psia  for  sustainers  and  500  psi 
for  boosters.   The  chamber  pressure  cannot  exceed  the  maximum 
expected  operating  pressure  of  the  missile  (MEOP) . 

The  thrust  of  the  sustainer  can  be  expressed  as  a  function 
of  the  thrust  coefficient,  chamber  pressure,  P  ,  and  throat  area, 
A  .   The  throat  area  required  to  deliver  the  required  sustainer 
thrust  can  be  determined  from 

T   =  C  ,  P   At 
s     Fd   c   t 

The  nozzle  area  ratio  is  a  function  of  the  pressure  ratio 

P./P  ,  where  P_  is  the  ambient  pressure  if  the  nozzle  expands 

the  flow  fully. 

Y+l 


*  -  -i  j  *  fi  +  (1+1)  m  21 

t   Me  |Y+1  L      2   Me  J 


|2(y-D 


where  P.     r        ,  .   -i    Y-l 


,  1  +  (-^)  M. 
c 
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From  the  previous  equation  the  exit  area  can  be  determined. 
From  the  throat  and  exit  diameters,  the  nozzle  length  is  found 
by 

D   -  D^ 
L  _  __e t 

n    2  tan  a 

At  this  point  in  the  design  the  nozzle  exit  area  and  length 
should  be  checked  to  see  if  they  exceed  any  specifications  on 
the  missile;  i.e.,  the  nozzle  exit  area  must  be  less  than  the 
base  area  of  the  missile. 

Finally,  the  delivered  specific  impulse  is  typical  of  the 
propellant  chosen  for  the  sustainer.   The  propell'ant  weight  is 
then, 

(Vs 

r     sp 

The  sustainer  combustor  volume  and  length  are  then  determined 
in  the  same  manner  as  the  booster.   The  equations  concerning 
the  sustainer  nozzle  also  apply  to  the  booster  nozzle  with  the 
appropriate  booster  areas,  pressures  and  thrust  substituted* 

C.   ROCKET  MOTOR  CHAMBER  PRESSURE 

As  indicated  earlier  the  delivered  thrust  increases  with 

increasing  chamber  pressure,  P  ;  however,  the  wall  thickness 

in  the  rocket  motor  depends  upon  the  expected  operating  pressure 

within  the  chamber.   If  the  wall  thickness  is  t  inches,  and  the 

yield  stress  of  the  casing  material  is  a  ,   it  can  be  shown  that 

Y 

the  thickness  required  of  the  motor  casing  is  given  by, 
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t  = 


P    r 
max 


(5) 


where  r  is  the  chamber  radius,  and  P    is  the  maximum  chamber 

max 

pressure,  which  is  taken  to  be  1 . 2  P   to  allow  for  variations 
within  the  propellant.   From  equation  (5)  the  weight  of  the  motor 
casing  can  be  determined  if  the  casing  is  cylindrical  and 
t  <  <  r 


W  =  p   V   =  2irri    tp 
c    c  c       c   c 


Substituting  for  t 


W  = 
c 


2.4-irr  I   p   P 

c  c   c 


(6) 


The  specific  impulse  of  the  rocket  motor  is  given  by 


V. 

I   =  -1 

sp  g„ 


-i   2g   (-M  R  T. 
gc  j  *c   y-1      0 


Y-l 
P0   Y 

l  -  (=2.) 

c 


1/2 


(7) 


From  equations  (6)  and  (7),  it  can  be  seen  that  a  compromise 
must  be  made  in  selecting  the  chamber  pressure;  since  increasing 
P   increases  the  specific  impulse  but  also  increases  the  casing 
weight.   An  attempt  must  be  made  to  find  an  optimum  operating 
pressure. 

If  the  rocket  motor  weight  consists  of  the  nozzle,  propellant 
and  casing, 


WM  =  WXT  +  w_  +  w_ 

M     N     P     C 


(8 
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and  the  total  impulse,  which  is  constant   is  given  by, 


I  =  Tt,  =1    W 
b    sp   p 


The  optimum  chamber  pressure  can  be  found  by  minimizing  the 
weight  and  maximizing  the  impulse. 


A  WM        A  WM 


dp   v  I'    dp   VI    W 
rc         rc        sp   p 


differentiating  yields. 


1     d   ,_M.  _  _M      sp   _  0 

I     dp   W  ;  ~  W    dp 
sp    c   p     p    rc 


,  -  VL,      VJ    d  I 

1  d  ,_M>       M       sp  _  n 

I  dp  lW  ;  ~   2~  dp 

sp  c  p  W  I      *c 

v  *  p  sp 


From  equation  (8) 


(9) 


wu   w  WXT 

Ji  =:  J=  +  -N  +  1                                  (10) 

w    w  w                                *  u' 

p    p  p 


Substituting  equation  (10)  into  equation  (9) , 


,  dW  W   dW  W„  dW       W           d  I 

1  C    C    p    N    p       M        sp    _  q 

W  dp  "  T7  2  dp  ~   TT  2  dp  "   W   I     dp 

p   ^c  W     ^c  W     ^c     p   sp     c 

^  p  p       *      * 
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Substituting  equations  (12)  and  (13)  into  equation  (11) 


Y-l 
2.4p   W       W  RTn      p    y 

c   P E Q (12) 

pa      g   I  A   p    p 
p   y     rc   sp   rc   rc 


Substituting  for  g   I     in  the  above  equation  yields , 


Y-l 
2Y 


(f)  -i 

pc 


p   y 


(14) 


If  a  propellant  and  casing  material  are  chosen,  equation  (14) 

can  be  solved  for  the  chamber  pressures,  p  ,  which  minimizes 

W„/I  .   One  interesting  feature  of  equation  (14)  is  that  there 
M  ^ 

is  a  minimum  yield  strength  to  density  ration  -%-   that  will  yield 

pc 

a  usable  chamber  pressure,  furthermore  this  optimization  yields 
a  very  shallow  minimum.   Therefore,  if  the  thrust  or  exit  dia- 
meter requires  an  increase  in  p   the  penalty  paid  in  additional 
casing  weight  will  be  small. 

D.   TYPICAL  PROPELLANT S 

From  the  preceding  discussion  it  can  be  seen  that  some  know- 
ledge of  propellants  to  be  used  in  the  rocket  motor  is  needed. 
The  thrust  developed  by  a  rocket  motor  depends  directly  upon  the 
pressure  in  the  chamber. 


T  =  Cp  AT  Pc 
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Since 


PI  I 

r  sp  sp 


dW  I         dl  W„    dl 

C    +      (w  }     _J3£  SJO  .  M  SJO    m 

dp  c  N      _      2      dp  I  dp 

rc  I  *c  sp        ^c 

sp  r 


From  equation  <8) 


w   +  w„  =  w„  -  w 

C      N     M     p 


Substituting   for  W     +   W„ 

c  N 


dW  W      dl 

_£. E        SP   =   o 

dp  I  dp  u 

rc  sp        ^c 


From  equation    (6) 


dW  2.4iTr   A    p 
c  _  c   c 

dp      -  a 

*c  y 


Since 


dp  p      a 

c  Kp      y 


from  equation    (3) 


(11) 


W      =    irr    I    p 
P  c   p 


dW  2.4ttp      W 

C  C     ^  (12) 


Y-l 
dl  RT  P         T 

-3E  =    2_    <-°)  (13) 

dp  g    p    I  vp 

^c        rc    c   sp      ^c 
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Because  of  this  relationship,  high  pressures  are  needed  for  the 
boost  phase,  and  relatively  low  pressure  are  needed  for  the  sus- 
tain phase  of  flight.   The  pressure  developed  in  a  rocket  motor 
chamber  is  a  function  of  the  burn  rate  and  burn  surface  area  of 
the  propellant.   In  order  to  provide  more  burn  area  for  the  boost 
phase,  the  grain  normally  has  an  internal  star  perforation;  while 
a  sustain  motor  is  normally  an  end  burning  grain  of  solid  cross- 
section.   The  volumetric  loading  of  a  rocket  motor  is  defined  as 
the  ratio  of  the  propellant  volume  to  the  rocket  motor   chamber 
volume. 

_  Grain  volume 
p   Chamber  volume 

Due  to  erosive  burning  effects  the  volumetric  loading  of  a  booster 
is  normally  limited  to  less  than  0.9.   For  efficient  packing 
the  ratio  varies  from  0.7  to  0.9.   The  volumetric  loading  of  an 
end  burning  sustainer  engine  is  normally  1.0.   The  range  of 
propellant   characteristics  is  shown  in  Table  (8-1). 

TABLE  (8-1).   Typical  Propellant  Properties  [14,15]. 


Sustainer 

Booster 

I    (sec) 
sp 

180-210 

210-260 

p   (lbm/in3) 
P 

.059-. 062 

.062-. 065 

Y 

1.24  -  1.27 

1 

.22  -  1.26 

n 

1.00 

0.85 
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E.   DESIGN  EXAMPLE  (BOOSTER) 

From  the  previous  analysis  the  launch  speed  is  M  =  1.5,  and 
the  cruise  speed  is  M  =  3.0  at  an  altitude  of  10,000  feet.  The 
boost  incremental  velocity  is  then, 


AV  =  1616.10  ft/sec 
If  the  maximum  acceleration  during  boost  is  a  constant  30  gs 


(assume  sea  level  g  ),  the  time  of  burn  is, 


ATT    ,                  1616.10    ft/sec 
t      =    AV  /a  =   - 0 

b      B   max    30(32.2  ft/sec2) 


t,  =  1.6  7  sec 


From  equation  (4)  with  yt  =  0 

Li 


«L     AVB 


1     rr-  =  -  - 

n  W_    g  I 

E    *   sp 


Assume  I   =  250  sec  from  Table  (I) 
sp 

WL 
I      ~  =    0.2008 
n  WE 

WE 
p     L       W_ 

Li 


For 


W_  =  794.83 

Li 


W   =  144.57  £bm 
P 
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The  thrust  provided  is 


I    W 

(T),  =  -2£ — E.  =  21,642  £b- 
b      tb  f 


Since  the  missile  flies  essentially  at  zero  angle  of 
attack  during  boost,  the  drag  at  the  end  of  boost  is, 


D  "  (CD0'  «  Sref 


D  =  2006  Jib- 


Therefore  if  an  average  drag  of  1000  lbs  is  assumed,  equation  (3) 
can  be  used  to  find  a  new  W./W^ 

Jj  ti 


WT  AV„    W 

I±  -  B       p 


n   WE         g(T    -    D)^ 


WL 

-^  =    .2270 

WE 


W      =    161.40    lb 
p  m 

This  gives  a  thrust  of  24,162  lbs.   One  more  iteration  of 
equation  (4)  gives, 

W   =  159-46  £& 
p  m 

(T)   =  23,872  Jib, 
s  f 
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Continued  iteration  of  equation  (4)  does  not  change  the  pro- 

pellant  weight. 

3 

From  Table  (I);   p   =  .062  lb   An 

p  m 

n   =  .85 
P 

The  booster  volume  then  becomes 
VD  =  3025.81  in3 

a 

Allowing  one-half  inch  for  structure,  the  booster  length  becomes. 

4  VB 
L_  =  5-  =   42.69  inches 

ttD 
c 

F.   DESIGN  EXAMPLE  (SUSTAINER) 

The  maximum  range  of  the  missile  was  determined  to  be  35  nmiles 

R    =  212,800  ft 
max 

The  distance  traveled  during  boost  is 


s„  =  vM  t, 

B     Mb 


-        Vl  +  V2 

V  =  -=-* =  2424.15  ft/sec 


SD  =  4048.33  ft 

Therefore,  at  a  cruise  Mach  number  of  3.0,  the  sustainer  time 
of  burn  is 

t   =  64.59  sec. 

s 
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The  cruise  drag  consists  of  zero-lift  drag  and  induced  drag 
The  lift  coefficient  of  the  missile  minus  the  booster  is; 


c.  =   w 


'L    q  S 


ref 


CT  =  0.1270 


From  Chapter  6 


Ctv  =  C„   -(-  CT  2/C.. 
&  D-     L    Na 


CD  =  .4016 


The  cruise  drag  then  becomes 
D  =  2009  £bf 

The  total  impulse  required  is  then; 


(Im)   =  Dt   =  129,761  £b£  -  sec 
T  s      s  f 


From  Table  (I)  the  specific  impulse  is  210  sec;  therefore, 
the  weight  of  propellant  needed  for  cruise  is 

(Vs 
(W  )   =   ./  s  =  618  £b 
p  s     I  m 

e  sp 

The  volume  required  is 

V   =  9967.7  in3 
s 
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The  sustainer  length  becomes, 

4V 


L   =  1  =  140.62  inches 

S        ~   * 

ttD 
c 


The  total  rocket  motor  length  is,  L_„; 

3         RM 


LRM  =  LB  +  LS 


L„M  =  183.31  in 
RM 


The  rocket  motor  alone  exeeds  the  total  length  of  the  base- 
line missile.   This  will  lead  to  large  length  to  diameter  ratios 
for  the  entire  missile;  therefore,  the  rocket  motor  length  must 
be  decreased.   That  is  the  subject  of  the  next  section. 

G.   REDUCING  ROCKET  MOTOR  LENGTH 

If  the  rocket  motor  length  of  the  last  section  leads  to 
excessive  length  to  diameter  ratios  for  the  missile,  the  motor 
length  must  be  decreased.   This  may  be  accomplished  in  one  of 
two  ways:   (1)  If  the  mission  profile  has  a  long  cruise  segment, 
the  cruise  altitude  may  be  increased  thereby  decreasing  the  drag 
and  the  total  impulse  required.   (2)  If  the  cruise  altitude 
cannot  be  varied,  the  missile  diameter  must  be  resized  for  pro- 
pulsion considerations.   In  this  design  example  option  (2)  will 
be  selected. 

The  missile  length  can  be  expressed  as 

L  -  Ln  +  LG  +  Lc  +  LWH  +  LRM  (15) 
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L   =  Nose  length 

L   =  Guidance  section  length 

L   =  Control  section  length 
c  3 

LT,TT  =  Warhead  section  length 
WH  3 

I»««  =  Rocket  motor  length 
RM  3 


The  nose  and  warhead  sections  have  been  previously  defined. 
Since  the  propellant  weight  for  the  range  requirement  is  known, 
the  diameter  necessary  can  be  determined  if  a  maximum  missile 
length  is  specified.   The  drag  will  increase  slightly  due  to  the 
increased  diameter  but  the  total  impulse,  and  therefore  the  pro- 
pellant weight,  will  change  only  slightly. 

The  rocket  motor  length  can  be  expressed  in  terms  of  the 

diameter  by 

4W 
L    =  E 


RM       n2 
p  ttD 

P 

Equation  (15)  then  becomes 

4Wp 
(L)  =  Ln  +  LG  +  Lc  +  LWH  +  — I 

0  TTD 

P 

The  guidance  and  the  control  sections  are  each  normally  about 

10  percent  of  the  missile  length. 

L  LT7H  4W 

(L)    =    (— )D    +    .11*   +    .1L    +     (-g±i)     D    +    — 2-j 

D  p    ttD 


L  L 

_£  WH 

D  D 


,  9  4W 

DJ    -    .8LD      +    — 2_  =    o  (16) 

P     TT 
P 
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As  mentioned  in  Chapter  2  the  missile  is  normally  designed  for 

an  existing  launcher.   This  launcher  will  have  a  maximum  length 

that  it  can  accomodate,  L    .   If  this  value  is  substituted  into 

max 

equation  (16),  the  required  diameter  can  be  found. 

1.   Design  Example  (Resizing  for  Propulsion) 

From  previous  analysis;   L  /D  =  2.2  3 

n 

WD  -  2-50 

The  propellant  weight  from  before  was,   W   =  777.46  lbs.   If 
the  maximum  length  of  the  missile  is  taken  as  210  inches,  equa- 
tion (16)  becomes, 

4.73D3  -  168D2  +  15,467.08  =  0.0 

The  required  diameter  is  then, 

(D)     =  11.75  inches 
req 

Allowing  0.25  in  for  structure,  the  missile  diameter  becomes 
12  inches.   From  the  equations  for  the  rocket  motor  lengths, 

L,  =  26.62  in 

L   =  90.47  in 

s 

L^M  =  117.09  in 

RM 

H.   DESIGN  EXAMPLE  (CHAMBER  PRESSURE) 

From  equation  (14),  the  strength  to  density  ratio  can  be 

solved  for  in  terms  of  p  . 

^c 
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a 


_  JlL 


(y-i) 


mm 


1-Y 

P    Y 

(P2) 


-1 


,2.4, 

P 


If  the  minimum  usable  chamber  pressure  is  P   =  500  psi,  and 
using  the  propellant  previously  selected. 


Y  =  1.24 


p   =  .062 

P 

P0  =  io.il 


the  minimum  strength  to  density  required  can  be  found 


a 


=  225,542.87  in 


min 


Inconnel  718  is  a  common  metal  used  in  combustion  chambers 
and  has  the  following  properties  [4]: 


a      =  180,000  psi 


P_  = 


.2662   b  /in' 
m 


y  _ 


676,183  >  (-£) 

pc 


min 


The  qptimum  chamber  pressure  is  then  froit  equation  (14) 


p   =  1132.55  psi 


The  rocket  motor  chamber  wall  thickness  is  given  by  equation  (5) 
as 
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t  = 


p    r 


t  =  0.05  inches 


The  wall  weight  is  given  by  equation  (2) 


W  =  2tt  r£„„tp 
c        RM  c 


W  =  59.5  lb 
c  m 


I.   DESIGN  EXAMPLE  (SIZING  THE  NOZZLE) 

For  initial  analysis  it  is  assumed  that  the  flow  is  fully 
expanded,  in  a  nozzle  with  a  half  angle  of  15  degrees.   With 
the  chamber  pressure  of  1132.55  psi  the  thrust  coefficient  at 
altitude  is 


Y+l 

2y2         2      y_1 
C_      =   C,A\/    £t-    (_£_) 
Fd  d     V     y-1      Y+l 


Y-l 

P0      Y 
1    -     (— ) 


C         =    2.3023 


The  throat  area  and  diameter  then  become 


A.  =  =^-  p   =  9.1552  in2 
t    C„   ^c 


d   =  3.4142  in 
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From  isentropic  tables  for  Pn/p   =  .0089 


gi=  (\)     =  .3342 


The  exit  diameter  becomes, 


d   =  10.22  in 

e 


The  nozzle  length  for  a  15  degree  half  angle  is, 


I      =  12.70  in 
n 


The  nozzle  is  shown  in  Figure  (8-2).   Figure  (8-3)  is  the  design 
with  resized  diameter  and  length  for  propulsion  and  the  optimum 
wings  determined  in  Chapter  7. 

J.   TAIL  SIZING 

As  mentioned  in  Chapter  7  the  tail  (canard)  sizing  is  accom- 
plished for  zero  static  stability  at  launch.   Now  that  the  missile 
length  has  been  fixed,  the  canard  can  be  sized.   Figure  (8-3) 
is  a  drawing  of  the  design  to  this  point.   The  total  lifting  sur- 
face required  is 

STT  +  S   =6.57  ft2 
W     c 

"  *  =  Sc/Sw  (16) 


From  equation  (5)  in  Chapter  7,  for  C„   =  0 

Ma 
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210 


s 

0  =  (P   )   +  (C7  )   XT7  5-2L   (K-.     +  Kn    ) 

Na  „      Na  „   W  S   -    W,_»     B  ,.7v 
N        W     ref     (B)      (W) 


(17) 
S. 


+  (CNa)   XC  S^   (KC     +  KB    } 
Na  C   C  Sref    C(B)     B(C) 


Equations  (16)  and  (17)  above  can  be  solved  for  ST7  and  k 

w 

Sw  =  4.70  ft2 


k  =  0.40 


Therefore,    S   =1.87 
c 


If  the  same  aspect  ratio  as  found  in  Chapter  7  is  used,  the  wing 
and  canard  planforms  become, 


Wing:   b   =1.86  ft      Canard:    b   =1.36  ft 
3    w  c 

(c  )   =  5.05  ft  (c  )   =  2.75  ft 

r  r 

w  c 

(c)   =  3.60  ft  (c)   =  2.63  ft 

w  c 
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IX.   NONLINEAR  AERODYNAMICS  AND  AEROl 

As  mentioned  in  Chpater  5  the  analysis  thus  far  has  been  for 
one  speed  and  small  angles  of  attack.   The  aerodynamics  change 
dramatically  from  subsonic  to  supersonic  flow  and  with  increasing 
angle  of  attack.   This  chapter  presents  the  methods  used  primarily 
by  USAF  Stability  and  Control  DATCOM  and  covers  all  configura- 
tions and  flight  speeds.   As  will  be  shown,  the  method  becomes 
very  involved  and  therefore  lends  itself  to  programming  on  a 
digital  computer. 

The  lifting  characteristics  of  both  wings  and  bodies  become 
nonlinear  as  the  missile  angle  of  attack  increases  above  10  . 
Up  to  this  point  in  the  thesis  development  only  the  linear 
contribution  has  been  considered.   For  large  angles  of  attack, 
the  nonlinear  term,  which  is  due  to  flow  separation,  must  be 
considered.   The  relative  effect  of  these  terms  on  the  normal 
force  and  pitching  moment  coefficients  is  shown  in  Figure  (9-1) . 

A.   VISCOUS  CROSS-FLOW 

As  can  be  seen  in  Figure  (9-1) ,  at  large  angles  of  attack 

the  forces  acting  are  primarily  nonlinear.   The  nonlinear  term 

is  normally  described  through  the  use  of  the  viscous  cross-flow 

coefficient,  C, 
d 
c 

f 

C,   =  — ^—  S   =  planform  area, 

d     q  S  P 

c    Mn   p  ^ 
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Figure  (9-1) .   Relative  contributions  of  linear  and  nonlinear 

terms  [19 ] . 


For  an  infinitely  long  circular  cylinder,  the  viscous  cross- 
flow  force  per  unit  length  is 


f   =  2r  C ,   q 
v       d    n 
c 


(1) 


The  viscous  cross-flow  coefficient  is  determined  experimentally 
and  is  a  function  of  normal  Mach  number  and  Reynolds  number. 
Figure  (9-2)  gives  the  cross-flow  drag  coefficient  as  a  function 
of  cross-flow  Mach  number. 
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M   =  Mw  sin  a 
n    M 

R    =  R   sin  a 
e     e 
n 

V   =  VM  sin  a 

n    M 

Since  in  real  flow  there  is  spillage  around  the  ends  of  a  finite 
length  cylinder,  the  drag  is  less  than  that  for  an  infinite,  2-D 
cylinder.   The  cross-flow  drag  proportionality  constant,  n,  is 
the  ratio  of  the  drag  coefficient  of  a  finite  cylinder  to  the 
drag  coefficient  of  an  infinite  cylinder. 

Cdc 

n  =  (c,  )° 

dC  2-D 

The  proportionality  constant  is  also  determined  experimentally 
and  is  given  in  Figure  (9-3) .   Notice  that  the  cross-flow  drag 
proportionality  constant  approaches  one  as  the  length  to  dia- 
meter ratio  becomes  large,  or  the  2-D  situation  is  approached. 
The  viscous  cross-flow  force  per  unit  length  for  a  finite 
cylinder  then  becomes 

pv2 

fv  -  2rncd  -r-  (2) 

c 


Since 


pVn2    pVM2   .2        .2 
— 7T-   =  — z—  sin  a  =  q  sin  a 


equation  (2)  becomes, 


2 
f   =  2nC,  q  r  sin  a 
v      d  ^ 
c 
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The  viscous  contribution  to  the  normal  force  acting  on  the  cylinder 

then  becomes 

2 

F        2nCd  sin  a  ~l 

(cn'    ;  ^t—  -      sc  :      /    rdx 

y   ^   ref       ref    JQ 

For  a  uniform  cylinder, 


s     2 

(C )   =  n  C ,   =  P   sin  a 

N  TT       d   S   - 
V       c   ref 


If  this  force  is  added  to  the  term  predicted  by  slender  body 
potential  theory,  the  total  normal  force  acting  on  the  body  can 
be  found.   Reference  [19]  gives  this  total  force  as; 

i~    \       b  .  2       a  „      d    .  2 

(CLT)   = sin  a  cos  ■=■  +  n  C ,   =  K   sin  a 

N_S£  2       d   S   ^ 
B    ref                  c   ref 

Similar  terms  can  be  developed  for  the  wing  and  tail  since 
they  also  exhibit  nonlinear  behavior  at  high  angles  of  attack. 
Care  must  be  taken  to  separate  the  lift  of  Chapter  5  into  the 
lift  acting  on  the  wing  panel  and  the  additional  lift  on  the  body 
before  the  cross-flow  term  is  applied.   The  nonlinear  cross-flow 
component  of  wing  lift  is  caused  by  flow  separation  and  the  for- 
mation of  vortices  on  the  upper  surface  of  the  wing.   This 

viscous  term  is  given  by,  C 

NW(V) 

C„     =  C^,   sin2a 


NT7,„.     d         S   - 
W(V)      n        ref 


Where  C,   is  the  cross-flow  drag  coefficient  for  a  flat  plate  and 
n 


is  given  in  Figure  (9-4)  . 
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B.   TOTAL  MISSILE  LIFT 

The  remainder  of  this  chapter  is  a  summary  of  references  [17] 
and  [18].   With  the  addition  of  the  force  due  to  viscous  cross- 
flow,  the  total  lift  on  the  body  becomes 


C__   =  (CM)   +  CM     +  C„       +  C„ 
NB      N  B     NB(W)     NB(T)-a     NB(T)-6 


where, 

(C\J   =  linear  and  nonlinear  lift  force  on  the  body 
N  B 

C  .     =  additional  lift  on  the  body  due  to  the  wing 
NB(W) 

C        =  additional  lift  on  the  body  due  to  the  tail 
NB(T)-a 

(due  to  angle  of  attack) 

C        =  additional  lift  on  the  body  due  to  the  tail 
NB(T)-5 

(due  to  control  deflection) 

The  additional  lift  on  the  body  in  the  presence  of  the  wing  and 
tail  can  be  found  as  described  in  Chapter  6 . 

The  wing  lift  is  now  composed  of  a  linear  and  a  nonlinear 
component  and  becomes 


NW  NW(B)  NW(V) 

s  s 

CNM   "    Vb)      <CWb   S-^T   Sina  +    Cd       Sin    °    S^- 

W                                  W      ref  n                    ref 


The  tail  lift  is  found  in  a  similar  manner  only  with  the 
additional  term  for  any  control  surface  deflection,  6. 


C    =  C        +  C        +  C 
NT     NT(B)-a     NT(B)-6     NT  (V) 
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Notice  that  the  linear  terms  in  the  above  equations  are  multiplied 

by  sina   instead  of  a  as  in  Chapter  6.   The  linear  theory  from 

which  the  Chapter  6  equations  were  derived  assumes  small  angles 

of  attack  and  therefore  sina=:a  .   This  is  not  true  at  higher 

angles  of  attack  and  is  therefore  included  here. 

The  total  normal  force  then  becomes ,  C„ 

N 

C=C    +  C    +C    +C 
N     NB     Nw     NT     L. 

The  last  term,  CT  ,  is  the  lift- loss  due  to  downwash  and  is  given 

J_i  • 

1 

by 

^wVT[\(B)^-kw(B)^]i  (»-r)TSw 

i  2tt  AR_  (f   -  r  )  S   . 

T   w    w   ref 


where, 


i  =  interference  factor 


f   =  votex  location 
w 


the  above  terms  are  found  by  the  methods  of  reference  [13]. 

C.   DRAG  CHARACTERISTICS  OF  A  MISSILE 

The  total  drag  acting  on  a  missile  is  the  sum  of  the  zero- 
lift  drag  and  the  induced  drag  (due  to  angle  of  attack  and/or 
control  surface  deflection) .   The  zero  lift  drag  of  bodies  and 
wings  is  highly  dependent  on  the  missile  speed. 
Three  speed  regimes  are  normally  considered. 
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Subsonic  ►  0  <  Mw  <  .8 

M 

Transonic  ►  .8  <  M   <  1.2 

-   M  - 

Supersonic—; — *  M  >  1.2 

Because  of  the  empirical  nature  of  the  formulas  for  zero-lift 
drag,  the  procedures  followed  are  essentially  those  of  DATCOM  [16] 
1.   Zero  Lift  Drag 
a.   Subsonic 

At  subsonic  speeds  the  zero-lift  drag  consists  of 
skin  friction  (incompressible)  and  pressure  drag.   The  pressure 
drag  at  subsonic  speeds  is  usually  small  compared  to  the  drag 
due  to  skin  friction.   The  entire  zero-lift  drag  of  a  missile 
at  subsonic  speeds  is  given  by; 

C    =  C     +  C     +  C 
D0     °0B     D0W     D0T 

Where  the  components  are  the  body,  wing  and  tail  contributions 


(1)   Body  drag,  C    .   The  body  zero-lift  drag  is 

OB 


given   by; 


=   1.02    :       ;      1  +  -A7I  +  -4 

D0B  fBl  fV2         f3 


+  c 


S  4=  AT3 

ref  B 


where  C    is  found  by  the  methods  of  Chapter  5. 
AB 
The  first  term  is  the  skin  friction  contribution  and  the 

next  two  terms  are  the  pressure  contributions.   f  is  the  body 

fineness  ratio  and  is  given  by; 
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(2)   Wing  drag,  C    .   The  wing  zero-lift  drag  is 

OW 


given  by , 


C     =  C 
OW     w 


4 

1  +  2(-)  +  100 (-) 
c        c 


Wet 


w 


ref 


(SrT  .)   is  the  entire  wetted  area  of  the  wings.   The  thickness 
Wet  Ty 
W 

to  chord  ratio  is  given  by  (t/c) . 


(3)   Tail  drag,  C    .   The  tail  zero-lift  drag  is 

D0T 


given  by, 


CD    "  Cf 

0T      T 


1  +  2(-)  +  100 (-) 
c        c 


Wet 


'ref 


The  above  analysis  assumes  fully  turbulent  flow  along  all  sur- 
faces.  C_  ,  C_   and  C-   are  the  local  average  skin  friction 

tB    If        T 
coefficients  based  on  the  local  Reynolds  number,  R  . 

The  reference  lengths  are  the  body  length  and 
wing/tail  mean  aerodynamic  chords.   The  skin  friction  coefficient 
is  then  given  by 


cf  - 


4.55 


(log1Q  R  ) 


2.58 


b.   Transonic 

The  total  transonic  zero-lift  drag  is  composed  of 
skin-friction  drag,  transonic  wave  drag,  pressure  drag  and 
base  pressure  drag.   The  compressible  skin  friction  drag  is 
found  by  using  a  correction  factor  on  the  incompressible  co- 
efficient, C-  ,  found  for  subsonic  flow.   The  compressibility 
B 
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correction  is  found  by  the  methods  of  Chapter  6.   The  skin 
friction  drag  then  becomes 

(WB 

CD    =  1'02  Cf   — S " 

UfB         rc    ref 

The  subsonic  pressure  drag  is  the  same  as  before 


1.5  7 

PB  ^B     L?71  ?j 


CD         =    1.02    Cf 


Wet'B 
ref 


The  above  equation  applies  for  Mach  numbers  in  the  range  of 

.8  to  1.0.   The  pressure  drag  then  decreases  linearly  to  zero 

at  a  Mach  number  of  1.2.   The  transonic  wave  drag  C    is  deter- 

rB 
mined  using  Figure  (9-5) ,  which  is  a  function  of  nose  fineness 

rati°  LN/D- 

The  total  transonic  zero-lift  drag  of  the  body  then  becomes 

SN 

c    =  c   +  c    +  c    — - —  +  c 
D0B    Df    DPB    °VB  Sref    AB 

where  S„  is  the  cross  sectional  area  at  the  nose-body  junction. 

The  total  transonic  zero-lift  drag  of  the  aerodynamic  surfaces 

is  composed  of  the  skin  friction  drag  and  a  drag  increment, 

AC   ,  which  is  the  transonic  wave  drag. 
0 

Experimental  results  show  little  increase  in  the  viscous 
drag  of  aerodynamic  surfaces  from  the  subsonic  to  the  transonic 
regimes;  therefore  the  skin  friction  and  pressure  drag  is  the 
same  as  for  subsonic  flow  and  is  given  by 
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c    =  c 

D0W     fw 


t         t  4 
1  +  2(-)  +  100  (-) 
c        c 


SWet 
ref 


The  wave  drag  of  the  aerodynamic  surfaces  if  found  from 

Figure  (9-6)  and  is  a  function  of  — ,  AR  and  M.   For  swept  wings 

the  Mach  number  used  in  Figure  (9-6)  is  given  by, 

M»    =   M    [cos  A  c/4]1//2 

The  transonic  wave  drag  increment  is  then  given  by 

2.5    S 
AC     =  AC  '    [cos  A  c/4] 


D0W     D0W  Sref 

where   AC    is  obtained  from  Figure  (9-6) .    The  tail  zero-lift 

D0W 
drag  is  found  in  the  same  manner  as  the  wing.   The  total  zero-lift 

drag  (transonic)  is  then  given  by; 

D0     D0B     D0W      D0W     D0T      D0T 


c.   Supersonic 

The  supersonic  zero-lift  drag  of  a  missile  is  deter- 
mined empirically  by  assuming  a  parabolic  variation  of  C   with 
Mach  number  between  1.2  and  3.0.   The  resulting  equation  is  for 
the  entire  missile  and  is  given  by 


C         =    C    ■     -    1.7209     (C    "    -    C     •)     +    1.5708     (C    "    -    CM     /m 
0  0  U0  0  0  0 
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Figure  (9-5)  .   Transonic  wave  drag  for  ogival  and  blunted 

conical  forebodies  [18]. 
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where 


U        U  M  =  1.2 


U        U  M  =  3.0 


C  '   is  determined  by  using  the  methods  for  transonic  flow. 

The  magnitude  of  the  supersonic  wave  drag  is  highly 

dependent  upon  the  nose  shape.   For  this  reason  two  methods 

are  used  to  determine  C  " 

0 
Method  I;  This  method  is  for  blunted  ogives,  pointed  ogives  and 

blunted  cones.   In  this  method 


F°r  LN/D  <  '5  %    '    Cd'q 

For  L.T/n  >  8.0        C  "  =  (C_  )        +  AC_ 

7  "  o      0  M         Dow 

M  =  .  o 


+  AC     +  C 

D0T     AB 


The  above  values  are  determined  from  the  methods  of  transonic 
flow. 


For  L.. ,_  =  .5  to  8.0 
N/D 

C  "   =  K,  C  ' 

D0     x  Do 

K,  in  the  above   equation  is  derived  empirically  and  is  given  in 
Figure  (9-7) . 

Method  II:   This  method  is  for  pointed  conical  noses.   In  this 
method  C  "  is  determined  by 
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C  "   =  (C_  )        +  AC_    +  AC_    +  (C_   )       +  C 

D0M=.8      °0W      D°T      D™„=3     ^ 


The  first  three  terms  in  the  above  equation  are  found  by  transonic 

flow  methods.   The  forebody  wave  drag,  C    ,  for  M  =  3.0  is  found 

W 
using  Figure  (9-8).   The  nose  semi-vertex  angle,  0  ,  is  the  same 

as  a  in  Chapter  6.   In  all  flow  regimes  the  base  pressure  drag 

is  found  as  in  Chapter  6 . 

2.   Induced  Drag 

The  induced  drag  due  to  angle  of  attack  depends  upon  the 

flight  regime  the  missile  is  in.   For  M  <  .85  and  AR  >  3.0  the 

induced  drag  is, 

<cl'2 

C 


D .    TrARe 


where  e  is  the  Oswald  efficiency  factor  =.  0.7. 

The  tactical  missile  normally  has  an  aspect  ratio  of  less  than 

3.0.   For  all  flight  speeds  with  AR  <  3.0  the  induced  drag  is 

C„   =  CT  tan  a 
D.     L 

l 

The  component  induced  drag  coefficients  are  found  in  the  same 
manner  as  above  using  component  lift  coefficients. 
3.   Total  Drag 

From  the  zero-lift  drag  and  induced  drag  the  missile 
total  drag  becomes, 


D     DQ     D. 
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D.   PITCHING  MOMENT  CHARACTERISTICS 

The  pitching  moment  of  a  missile  is  equal  to  the  sum  of 
the  pitching  moments  due  to  lift  and  drag  forces  acting  on  the 
body,  wings  and  tails.   If  only  small  angles  of  attack  are  of 
interest,  the  pitching  moment  is  due  primarily  to  the  lift 
forces.   The  methods  presented  here  are  for  all  angles  of  attack. 

1.   Body 

The  body  alone  pitching  moment  about  its  center  of  gravity 
is  given  from  slender  body  theory  and  viscous  cross- flow  theory 


as 


"M 


B 


VB  -  SB(L  -  *cg' 

ref  ref 


.2      a 

sin  a  cos  *■  +  n  C , 

2       d    S 

c    ref 


S 

(7T-E-)  ( 


-  x     ~ 
cG     d»  .2 


-£)sin 


li 


'ref 


S   = 


x  = 


VB  = 

SB  " 


Planform  area 

Centroid  of  planform  area 

Body  volume 

Base  area 


With  C   given  the  center  of  pressure  of  the  body  becomes 


M 


B 


x 


(xc  )   - 

P  B 


cG 


M 


B 


ref 


'N 


'ref 


B 


2 .   Wing  (Fixed  Surface) 

The  center  of  pressure  location  for  the  wing  must  be 
found  before  the  pitching  moment  can  be  specified.   The 
center  of  pressure  of  the  various  lift  components  are  found 
by  the  method  of  Pitts,  Nielsen  and  Faatiari  [13].   The  lift 
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components  are  the  same  as  in  the  lift  section.   The  center 

of  pressure  of  the  lift  on  the  wing  in  the  presence  of  the 

body,  (x)__    ,  is  found  from  Figure  (9-10)  if  the  flow  is 

W(B) 
subsonic  and  Figure  (9-11)  if  the  flow  is  supersonic. 

The  center  of  pressure  of  the  additional  lift  on  the  body 

in  the  presence  of  the  wing,  (x)     ,  is  found  using  Figure 

B  (W) 
(9-12)  if  the  flow  is  subsonic.   If  the  flow  is  supersonic 

the  center  of  pressure  is  found  from  Figure  (9-13)  or  Figure 

(9-14)  depending  on  the  parameters; 

3AR(1  +  X)   (1  +  -i) 

mp 

X  =  taper  ratio 

m  =  co  +  A 

If  the  center  of  pressure  reference  is  moved  to  the  nose,  the 
following  expressions  result. 


(x   ) 


=  t-T^)  (C  )   +  x 


cp  W._.    "c  '- ,      r'  T7    w 
*       (B)      r  W,B,      W 


(x   ) 


.  X 


=  (— )      (C  )   +  x 


cp  B  ,„%     *C._'  - 


(W) 


'r  B(W)    '  w 


cp 


B(W) 


CP  W(B) 


AC 


D 


W-a 
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Figure  (9-10) . 


Wing  center  of  pressure 
(subsonic)  [18] . 
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Wing  center  of  pressure 
(supersonic)  [18]. 
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Figure  (9-12) .   Body  center  of  pressure  [18] 
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Figure  (9-13).   Body  center  of  pressure  [18]. 

237 


CD 


Ix 


U 


en 
O 


CO 


Ix 


CO 


U 


O 


1.6 
\2 

0.8 

0.4 

0.0 
1.6 

1.2 
0.8 


CD 


l«|u"       °-4 

0.0 
2.0 

1.6 


CD 


£  I  <-£  1 .2 


o 


CO 


'"I 


(~L 


0.8 


0.4 


0.0 


A -0.0 

.  r   b  -0.6 

1/ 

^^ 

V 

r/b  -0.4 

^^** 

/yV^*- "* 

rUdA                      r/i,-ni 

A  =  0.5 

/  ''S\**^~~ 

r 

_   /L            rt     *                                                                       r/U    _    ft    /< 

X-1.0 



^^^^^ 

1  /jS"^ 

is 

2  4  6 

EFFECTIVE  ASPECT  RATIO,  /MR 


r/b  -  0.2 


r/b  -  0.0 


r/b -0.2 


r/b  -  0.0 


r/b  -  0.2 


r/b  -  0.0 


Figure  (9-14) .   Body  center  of  pressure  [18] . 
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From  the  above  figure  the  entire  pitching  moment  about 
the  nose  of  the  missile  due  to  the  forces  acting  on  the  wing  is, 


CM   =  [  (C_      +  C_      +  CT  )cosa  +  AC      sina]  (■ 


^l 


'"w    L^LW(B)     LW(V)     Li  'D" 


W-a  ref 


(xc  ) 

+  EC.     cosa]  (   -   B(W)) 

^BCW)  ref 


The  viscous  lift,  C      and  lift-loss  due  to  downwash, 

LW(V) 
C   ,  are  not  shown  in  the  figure  but  are  the  same  forces  as 

Li  ■ 

1 

found  in  the  lift  section.   These  forces  are  assumed  to  be 

acting  through  the  center  of  pressure,  (x   ) 

cp  W(B) 
An  average  center  of  pressure  of  the  wing  due  to  all  forces 

acting  on  it  can  now  be  found  as 


.         ,       _  ^^  ref 

CP  „        Tel       +   CT     )  cosa   +    AC  sina 

w  LT,  L.  Dtt 

W      1  uW-a 


where 


C    =  C      +  C       +  C 
LW     LW(B)     LB(W)     LW(V) 


3.   Tail  (Control  Surface) 


The  tail  pitching  moment  about  the  nose  of  the  missile 
is  found  in  the  same  way  except  now  a  control  surface  deflec- 
tion must  be  included.   The  equations  now  become 
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'    =  [  (C_      +  Cl      +  C   )cosa   +  C 
MT       LT(B)      T ,  x      l  LT(B)-5 


(V) 


(x   ) 

CP  T(B) 


(AC_     +  AC_    )sina]  ( 

DT-a     DT-5  Lref 


]  (■ 


(x   ) 
cp 


+  [c_      cosa  +  C 

H3(T)  ijB(T)-6       ref 


B(T)  \ 


C    L   - 
M   ref 

fx    )    a   i 

cp  \_     (CT     +  C_  )cosa  +  CT     +  (AC_     +  AC_    )  sina 

*        "m        Li .  -i-<m  r        D—.         Dr 


T-a 


T-5 


T-a 


'T-6 


where, 


C     =  C      +  C      +  C 
LT-a    LT(B)     LB(T)     LT(V) 


c    =  c      +  c 

LT-6     LT(B)-<5     LB(T)-S 


The  wing  and  tail  pitching  moments  above  were  found  about 
the  nose  of  the  missile.   Transferring  the  pitching  moments 
about  the  center  of  gravity  the  complete  missile  pitching 
moment  becomes 


C   =  C     +  C 

M     MB      MW 


CG      cp 


W 


<x™} 

cp 


w 


+  c 


M. 


x  _  -  (x   ) 

cG      cp 


(x   ) 

cp 


E.   AER01  DESCRIPTION 

It  can  be  seen  from  the  preceding  description  of  the  aero- 
dynamics of  a  missile,  that  the  process  of  obtaining  a  complete 
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description  of  the  aerodynamic  coefficients  becomes  an  involved 
task.   This  was  the  justification  for  initially  using  linear 
aerodynamics.   To  complete  an  analysis  a  fast  method  of  pre- 
dicting the  static  aerodynamic  characteristics  is  needed. 
The  method  needs  to  include  both  linear  and  nonlinear  contri- 
butions as  well  as  interference  factors  and  must  be  applicable 
to  all  speed  regimes.   This  analysis  lends   itself  to  program- 
ming on  a  digital  computer.   AER01  is  a  modification  of  the 
program  AEROCF  which  was  developed  at  the  Naval  Air  Development 
Center  by  Mr.  F.A.  Kuster,  Jr.  [17] .   The  program  is  essentially 
as  written  except  for  calculation  of  planform  areas  and  centroids. 
The  program  was  also  modified  for  use  on  the  Naval  Postgraduate 
School  CP/CMS  system. 

AER01  consists  of  a  main  program  and  three  subroutines.   The 
inputs  to  the  program  are  the  geometric  characteristics  of  the 
missile,  the  flight  conditions,  engine  and  inlet  type  and  the 
protuberance  drag.   The  output  consists  of  the  aerodynamic  co- 
efficients for  lift,  drag  and  moment.   The  component  contribution 
to  these  coefficients  are  also  given.   The  component  and  overall 
center  of  pressure  are  also  determined. 

Subroutine  GEOSUB;  This  subroutine  calculates  the  missile 
wetted  area  and  the  Reynolds  number  per  foot  based  on  the  flight 
altitude . 

Subroutine  CLASUB;  This  subroutine  calculates  the  aerodynamic 
surface   lift-curve  slopes. 
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Subroutine  CATSUB;  This  subroutine  calculates  center  of 
pressure  locations,  cross-flow  drag  coefficients,  and  interference 
factors . 

Main  Program;  The  main  program  assumes  the  control  surface 
is  the  tail.   This  is  regardless  of  the  method  of  control  (Wing, 
Tail,  Canard).   Because  of  this,  care  must  be  taken  to  input  the 
right  data  for  the  tail.   For  instance,  if  the  missile  is  a  wing 
control  missile,  the  wing  data  is  input  as  the  tail  and  the 
tail  data  as  the  wing.   Figure  (9-15)  and  Table  (9-1)  give  a 
complete  listing  of  the  input  data.   Table  (9-II)  is  a  list  of 
the  output  data.   Appendix  A  is  a  listing  of  the  program  as  modi- 
fied for  use  on  the  Naval  Postgraduate  School  IBM  360  computer. 

1.   Verification  of  AEROl 

Before  using  the  program  an  attempt  was  made  to  verify 
its  prediction  techniques  and  find  any  pitfalls  in  its  use.   To 
accomplish  this  the  program  was  run  for  various  configurations 
for  which  experimental  data  were  available,  and  the  results  were 
compared.   The  comparisons  are  shown  in  Figures  (9-16)  to  (9-25) 
from  references  [19]  -  [22]  which  are  NASA  technical  notes  and 
memoranda  which  report  results  of  wind  tunnel  tests  for  various 
body-wing-tail  combinations.   In  Figures  (9-16)  to  (9-25)  the 
solid  lines  are  AEROl  predictions. 
a.   NASA  TN  D-6996 

This  technical  note  presents  a  method  of  predicting 
aerodynamic  characteristics  for  bodies  alone  at  angles  of  attack 
from  0  to  180  degrees.   Several  nose-body  combinations  are  given. 
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Variable 


Name 


ICSC 


INOSE 


IDT 


TABLE  (9-1) 
Input  Data 

Format 
12 

12 


12 


4 

IM 

12 

5 

IAL 

12 

6 

NBODY 

12 

7 

ISWPW 

12 

8 

IAFBW 

12 

9 

IWEPW 

12 

10 

NWING 

12 

11 

ISWPT 

12 

12 

IAFBT 

12 

13 

ISWEPT 

12 

14 

NTAIL 

12 

15 

XLAMW 

F10. 5 

16 

CLAMW 

F10. 5 

17 

BW 

F10. 5 

18 

CROOTW 

F10.5 

19 

SW 

F10.  5 

20 

XMACW 

F10. 5 

21 

XWING 

F10.  5 

Meaning 

Type  of  control  (wing, 
tail,  canard) 

Nose  shape  (ellipsoidal, 
cone,  ogive) 

Number  of  control  surface 
deflections 

Numer  of  Mach  numbers 

Number  of  angles  of  attack 

Number  of  configurations 

Wing  shape  (delta, nondelta) 

Missile  body  after  wing 

Leading  edge  sweep  indicator 

Number  of  wings 

Tail  shape  (delta, nondelta) 

Missile  body  after  tail 

Leading  edge  sweep  indicator 

Number  of  tails 

Wing  taper  ratio 

Wing  leading  edge  sweep 

Wing  span 

Wing  root  chard 

Wing  exposed  area 

Wing  mean  geometric  chord 

Distance  to  wing  leading 
edge 
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Variable 

Name 

Format 

22 

TOVCW 

F10.5 

23 

XLAMT 

F10. 5 

24 

CLAMT 

F10.  5 

25 

BT 

F10.  5 

26 

CROOTT 

F10.  5 

27 

ST 

F10.  5 

28 

XMACT 

F10.5 

29 

XTAIL 

F10.  5 

30 

TOVCT 

F10. 5 

31 

HT 

F10. 5 

32 

D 

F10.  5 

33 

XL 

F10. 5 

34 

XLNOSE 

F10.5 

35 

XCG 

F10.5 

36 

AREA 

F10.  5 

37 

XREF 

F10.5 

38 

ENGINE 

F10.5 

39 

INLET 

F10.  5 

40 

BETA 

F10. 5 

41 

DBASE 

F10.  5 

42 

DJET 

F10. 5 

43 

XLABOD 

F10.5 

44 

CDPROT 

F10.  5 

Meaning 

Wing  thickness  to  chord  ratio 

Tail  taper  ratio 

Tail  leading  edge  sweep 

Tail  span 

Tail  root  chord 

Tail  exposed  area 

Tail  mean  geometric  chord 

Distance  to  tail  leading 
edge 

Tail  thickness  to  chord  ratio 

Altitude 

body  diameter 

Body  length 

Nose  length 

Center  of  gravity  location 

Reference  area 

Reference  length 

Engine  code 

Inlet  code 

Boattail  angle 

Base  diameter 

Nozzle  exit  diameter 

Boattail  length 

Proturberance  drag 
coefficient 
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TABLE  (9-II) 
Output  Data 


Variable  Name 
AL 

CLTOT 
CDTOT 
CLWP 
CLBW 

CLTP 
CLBT 

CLB 

CDI 

CNWP 

CNTP 

CLTD 

CDTD 

CN 

CA 

XCPW 

XCPT 

XCP 

CM 


Meaning 

Angle  of  attack 

Total  coefficient  of  lift 

Total  coefficient  of  drag 

Wing  panel  coefficient  of  lift 

Additional  lift  on  body  due  to 
wing  coefficient 

Tail  panel  coefficient  of  lift 

Additional  lift  on  body  due  to 
tail  coefficient 

Body  alone  lift  coefficient 

Induced  drag  coefficient 

Wing  panel  normal  force  coefficient 

Tail  panel  normal  force  coefficient 

Coefficient  of  lift  due  to  tail 
deflection 

Coefficient  of  drag  due  to  tail 
deflection 

Total  normal  force  coefficient 

Total  axial  force  coefficient 

Wing  center  of  pressure 

Tail  center  of  pressure 

Total  missile  center  of  pressure 

Total  pitching  moment  about  center 
of  gravity 
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Figures  (9-16)  to  (9-18)  compare  the  coefficients  predicted  by 
AEROl  with  those  obtained  for  body  number  9  in  the  NASA  report. 
The  normal  force  coefficient  is  predicted  well  throughout  the 
range.  The  moment  and  axial  force  trends  are  predicted  by  the 
program  but  large  errors  exist  throughout  the  range  of  angles 
of  attack. 

b.  NASA  TM  X-236  7 

This  technical  memorandum  investigates  the  aero- 
dynamic characteristics  of  various  cruciform  body-wing  combina- 
tions.  The  coefficients  for  this  configuration  agree  very  well 
with  experimental  values  up  to  10  degrees  angle  of  attack.   Beyond 
this  value  the  lift  and  moment  coefficients  predicted  by  AEROl 
exceed  the  experimental  values  by  as  much  as  10  percent.   Al- 
though the  exact  cause  of  this  error  was  not  investigated,  it  may 
be  partially  explained  by  the  nose  shape  of  the  body.   The  nose 
is  a  combination  ogive  and  cone.   For  purposes  of  AEROl  it  was 
assumed  an  ogive.   The  greater  presented  area  of  the  ogive  would 
contribute  to  a  higher  C   and  C   through  the  cross- flow  terms  in 
these  coefficients.   Figures  (9-19)  through  (9-21)  compare  the 
predicted  with  the  experimentally  determined  coefficients  for 
the  wing-tail  configuration  of  this  reference  at  M  =  .9. 

c.  NASA  TM  X-2780  and  NASA  TM  X-2289 

These  technical  memoranda  investigate  the  aerodynamics 
of  a  delta  wing  missile  using  tail  control  and  a  tail-less  cruci- 
form missile.   As  shown  in  Figures  (9-22)  to  (9-25)  there  is 
excellent  agreement  between  the  experimental  values  of  the  aero- 
dynamic coefficients  and  those  predicted  by  AEROl. 
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F.   COMPONENT  WEIGHTS 

As  with  the  gross  weight,  the  component  weights  can  be 
determined  through  the  use  of  parametric  regression  equations. 
If  the  dimensions  of  the  components  are  known,  the  following 
formulas  can  be  used  to  determine  the  component  weights  [7]. 
Aero  Surface  Weight 

The  weight  of  one  wing  panel  is  given  by, 

WAS  -  6-77483  (EAS)1-02  (ARAS>'56  (1) 

2 
E    Exposed  are  of  one  panel,  ft 

AR_C   Aspect  ratio  of  one  panel 


Body  Structure  Weight 

WBS  -  -0604  (LBS)-64  (Dgg)1-77  (2) 


L    Length  of  body  to  be  covered  by 

85 

structure.   This  does  not  include 
the  rocket  motor  unless  a  separate 
structure  surrounds  the  motor 

casing.   (inches) 

Diameter  of  body  structure 
DBS 

(inches) 


Internal  Systems  Weight 

W        =    .00485     (WG)-74    CLjgJ1-00     (DIS)-42       (3 


W^   Gross  weight  of  the  missile 
G 
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L_q   Length  of  subsystem  (inches) 
D_q   Diameter  of  missile  (inches) 


In  the  above  equation,  the  internal  system  includes  guidance, 
avionics  and  control. 

1.   Design  Example  (Component  Sizing) 

If  the  guidance  and  the  control  sections  are  kept  at 
10  percent  of  the  total  length,  and  the  remaining  components 
are  as  sized  previously,  the  design  is  as  shown  in  Figure  (9-27) 
Since  these  components  are  considered  internal  systems,  their 
weights  can  be  determined  from  equation  (3). 

W    .  =  .00485  (1809)*79  ( . 1L)  (12. 5)*42 
cont 

W     =  W    ,=83. 64  lb 
cont     guid 

The  wing  and  canard  weights  are  found  using  equation  (1)  and 
Figure  (9-26) . 


AR 


AS 


b'  =  b/2 


'AS 


Figure  (9-26)  .   Aero  surface  weight 
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The  wing  weight  is  four  times  the  weight  of  one  panel. 

Wr7  =  4  3.42  lb 
W 

Similarly  the  canard  weight  is 

W  =  9.8  8  lb 

c 

If  the  entire  body  is  covered  by  structure,  the  body  structure 
weight  becomes 

W  .  =  .0604  (232)'64  (12. 5)1*77 
st 

W  .  =  172.38  lb 
st 

The  following  component  table  can  now  be  constructed. 


Component 

Length  (in) 
23.2 

Weight  (lbs) 

Center  of  Gravity  (in) 

Control 

83.64 

39.48 

Guidance 

23.2 

83.64 

62.68 

Warhead 

31.25 

28.26 

89.91 

Sustainer 

87.51 

618.92 

149.29 

Booster 

26.31 

159.46 

206.20 

Motor  Casing 

113.82 

59.50 

162.44 

Wings 

43.42 

211.40 

Canards 

9.88 

44.08 

Structure 

232.00 

172.38 

116.00 
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The  complete  missile  weight  is  then  1264.62  lbs 
The  center  of  gravity  from  the  above  table  is 


x  _  =  139.44  in 
cG 


The  center  of  gravity  at  the  end  of  boost  is 


(xcG)  =  129.85  in 


G.   DESIGN  EXAMPLE  (FINAL  ANALYSIS) 

From  the  design  parameters  so  far,  a  complete  description 
of  the  missile  can  be  determined.   Figure  (9-27)  is  a  drawing 
of  the  missile.   The  launch  conditions  for  the  missile  were 
specified  as, 

MM=  1.5 

h  -  10,000  ft 

W_  =  1264.62 
G 

x  _  =  139.44  in 
cG 

From  these  conditions  the  input  data  for  AER01  is  shown  in 
Table  (9-III) .   The  format  is  the  same  as  that  of  the  output  of 
AEROl  and  is  printed  as  a  check  to  ensure  that  the  input  data 
was  entered  properly.   The  output  is  shown  in  Table  (9-IV) . 
Figure  (9-28)  is  a  plot  of  the  coefficient  of  moment  versus  angle 
of  attack  for  the  launch  condition  and  for  the  beginning  of 
cruise.   The  center  of  gravity  has  moved  forward  approximately 
10  inches  for  the  beginning  of  cruise  so  that  the  stability  has 
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increased  as  indicated  in  this  figure.   The  performance  objective 
from  the  beginning  of  the  design  was  to  provide  a  31  g  maneuver 
capability.   Table  (9-  V)  shows  the  output  of  AER01  for  Mach 
numbers  of  1.5  to  3.0  and  control  deflections  of  0.0  to  30  degrees. 
From  this  output  the  trimmed  normal  force  can  be  found  and  the 
corresponding  load  factor  is  then, 

NTR     ref 

n  =  w 

The  following  values  of  maneuver  load  factor  were  found  using 
AEROl. 


M 

Cntr 

18.80 

q (lb/ft2) 

n(g*s) 

6  req 

2.5 

6369 

85.00 

10  deg. 

2.0 

24.50 

4076 

70.99 

10  deg. 

1.5 

12.00 

2293 

19.56 

10  deg. 

As  indicated  in  the  above  table,  the  missile  meets  the  maneuvering 
specifications  of  the  operational  requirements. 
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X.   CONCLUSIONS  AND  RECOMMENDATIONS 

This  thesis  has  presented  the  methods  and  general  procedures 
for  the  conceptual  design  of  tactical  missiles.   As  mentioned, 
this  is  not  necessarily  the  best  design.   Now  that  the  design 
procedure  has  been  attempted  once,  the  reason  for  this  is  readily 
seen.   The  process  is  one  of  continuous  compromise.   As  was  shown 
in  Chapter  7 ,  the  optimum  wing  for  lift  is  not  the  best  wing 
for  minimizing  drag.   We  also  saw  in  Chapter  8  that  increasing 
chamber  pressure  increases  thrust  at  the  expense  of  increased 
weight.   This  thesis  has  tried  to  point  out  some  of  these  areas 
of  compromise  and  present  methods  to  deal  with  them. 

Areas  which  were  not  covered  which  need  to  be  investigated 
in  conceptual  design  are  structures,  radar  cross-section  and  cost, 
With  the  increased  emphasis  on  survivability  and  the  decreasing 
budget,  it  becomes  increasingly  important  to  define  the  effects 
of  these  areas  on  design  early  in  the  process. 

The  complexity  of  the  design  process  and  the  need  to  obtain 
timely  and  accurate  information  have  made  it  ideally  suited  for 
the  digital  computer.   The  AEROCF  program  used  in  this  thesis 
is  part  of  a  large  scale  computer  program  (MISSYN)  which  con- 
sists of  modules  for  each  section  of  the  design  analysis. 

With  a  good  understanding  of  the  theory  and  methods  used  in 
missile  design,  the  computer  aided  design  program  with  graphics 
capability  gives  the  designer  the  capacity  to  make  intelligent 
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design  interations  and  almost  instantaneously  see  the  effect 
of  the  change  on  all  areas  of  design.   A  limited  example  of 
this  can  be  seen  in  the  use  of  AER01.   A  change  in  the  per- 
formance requirements  for  maneuver  capability  would  require 
a  redesign  of  the  lifting  surfaces.   A  change  in  the  lifting 
surface  design  would  change  the  drag  characteristics  of  the 
missile  and  therefore ,the  propulsion  requirements.   The  AER01 
program  coupled  with  a  similar  propulsion  module  would  allow 
the  designer  to  make  the  changes  and  instantly  see  the  penalty 
or  savings  in  propellant  weight. 

One  pass  at  the  design  has  been  accomplished  in  this  thesis. 
As  was  seen  throughout  the  process,  decisions  in  one  area 
affect  the  design  in  others.   For  this  reason  the  design 
process  becomes  an  iterative  one.   The  final  design  of  the 
first  iteration  is- .the  baseline  missile  for  the  second  itera- 
tion  and  the  design  is  started  again.   By  making  several 
passes  through  the  loop,  the  solution  converges  on  the  final 
design. 
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